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ANALYSE DE LA RADIOACTIVITÉ DES MINÉRAUX FORMÉS 
AU CONTACT D'UNE VEINE GRANITIQUE AVEC UNE SÉRIE 
CALCAIRE 


STANISLAWA ZMYSLOWSKA 
Institut de Physique de l'Université de Varsovie 
(Manuscrit reçu le 3 mars 1954) 


En employant l'émulsion nucléaire Ilford C2, on a étudié le caractère de la réparti- 
tion des substances radioactives dans des échantillons de roche; on a déterminé avec 
quels minéraux elles ont été prélevées et on a calculé leur teneur générale en thorium 
et en uranium. On a relévé que ce sont des inclusions de surface différentes allant de 
quelques u? jusqu’à environ 1 mm? réparties inégalement qui consitituent les centres 
radioactifs. Les substances radioactives se sont concentrées avec la magnetite. La teneur 
générale en thorium et en uranium est du méme ordre de grandeur dans tous les échantil- 
lons étudiés. - 


1. But du Travail et Méthode Employée 


Le but du travail est d'effectuer, par la méthode de l'émulsion nucléaire, l'analyse 
de la radioactivité d'une roche provenant de Silésie et représentant le point de contact 
d'une veine de granit avec du calcaire. 

L'analyse envisagée devait consister dans les opérations suivantes: 

a) établir des corrélations entre les centres radioactifs et les minéraux contenant 
de l'uranium et du thorium à une concentration sensible; j 

b) étudier le caractère de la désintegration de la substance radioactive dans la 
roche donnée. 

c) évaluer la teneur totale en thorium et en uranium. 

J'ai étudié des surfaces polies, de la matiére en poudre versée dans de petits 
moules ainsi que des comprimés de poudre. J'ai examiné chaque échantillon à plusieurs 
reprises avec des temps de pose différents. Un long temps de pose m'a permis d'étudier 
le caractére de la répartition des groupes de trajectoires de particules alpha; un court 
temps de pose m'a permis de tirer des conclusions quantitatives. 

Pour mon travail, je me suis servie de plaques nucléaires Ilford C2, où l'épaisseur 
de l'émulsion était de 50 y. Les plaques ont été directement appliquées sur les échan- 
tillons de granit étudiés. Après exposition, elles ont été développées à l'aide du révélateur 
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D 19 —b, selon le procédé connu (Powell et Occhialini 1947, Beiser 1952). Vu 
l'épaisseur minime de l'émulsion, il etait possible de ne pas tenir compte de la distor- 
tion. Les plaques ont été séchées dans une position horizontale, à l'air et à l'abri de 
poussiére. Aprés les avoir séchées, je les ai appliquées encore une fois les échantil- 
lons étudiés, j'ai dessiné les contours de ceux-ci et je les ai examinées au micro- 
scope. 


2. Le Microscope 


Раі effectué ce travail à l'aide d'un microscope binoculaire Reichert еп me 
servant de l'agrandissement total 90 fois, 270 fois, 405 fois sans immersion, et en 
employant également l'objectif d'immersion donnant un agrandissement de 100 fois. 
Je me suis servie de l'agrandissement de 90 fois lors de l'examen préliminaire de la 
plaque, et lors du relevé de la répartition générale des groupes de trajectoires. Je me suis 
servie des autres agrandissements lors de la prise des mesures. Afin que les mesures 
soient comparables, bien que des agrandissements différents aient été employés, 
jai etudié la méme surface-témoin de l'échantillon en employant divers agrandisse- 
ments, et j'ai noté que le nombre de trajectoires de particules alpha observées avec 
un agrandissement de 405 fois est 4,2% moindre que le nombre obtenu si l'on employe 
l'immersion ; dans le cas où les observations se font avec un agrandissement de 270 fois, 
ce nombre est 17,4% moindre. Par des calculs, j'ai réduit le nombre de trajectoires 
observées à celui que l'on obtiendrait en employant l'immersion. 

On a fait les microphotographies en se servant du méme microscope. Les surfaces 
polies ont été étudiées dans un microscope. Reichert de type MEF muni d'un 
dispositif polarisant. La lumiére polarisée. permet de distinguer les cristaux, qui se 
détachent sur le fond. Les deux types de microscopes Reichert sont muni d'un 
appareil photographique. La platine porte-objet est munie d'une échelle au milli- 
métre avec un vernier permettant de lire la position avec une exactitude de 0,1 mm. 

J'ai étudié la plaque par bandes horizontales. Comme champ dans le microscope, 
jai choisi un rectangle délimité par un oculaire muni d'une grille. Le côté d'un carré 
de cette grille, qui avait été étalonné, permettait de mesurer la surface observée, 
et d'évaluer approximativement la longueur des trajectoires observées. 


3. La Plaque nucléaire 


Les traces enregistrées sur la plaque nucléaire employée sont hétérogènes. On y 
voit des trajectoires de particules alpha isolées, disséminées aussi bien en profondeur 
qu'à la surface de l'émulsion, des étoiles de particules alpha à deux, trois, quatre 
et cinq branches, des groupes de trajectoires de particules alpha, des protons, plus 
rarement des mésons et des étoiles cosmiques. La carte établie selon l'une des plaques 
en est l'illustration. (fig 1). 


Пе l'ensemble de ces données, on choisit uniquement les trajectoires de particules 
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alpha qui commencent sur la surface de l'émulsion (lorsqu'on les a observées sans 
immersion, elles semblent se trouver dans le plan de la surface de l'emulsion), car ce 
sont celles-ci seulement qui peuvent provenir de l'échantillon étudié. La chose est 
rendue plus compliquée, par le fait que l'on peut également observer un certain 
nombre de particules alpha ,,de surface“ sur la plaque non exposée. Il est impossible de 
distinguer ces trajectoires parasites des trajectoires de particules alpha ш sont sorties 
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Fig. 1. Carte de la répartition générale de différentes particules observées dans la plaque nucléaire exposée. 
Les axes des coordonnées représentent les mouvements verticaux et hirozontaux de la platine du micros- 
cope. La ligne en pointillé à gauche représente le bord de la surface polie exposée. Les chiffres expri- 
ment le nombre de trajectoires de particules alpha sur le champ d'un carré dont le côté est de 0,5 
mm; ces trajectoires proviennent de l'échantillon exposée. Les petits traits isolés et les signes p dési- 
gnent les protons, les lettres a indiquent l'emplacement des étoiles cosmiques, les lettres m, les mésons, 
les lettres k, les centres de contamination concentrés en un point; les petites étoiles de contamination 
n'ont pas été indiquées.. М.В. les longueurs des trajectoires de particules alpha de contamination sont 
agrandies environ 10 fois à l'échelle de la carte. 


isolément de la surface de l'échantillon. On a donc decoupé une partie de la plaque 
avant l'exposition et développé cette section de la plaque en méme temps que la 
plaque étudiée. Un examen de la surface-témoin permet d'évaluer approximativement 
le fond formé par les trajectoires de surface. On soustrait ce fond du nombre total de 
particules alpha observées sur la surface de l'émulsion. Le nombre obtenu représente 
le nombre de particules alpha émises par une surface donnée de l'échantillon pendant 
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un temps donné. Sur la base de ce nombre, nous déterminons la concentration des 
substances radioactives contenues dans un échantillon donné. 


Les particules alpha enregistrées dans l'épaisseur de l'émulsion proviennent 
généralement des impuretés qui y sont inévitablement contenues, ou du verre. Elles 
constituent la contamination. Elles peuvent parfois étre nombreuses et se concentrer 
en un seul point; nous en avons des exemples sur la carte de la figure 1. On donne 
sur la figure 2 la photographie d'un tel centre de contami- 
nation. Deux faits prouvent qu'is'agit là d'impuretés: 1) une 
autre plaque, exposée sur le méme échantillon, ne donne pas 
de groupes semblables de trajectoires de particules alpha 
aux mémes endroits; 2) le centre de contamination se trouve 
plutót sous la surface de l'émulsion. 


50м | 4. Esquisse de la Théorie 


A Les surfaces polies, les comprimés sont des sources 
Fig. 2. Centre decontami- |, 3 г 
nation dans. unes plaque Өргіззев. Г. Cure (1946), Н. Yagoda (1949), R. Coppens 
exposée (1947) et S. M. Tamburino (1952) ont fourni les bases thé- 
oriques permettant de tirer des conclusions quantitatives 
de l'étude par la méthode de l'émulsion nucléaire des échantillons de roche consti- 
tuant des sources épaisses. Les formules déduites par divers auteurs peuvent étre 
représentées sous une forme générale identique: 


N 


Ср и: 
CIE) 
су 


N exprime le nombre de trajectoires de particules alpha émises par l'échantillon 
étudié, d'une surfáce de 1 cm? pendant un temps d'une seconde. En pratique, c'est 
le nombre de trajectoires de particules alpha observées sur la surface de l'émulsion, 
d'une longueur plus grande que ou égale à 3 u, divisée par le temps de pose exprimé 
en secondes et par la surface exprimée en cm?. c;,/c, est le rapport de la teneur en 
thorium exprimé en poids par celle en uranium de la roche donnée. 


Le coefficient k est lié au freinage dans l'échantillon étudié des particules alpha 
que celui-ci émet. Les auteurs cités ci-dessus diffèrent pour ce qui est de la manière 
d'évaluer ce freinage, ce qui entraine des différences dans les constantes a et b. 


J'ai pris pour mes calculs la formule donnée par S. M. Tamburino (1952) ou 
Е = 11, et К s'exprime comme ZcS/4, c'est la teneur (exprimée en poids) de la 


roche en élément donné, S est le pouvoir d’arrét d’un atome, A est le poids ato- 
mique. | 


Dans cette courte esquisse, on voit immédiatement que pour déterminer par 
la méthode des plaques nucléaires la teneur en substances radioactives de la source 
épaisse donnée, il faut établir expérimentalement: 


O — а 
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1) le nombre de trajectoires de particules alpha observées dans l'émulsion nuclé- 
aire, dont j'ai parlé ci-dessus, 

2) le rapport de la teneur en thorium à la teneur en uranium dans les échantillons 
étudiés 

3) la composition minéralogique et, partant, le coefficient k’. 


5. Détermination de k’ 


Aprés avoir été exposés, les échantillons ont été réduits en poudre et étudiés 
du point de vue minéralogique par le Dr. H. Makowski, qui m'a donné le pourcentage 
de leur teneur en différents minéraux. Sur la base de cette teneur, j'ai calculé c pour 
chaque élément contenu dans l'échantillon. Le rapport S/A pour chaque élément 
respectivement a été etabli selon la courbe de S. M. Tamburino (1952). La somme 
des produits cS/A donne le coefficient №. J'ai employé quatre échantillons (de 
quatre fragments différents de la roche étudiée). Je les désignerai par Gr5, Gr5-a, 
G 5-b, Gr5-témoin. Le morceau de roche prélevé provient d'une localité oà l'on trouve 
les minerais de fer dans une série métamorphique calcaire et amphibolitique. Il repré- 
sente une veine de pegmatite entrelacée de calcaire. La composition minéralogique 
des échantillons est la suivante: 


Ст5 : 98,5% calcite 
1,3% substance argileuse 
0,1% magnétite 
0,1% chlorites (non déterminées) 


Gr5-a 75,4% feldspath 
19,7% calcite 
4,9% quartz 
_ 0,02% — magnétite 
Gr5-b 
67,5% calcite 
32,5% magnétite 


Gr5-temoin 
54,8% calcite 
43,8% silice ` 
1,4% oxydes de fer, - 
d’où j'ai obtenu les valeurs suivantes pour k’: 
Gr5 - 0,0579 
Gr5-a - 0,0568 
Gr5-b - 0,0574 
Gr5-temoin 0,0576 
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6. Rapport de la teneur en thorium à la teneur en uranium 


En général, pour definir c7,/cy, on emploie la formule déduite par I. Curie (1946), 
qui permet de déterminer ce rapport d'aprés le nombre de trajectoires longues; on 
peut distinguer deux espéces de trajectoires: l'une comprenant des trajectoires plus 
longues que 39 u, dont le nombre est ту, et que Гоп peut certainement attribuer au 
thorium С’, et l'autre comprenant des trajectoires d'une longueur située entre 32 
et 39 u, dont le nombre est N, et qui sont attribueés au radium С’. 


Мт, CTh 


En 1951 Bremner (1951) a émis l'hypothése qu'il serait possibie de déterminer 
ce rapport d’après le nombre d'étoiles à deux, trois, quatre et cing branches. Poole 
et Mathews (1952) font mention de la possibilité de le déterminer d'aprés les projections 
horizontales des trajectoires de particules alpha. Dans mon cas, les tentatives d'appli- 
quer ces deux méthodes ont donné un résultat négatif: aussi bien les nombres d'é- 
toiles à la surface que les projections des trajectoires des longueurs dont il a été question 
sont trop faibles pour qu'ils puissent constituer une base pour calculer la valeur 
Crhlcu- On a également donné (9) (Tournay 1952) une esquisse de la méthode à suivre 
pour determiner ce rapport à l'aide des rayons X. Cette méthode est maintenant en 


د 


cours d'élaboration à notre Institut. 


J'ai déterminé c суу par la méthode de I. Curie, en calculant la longueur des 
trajectoires des particules alpha, sur la base des mesures de leur projection sur la 
surface de la plaque, et des mesures de leur profondeur, et en tenant compte de l'affais- 
sement de l'émulsion. 


Comme on le sait, il existe toute une serie de méthodes pour déterminer l'affais- 
sement de l'émulsion (Vigneron 1949, Rotblat 1949, Ciok). Dans mon cas, j'ai 
appliqué la méthode qui consiste à mesurer l'épaisseur de l'émulsion aprés le déve- 
loppement de la plaque, en adoptant la valeur standard de l'épaisseur de l'émulsion 
avant qu'elle soit développée, telle qu'elle est donnée par la maison produisant les 
plaques. L'affaissement de l'émulsion, selon l'humidité de l'air, varie de 2,34 à 2,8; 
ce qui est en accord avec les résultats des autres mesures. 


J'ai étudié 2162 champs au microscope, d'une superficie totale de 9,61 m2, 
à l'aide d'un agrandissement (avec immersion) de 900 fois. J'ai mesuré la longueur de 
1165 trajectoires de particules alpha. Parmi celles-ci, j'ai établi la présence de 92 
trajectoires longues, et notamment 65 trajectoires RaC’ d'une longueur plus grande 
que ou égale à 32 u, et plus petites que 39 4, et 27 trajectoires ThC’, d'une longueur 
égale à ou plus grande que 39 u. En me basant sur ces données, j jai déduit de la 
formule de Curie que c7,/c, est égal à 2. Jai employé cette valeur plus tard pour 


determiner la concentration globale des corps radioactifs dans les échantillons de 
roches étudiées. 
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7. Les surfaces polies 


J'ai examiné 4 surfaces polies, que je désignerai par Gr5, Gr5-a, Gr5-b et Gr5-c. 
La surface polie Gr5 a été expoxée à trois reprises avec des temps de pose de 230, 


24 et 6 heures. La première plaque a été étudiée avec un agrandissement de 90 fois 
pour établir la répartition 


générale des groupes, comme 
le montre la carte. J'ai obser- 
vé quelques centaines de grou- 
pes de trajectoires de parti- 
cules alpha. J'en ai tiré la 
conclusion que dans cet échant- 
illon, il existe des substances 
radioactives en inclusions ou 
en petites inclusions dans le 
minéral avec lequel elles ont 
eté préleveés. Dans quelques 
groupes, il a été impossible Fig. 3. Centre radioactif provenant de l'échantillon Gi5. Micro- 
de compter. les trajectoires; photographie de la plaque nucléaire exposée pendant 230 heures 
nous en avons un exemple 

sur la figure 3. La plaque de 24 heures a été étudiée avec un agrandissement 
de 405 fois seulement dans la région oü il existait des groupes et des con- 


Fig. 4. Photographie de la өшіре polie Gr5-a; les chiffres de 1 à 5 désignent les fragments de forte radio- 
| activité mentionnés ci-dessus 
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centrations de trajectoires de particules alpha. On a obtenu une carte d'essai. Comme 
base permettant d'évaluer la concentration générale de substances radioactives dans 
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Fig. 5. Carte de la plaque nucléaire Gr5-a. Оп a indiqué les contours de l'échantillon étudié qui ont été 
dessinés pendant l'étude de la plaque nucléaire 


l'échantillon donné, on s'est servi de la plaque de 6 heures, qui a également été étudiée 
avec un agrandissement de 405 fois. La carte de la figure 1 a été établie selon cette 


Fig. 6. Photographie de la surface polie Gr5-b. A cóté de la fissure qui traverse tout l'échantillon, on 
remarque des taches blanches qui montrent la kaolinisation du feldspath 

plaque. Les chiffres sur cette carte donnent les trajectoires de particules alpha obser- 

vées dans les contours de l'échantillon et émises par celui-ci. La surface étudiée est 

de 11,625 cm?, ce qui constitue 42150 champs au microscope. J'ai observé 13005 
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trajectoires de surface de particules alpha. Aprés avoir déduit le fond de trajectoires 
de surface il reste 9880 trajectoires. Si nous tenons compte des pertes dues au fait 
que l'étude a été faite sans immersion, nous obtenons 10290 trajectoires de particules 
alpha d'un parcours plus grand que ou égal à Зи. 

La surface polie Gr5-a (fig.4) a été exposée à deux reprises pendant 134 et pendant 
12 heures respectivement. Un examen de la plaque photographique exposée pendant 
le temps plus long avec un agrandissement de 90 fois a montré des centres radioactifs 
peu nombreux, mais intenses, d'une répartition particuliérement caractéristique. 
La plaque photographique de 12 heures a servi pour faire la carte (fig. 5) et pour déter- 
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Fig. 7. Carte de la plaque nucléaire exposée sur l'échantillon Gr5-b; les groupes de trajectoires des 
particules alpha correspondent au trace des fissures dans l'échantillon étudié 


miner la concentration d'uranium. Je l'ai étudiée avec un agrandissement de 270 
fois. Dans le contour de l'échantillon, sur un champ de 5,29 cm?, j'ai observé 6318 
trajectoires de surface de particules alpha; aprés avoir déduit le fond — 4840 et 
aprés avoir pris en considération les pertes par rapport à l'étude avec immersion — 
5680 trajectoires.: . L^ 
La surface polie Gr5-b (fig.6) a éte exposée à deux reprises, avec un temps de 
pose respectivement de 134 et de 48 heures. Sur la plaque photographique qui a été 
exposée plus longtemps, j'ai pu déceler un arrangement différent de groupes de 
trajectoires de particules alpha, et notamment des chaines entiéres. la carte établie 
selon la plaque photographique exposée pendant 48 heures donne une image de cet 
arrangement de trajectoires de particules alpha (fig. 7). J'ai déterminé la concentra- 
tion d'uranium d'aprés cette plaque. Je l'ai étudiée avec un agrandissement de 270 
` fois. Le- nombre de trajectoires de particules alpha venant de l'échantillon s'éléve à 
11595; aprés avoir déduit le fond et tenu compte des pertes subies par suite de Г exa- 
men fait sans immersion, il reste 12400 trajectoires de particules alpha. 
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La surface polie Gr5-c a été exposée pendant 28 heures et étudiée avec un agran- 
dissement de 90 fois uniquement en ce qui concerne la répartition générale des tra- 
jectoires de particules alpha. Cet échantillon a montré l'arrangement en chaine des 
groupes de trajectoires. 

En établissant la carte des plaques étudiées au microscope, il a fallu effectuer 
les calculs d'aprés les bandes de la plaque observée, afin de pouvoir placer les données 
relevées à une échelle appropriée en millimétres. J'ai fait les cartes à l'échelle 10 : 1. 


Fig.8. Microphotographie du fragment (1) de l'échantillon poli Gr5-a (Voir les figures 3 et 4) 


8. Corrélation entre les centres radioactifs intenses et les fragments correspondants 
des surfaces polies 


Cette correlation se fait à l'aide d'une carte de la plaque photographique et du 
contour d'un échantillon tracé sur la plaque. Là où il s'agit de groupes d’une surface 
minime, méme si les groupes sont nombreux, comme nous le voyons sur la figure 3, 
il est presqu' impossible de faire la corrélation de cette facon; on peut la faire seule- 
ment à l'aide d'une plaque mince transparente. C'est precisément dans le cas de la 
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figure 3 que la corrélation n’a pu être faite. La cause de cet insuccès était la structure 
finement granulée de l'entourage de ce centre dans l'échantillon. Cependant, aussi bien 
dans la surface polie Gr5-a que dans la surface Gr5-b, on a réussi à faire la corrélation. 

L'échantillon Gr5-a donne 5 centres trés intenses, où l’arrangement des tra- 
jectoires de particules alpha est tellement caractéristique qu'il permet d'identifier 
sans aucun doute l'emplacement des minéraux qui correspondent à ces groupes et 
qui ont la méme forme. Sur les figures 4 et 5, ces centres sont numerotés de 1 à 5. 
Comme exemples, je donne la microphotographie de l'endroit (1) de l'echantillon 
faite à la lumiére polarisée (fig. 8) et la répartition des trajectoires de particules 
alpha dans la plaque photographique qui lui correspond, plaque qui a été exposée 
pendant 134 heures (fig. 9). La figure 10 donne un fragment du méme groupe de la 


Fig. 9. Répartition des trajectoires des particules alpha à l'emplacement de la plaque: photographique 
correspondant au fragment (1). La microphotographie montre que la petite inclusion (1) est responsable 
de l'émission des particules alpha (Voir figures 3 et 4) 


plaque photographique avec un temps de pose plus court. Aussi bien le fragment 
d'échantillon representé sur la figure 8 que les autres endroits de 2 à 5 ont été reconnus 
par le Dr. H. Makowski comme étant de petites inclusions de magnétite dans les- 
quelles se sont concentrées les substances radioactives. ; 

Dans l'échantillon Gr5-b, la répartition radioactive est liée de facon frappante, 
qui ne laisse aucun doute, aux petites fissures, се que по pouvons facilement 
remarquer en comparant les figures 5 et 6. On voit que les oxydes de fer, qui contiennent 
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de l'uranium et du thorium, ont rempli les fissures de la roche. Il est à supposer que, 
dans le voisinage des endroits où la radioactivité est la plus intense, le réseau cristallin 
de feldspath a été détruit, et le processus de kaolinisation a commencé. On voit 
clairement des taches blanches dans l'entourage de la fissure (fig. 6). Il est aussi 
possible que la cause de la kaolinisation soit différente. 


9. Poudres et Comprimés 


J'ai étudié deux échantillons de la poudre-étalon Gr5 avec des temps de pose 
de 337 heures et de 50 heures respectivement. Les surfaces étudiées sont respective- 
ment de 1,04 cm? et de 0,67 cm?. Le 
nombre observé de trajectoires de particu- 


les alpha provenant du premier échantillon 
T з, est de 39192 (compte tenu du fonds et 
Xo | des pertes subies lors de l'examen à un 
agrandissement de 270 fois, il s'éléve à 
45750 trajectoires); pour le deuxième 
échantillon, étudié également avec un agran- 
dissement de 270 fois, on a compté 3760 
trajectoires de particules alpha (compte 
tenu des corrections, 4250). 

Le comprimé-étalon Gr5 de poudre, 
dont les grains étaient plus petits que 50 u, 
aprés avoir été exposé pendant 20 heures, 
a été étudié avec un agrandissement de 
270 fois; sur un champ de 1 cm?, j'ai 
observé 1968  trajectoires de particules 
alpha. Un autre champ de 9,61 mm? a été 
ensuite étudié avec un agrandissement de 


Fig. 10. Microphotographie d'une partie du 900 fois, avec immersion, ce qui a donné 
groupe (1) de la plaque nucléaire soumise à un 1165 trajectoires d'une longueur plus 
temps de pose plus court (Voir figure 8) grande que ou égale à 34. 


10. Résumé des Conclusions 


a) Les centres radioactifs correspondant à des fragments convenables de roche 
permettent de tirer la conclusion que l'uranium et le thorium se concentrent avec 
la magnétite dans les échantillons étudiés. 

b) Le caractère de la répartition des substances radioactives dans les échantillons 
étudiés est le suivant: les substances radioactives entrent dans la composition des 
inclusions microscopiques qui forment des petites taches noires et dans la composition 
des veinules de magnétite remplissant les fissures de la roche. 
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c) Teneur en thorium et en uranium: les résultats de la statistique, en prenant 
pour k’ et cz,/cy les valeurs données ci-dessus permettent de déterminer la teneur 
moyenne en uranium, et ensuite, en thorium. Elle est du méme ordre pour tous les 
échantillons donnés. 

En ce qui concerne l'évaluation des erreurs, j'ai estimé l'erreur totale commise 
lors du calcul de la moyenne en uranium à mois de 10%. 
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КРАТКОЕ СОДЕРЖАНИЕ 

Змысловская, Анализ радиоактивности минеральных творений на кон- 
такте пранитной жилы с кальцитной серией. 

Используя ядерную эмульсию Илфорд С 2 исследовано характер распределения 
радиоактивных веществ в образцах скалы; установлено, с какими минералами 
они собираются, и определено общее содержание тория и урана. Установлено, 
что радиоактивные центры образованы врастающими областями с различной 
поверхностю, от нескольких и? до 1 мм?, появляющимися неравномерно. Радио- 
активные вещества накоплены вместе с магнетитом. Общее содержание тория 
и урана имеет такой же порядок величины для всех исследованных образцов. 
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THE INFLUENCE OF HINDERED ROTATION ON THE 
SCATTERING OF SLOW NEUTRONS BY BOUND PROTONS II* 


By WropnziwrERZ Koros 
Institute of Physics, Polish Academy of Sciences, Warsaw 
(Received June 21, 1954) 


The influence of rotationo f a rigid molecule on the cross-section of a proton bound 
in the molecule is considered. The molecule is assumed to be a symmetric rotator. Using 
the »static approximation« of Placzek, contributions to the cross-section due to various 
changes of the rotation about the axis of symmetry are calculated. !t was shown in Part I 
that this enables us to calculate the influence of hindered rotation on the proton cross- 
section for collisions elastic with respect to the hindered rotation. 

The method is applied to the methyl alcohol molecule. It is shown that the influence 
of hindered rotation makes the cross-section of the molecule about 13 per, cent smaller 
compared to that of the molecule assumed to be rigid. A rough estimate shows that the 
remaining experimentally observed effect may be due to vibrations of the protons. 


1. Introduction 


It has been shown in Part I of this work! (Kolos 1954) that the cross-section of 
a proton bound in a molecule having an additional degree of torsional freedom is 
given by : 
бор тұр; = bo Omg, mi? (1) 
where m denotes the three rotational quantum numbers (j, T’, т), о is the quantum 
number of the hindered rotation, оњу, т; is the cross-section of à proton bound in 
a rigid rotating molecule associated with the transition from the initial state m; to 
the final'state m,, and b, is given by 
oo 
b, = У а, (05, тда, (67у), Q) 
n=0 
where a, (0,т) are the coefficients of the Fourier expansion of the periodic part of 
the Mathieu function describing the hindered rotation. 


* Supported by the Polish Academy of Sciences. 
1 Henceforth designated by I. | 
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For collisions in which neither the hindered rotation nor the rotation of the 
whole molecule about the axis of the hindered rotation is excited, we have b, = 1, 
and for collisions in which the latter is excited, we have by < 1. 

To get a result which may be compared with experiment, we have to calculate 


o= у Py У b} он, (3) 
2 у 


where Ри; is the probability of the molecule to be initially in the state m;. The summa- 
tion over m, is only limited by the conservation of energy, and as the differences 
between the rotational levels are small compared to the neutron energy we have 
to calculate a very great number of contributions associated with all possible transitions. 
This makes the method impracticable in general. However, in the next section we 
shall show that in the case when the neutron energy is large compared to the mean 
level spacing of the initial rotational states of the molecule, we may greatly simplify 


the calculations. 
2. The influence of rotation 


In Fermi's approximation, i.e. assuming a Ó-type neutron — proton interaction 
and using Born's first approximation, the differential cross-section of the proton being 
initially in the state 7 is given by the well known poss 


> f 9 exp (ikr) P} (с) de, | 40, (4) 


where f denotes the final state of the proton, Ÿ. Ar, p) and Wilts ) are the wave functions 


of the proton in _the initial and final states respectively, r, is the position vector of 


the proton, k= p Ee Р) ћ where Р is the momentum transferred to the proton. 
If the Shear energy of the scattering medium is large compared with the 
rotational constant of the molecule, the initial rotational wave function can be replaced 
by a rotating wave packet with given orientation and angular velocity at the beginning 
of the collision. If in addition the energy of the incident neutrons is large compared 
with the mean level spacing, the wave function of the incident neutrons can also 
beapproximated by a wave packet of such a small extension that the time of collision 
is small compared with the period of rotation. In this case the unique role of the 
matrix elements on the right hand side of (4) is to guarantee the conservation of linear 
and angular momentum i.e. to require that the energy transferred to the translational — 
rotational motion of the molecule has the same value as in a classical collision. This 
latter according to Sachs and Teller (1941) is given by 
Е= Ў ЪЪ РМР (5) 
7 s 2 ? 
where V, is the initial velocity of the scattering nucleus and M is the mass-tensor 
introduced by Sachs and Teller. 
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Thus the transition probability for the neutron to be scattered into the solid 
angle dQ will be proportional to 


dQ f 8 (E, — E; + E) P?aP,, 
0 


where E, is the initial energy of the neutron, Е and P ‚are the energy and momentum 
of the scattered neutron. 

It was shown by Messiah (1951) that in this way using the method of Placzek 
(1952) and averaging the cross-section over random orientation of the molecule 
and over the initial values of Рр. we get the same result as in the classical theory of 


Sachs and Teller (1941), namely 


y ЕТ 
о = 46; (Up Mes)" + b TC (6) 
where u;s(i = 1,2,3) are the principal values of the reduced mass tensor and 
3 
I 
= 3 D Ш; - 


Or, neglecting the thermal energy of the scattering molecules (which is equivalent 


to puting Г, == 0 in the preceding steps), we get 
б = 40% (ut Mots) ^. | (7) 


However, we have to calculate the contributions to the cross-section due to 
various transitions associated with the rotation of the molecule about the axis of 
the internal rotation. So we must go back to formula (4) and calculate first the pro- 
portionality constant c. 

When thermal neutrons are scattered by protons bound in heavy molecules, 
then, for low temperatures of the scattering medium, the neutron energy is sufficiently 
large to excite high rotational transitions. But it may be shown by numerical calcula- 
tions, (see also Placzek 1952), that the transition probabilities decrease rapidly for 
higher transitions and the energy transferred to the protons will be small for all 
transitions with appreciable matrix elements. So we may neglect approximately 
the changes of the neutron energy associated with the transitions of the system, 
and the summation in (4) may be carried out by clusure. This is the static approxima- 
tion of Placzek, which in our case gives for the total cross section 


с == Але (8) 


Comparing (8) with (7), we get 


[07 
с = — up papa). (9) 


We see that the callculation of the influence of rotation requires a generalization 
of the reduced mass factor u? in (2,1), which must be replaced by (uj us us) ^ 
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Therefore, using (4) with (9), taking for (r „) the eigenfunctions of the symmetric 
rotator, and choosing a coordinate system with z-axis in the direction din k, — kes 
we get for the differential cross section 


dom, = 29 (ишигиз) * У 
л 


ІШ .vee [ik(c cos 0 


n, 
+ a sin 0 sin ¢)| Umf* (0, v,9) sin 0d0dydg d2, (10) 
where 
U, (0, Y, — @ (бейе v vo) | 
( Virorum Wm m(0) 


6,0-4 (1-1)? F(n;1) 


1 


e 1 p!(d!)? (s +p)! 
„= | O — Pn: dt = 1 + d +s + 2р ‘(d+s+p)! )! (d + p)! (11) 
0 
t= 5 (l — cos 6) 


d= vett; Sese ser 
F(m;t) = F(—p, 1 +d +s +p, 1 + d; ©, (р=0,1,...) 
Using now the static approximation, we shall neglect the changes of the neutron 
energy due to transitions of the system by puting 
k = kky, (12) 


where x is the unit vector in the direction of scattering, and we shall replace the 
restricted sum over j and v, by an unrestricted one. Carring out the summation 
оуег:) and т by closure, we get 


_ 04123) * D ; 
dom; = Ber - 2 dQ [ e exp (Кт; = т/ф 

— i(7; — 1j)g' + Ша sin 0 sin p — Ша sin Ө sin y } sin б40фрарар. (13) 
Integrating in (13) over y, y and ф’, we get 


о Е Р | 
ы cuu > 40 ii Om? (ka sin 6) sin 828, ca 


where, for abbreviation, we have written т, = T; — Ty, and Л, (ka sin 0) is the 
usual Bessel function of order Tip. 


Expanding the square of the Bessel function in a power series by the aid of the 
formula 
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© (—1)" (2k + 2n)! Be 


(15) 
I) = п! (Zk + n)! [(k + n)!]? 
and writing 
2p 
2 (тт) = У ea (16) 
the integration over £ in (14) gives for the differential cross section 
y кісе E 3 Ves nw + 2n) 1( + vi, + n) ka)?rig t? 
| buds nel, F TT 
2p (d 3- m 4- т; ae n)! 3 (17) 
га (4 +s + m +27, 2n + 1)! ок 
and for the total cross section 
gu’ E.a2\ tit” 
„ ger, + 216 + r, + (E) 
a 400 (H pa Hapa) ۶ iege AT h 
"GE n! (2т;, + п)! (ту, + л) 1] 
UF n=o 
(d +m + tytn)! ud 
á Dp (HOME 21, H2n + 1)!” 


Formula (18) may be written in a form more suitable for numerical calculations, 


namely 
409 (шилизиз)”* tat eur ДЕШ; 
nr) `> Y erra n! (Ory +m)!’ (19) 
where 
8u’ E; 
Gy) m = азои («+ m 4 n)! 
en (n!) NEES ETE ROSEN 


т=о 


3. Comparison with experiment 


Now we want to apply our results to the calculation of the cross-section of the 
methyl alcohol molecule. The molecule will be taken to be approximately a symmetric 
rotator, the axis C-O being the axis of symmetry. Using recently published spectroscopic 
data Ivash and Dennison (1953) have found the height of the potential barrier in 
the methyl alcohol molécule to be 384 cm~. 
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The coefficients 62 for protons bound in the CH; group may be calculated by 
means of formula (2). As it was shown in I, for the proton bound in OH the coefficients 
b, are given by 

b, = У а, (0p Ti) ан-т,, (оу, ту), (21) 
п=0 

But Koehler and Dennison (1940) have shown that for a given quantum number 
of internal rotation the coefficients of the Fourier expansion of the periodic part 
of the Mathieu function satisfy the following relation 


Gy, — kh (т) = dn (= 35 3 , (22) 
В, 
where 


4; 


Mn TE 


(7 = 50) (23) 
and A,, A, are respectively the moments of inertia of the ОН or СН, group in respect to 
the axis of hindered rotation. Using (22), we get from (21) for the proton bound in OH 


oo 
NEE = 
2 À A 

We have not indicated in the formulas given above the dependence of a, on 
the quantum number @ because we are interested only in elastic scattering with 
022202220 ; 

The numerical values of the coefficients b? for some transitions are given in 
Table 1. They are calculated on the assumption that the height of the potential 
barrier is 384 cm”! and are average values for the three states which arise owing 
to the threefold symmetry of the hindering potential. (In the table of the coefficients 5? 
published in Part I there are some numerical errors) 


0-1 0++3 0++5 0-10 | O++1 и 


For the proton bound іп the CH, group 


Table 1 


Transition 


Tj ++ Tf 


For the proton bound in OH 


Considering only the influence of rotation of a rigid molecule, we obtain for 
the proton cross-section 40 (u и; из U3) ^, but using formiula (19) we are able to calculate 
the contributions to that cross-section due to various changes of the quantum number r. 
The results of the numerical calculations for two initial states Орта OR = 0 
and j = 15, т” = 15, т = 0 are shown in Fig. 1. 
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We must remember that in our approximation the proton cross-section depends 
only on the absolute value [т 


‚ во that the proton cross section is given by 


On; = Om, (т, = 0) +2 V o, (т, = n). (25) 


Taking into account the influence of the hindered rotation, we have to multiply 
every term in (25) by a corresponding coefficient DA | 


08 
06 
04 


04 


02 02 


D үш = КА з Eile) 2m 002 003 Ele) 
a) initial state j=0, т'=0,т=0 b) initial state j=45, т'=15 т=0 
Fig. 1. The contributions to the total cross-section due to various changes of the quantum number T 


Curve 0 gives the elastic cross section (r;, — 0), curve p (p — 1, 2,...) gives 


P 
Om; = Om (Cy 0) + 2 33 Om; (Tig = п) 


n= 


The resulting cross sections of the methyl alcoħol molecule as a function of 
the initial neutron energy, are shown for some initial states, in Fig. 2. The great 
. differences between the cross-sections for the initial states such as e.g. / = 2, т = 0, 
t = 2 and j = 2, т = 0, т = 0 are due to various orientations of the molecule with 
respect to the initial momentum of the neutrons. 

To calculate a cross-section which may be compared with that measured experimen- 
tally we must apply the formula (3). The distribution of the molecules over the initial 
rotational states depends on the temperature T' of the scattering medium according 
to Boltzmann's statistics. Averaging the cross section over the initial rotational states, 
we get the final result which is shown in Fig. 3. and compared with experimental 
results of Janik (1954). In the calculations a Maxwell-like distribution is assumed 
for neutrons with an effective energy of $ kT. The calculated cross section is only 
an elastic one with respect to the hindered rotation. | 
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Fig. 2. The cross-section of the methyl alcohol 
molecule in units of that of the rigid molecule as 
a function of the initial neutron energy, for some 
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Fig. 3. The cross-section of the methyl alcohol 
molecule averaged over the initial rotational states. 
The vertical lines are the experimental results 
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We see from Fig. 3. that the agreement with experiment is not satisfactory. It must 
be remembered that the experiments of Janik on methyl alcohol were performed at 
room temperature, and the static approximation used in this paper holds better for lower 
temperatures of the scattering medium. But the principal reason of the discrepancy 
between the theoretical and experimental results may be found in our simplified 
model in which the vibrations of the protons have been neglected. The lowest vibra- 
tional frequency of the protons in the CH, group of the methyl alcohol molecule is 
1034 cm! (Volkenshtein, Eliashevitch, and Stepanov 1949). For the methane molecule 
the lowest vibrational frequency of the protons is 1304 cm ! (Volkenshtein et al. 
1949) and for that molecule Messiah (1951) has calculated the influence of vibration, 
rotation, and translational motion on the cross-section. Neglecting the influence of 
translational motion his results in the interesting for us interval of energy are 
shown in Fig. 4. 

If we now assume the influence of vibrations on the cross section of the methyl 
alcohol molecule to be the same as for a methane molecule, we get a better agreement 
with experiment, as may be seen from Fig. 5. 

However, the agreement with experiment must be taken with some reserve 
on account of the rough estimation of the influence of vibrations, and of our 
simplified model in which the influence of translational motion and the association 
of the molecules in the liquid state have been neglected. 


I would like to express my gratitude to Professor L. Infeld and to Doc. J. Werle 
for helpful and valuable discussions. 


КРАТКОЕ СОДЕРЖАНИЕ 


В. Колос, Влияние заторможенного бращения на рассеяние медленных ней- 
троноб сбязанными протонами П. 


Рассмотрено влияние вращения жесткой молекулы на эффективное сечение 
рассеяния медленных нейтронов протонами, связанными в этой молекуле. Пред- 
полагается, что молекула является симметричным ротатором. Используя ,,/стати- 
стическое приближние” Плачека вичислено вклады в эффективное сечение, про- 
исходящие от различных изменений вращения вокруг оси симметрии. 3 части 
Т показано, что это дает возможность вичисления влияния заторможенного вра- 
щения на протонное эффективное сечение для столкновений упругих по отно- 
. шению к заторможенному вращению. 

Метод применяется к молекуле метилового спирта. Показывается, что влияние 
заторможенного вращения дает эффективное сечение молекулы на 13 процент 
меньшим, чем в случае жесткой молекулы. Приближенная оценка показывает, 
что остаюшаясь часть экспериментально наблюдаемого эффекта может происхо- 
дить от вибраций протонов. 
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ÜBER EINE VERALLGEMEINERUNG DES 
REZIPROZITATSTHEOREMS FÜR LOSUNGEN DER 
SCHWINGUNGSGLEICHUNG MIT MULTIPOLQUELLEN 


Von А. RUBINOWICZ 
Staatliches Mathematisches Institut, Warschau 
( Eingegangen am 28. Juni 1954) 


Der Satz von der Vertauschbarkeit von Quell- und Aufpunkt einer Greenschen 
Funktion der Schwingungsgleichung Au + k? u = 0 wird auf den Fall zweier Lösungen 
ul) und u(2) dieser Differentialgleichung verallgemeinert. Dabei wird vorausgesetzt, 
daß jede der beiden Funktionen u(U bzw. u(2) nur eine einzige in dem Raumpunkte 1 
bzw. 2 gelegene singuläre Stelle aufweist, die durch beliebige Multipolquellen ausdrückbar 
ist. Die angegebene Verallgemeinerung des Reziprozitätstheorems stellt dann eine bilineare 
Beziehung zwischen den nachstehenden beiden Sätzen von Koeffizienten dar: 

1) den Koeffizienten, die die Stárken und Phasen der Multipole der beiden Funk- 
tionen u(Dund u(2 in ihren Quellpunkten 1 bzw. 2 festlegen und 

2) den Entwicklungskoeffizienten der nach Kugel- und Besselschen Funktionen 
fortschreitenden Entwicklungen von u(t) und u(2 in den Raumpunkten 2 bzw. 1, in 
denen diese Funktionen regulär sind. 


In der angewandten Physik z.B. in der Festigkeitslehre bietet bei der Behandlung 
spezieller Probleme die Verwendung eines Reziprozitätstheorems oft erhebliche 
Vorteile. Für gewisse Fragestellungen dürfte das gleiche auch in der angewandten 
Akustik gelten.. Bei Vorhandensein einer einfachen kugelsymmetrischen Quelle 
wird das in Rede stehende Theorem mit dem Satze von der Vertauschbarkeit von 
Quell- und Aufpunkt in der Greenschen Funktion identisch. In der angewandten 
Akustik scheint aber dieses spezielle Reziprozitütstheorem praktisch wertlos zu 
sein, da das Strahlungsfeld einer technischen Schallquelle nicht einer einfachen 
Quelle entspricht, sondern der Superposition einer Gesamtheit von Multipolquellen. 
Im folgenden wird daher die Frage untersucht, wie das и für 
solch einen allgemeinen Fall zu formulieren ist. 

81. Die Lósung der Schwingungsgleichung. Wir betrachten eine stationäre 
mit einer gegebenen Frequenz erfolgende Schallausbreitung in einem sich eventuell 
auch ins Unendliche erstreckenden Raume R, der im Endlichen durch irgend welche 
Flächen F begrenzt wird. Die Schallausbreitung werde dabei durch eine im Punkte Q 
befindliche Schallquelle vom gegebenen Multipolcharakter bestimmt. Falls der 


(183) 
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Raum R endlich ist, soll vorausgesetzt werden, daß die der Konstanten k entsprechende 
Schwingungsfrequenz mit keiner der diskreten Eigenfrequenzen des Raumes R 
zusammenfallt. 

Das mathematische Problem, dem wir in diesem Falle gegenüberstehen, besteht 
dann in der Aufsuchung einer mit Ausnahme des Quellpunktes Q überall im Raume R 
regulären Lösung u (x, у, z) der Schwingungsgleichung Au + k?u = 0, die im Punkte Q 
eine vorgeschriebene Singularitát besitzt. In den Flüchen F soll sie ferner eine der 
beiden Randbedingungen 


= о (1) 
on 


und schlie?lich — falls sich der Raum R ins Unendliche erstreckt — im Unendlichen 
die Sommerfeldsche Endlichkeits- und Ausstrahlungsbedingung (Sommerfeld 1913, 
Rellich 1943) erfüllen. Diese letztere Forderung schlieft für einen unendlichen 
Raum R alle Eigenschwingungen aus. Unter den gemachten Voraussetzungen ist 
somit die Lösung u für jeden Raum R eindeutig bestimmt, sei es daß er endlich 
oder unendlich ist. 

Um nun das Verhalten von u im Quellpunkte Q näher zu charakterisieren, 
bemerken wir daß Funktionen, die Multipolquellen im unendlichen Raume darstellen, 
erhalten werden, wenn die Schwingungssgleichung durch Separation der Variabeln in 
räumlichen Polarkoordinaten г, 9, y gelöst wird. Setzt man о = kr und bezeichnet mit 


Yr(d,g) = е” Pr (cos 8) (9) 
eine nicht normierte Kugelflächenfunktion so wird im unbegrenzten, unendlichen 
Raume eine punktfórmige Multipolquelle durch 

Ү"(9, $) 6„(0) (3) 
0 


gegeben. Hier ist 


л 
o = (2) 42,00. | e) 
wo Ho die Hankelsche Funktion bedeutet. Da £,(e) asymptotisch durch 
GC) enr езй 
gegeben wird, so stellt (4) eine vom Koordinatenursprung aus sich fortpflanzende 
Welle dar, falls der Zeitfaktor durch е! gegeben wird. 


Im folgenden nehmen wir, daß die in Q vorgeschriebene Singularität in der 


Umgebung dieses Punktes durch eine nach (3) forstchreitende Reihe charakterisiert 
wird. 


Für die durchzuführenden Überlegungen ist es notwendig die Definition der 
Kugelfunktionen (2) auch für negative m festzulegen. Wir setzen daher voraus, daß 


P™(cos 0) = Р!"\(соз 9) 
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ist. Unter dieser Annahme gilt dann die Relation 


S Yz (9.9) YZ (дур) dQ = М, б (5) 


nym n,n' om a 


wo d ein Raumwinkelelement bezeichnet. Die Integration ist über die ganze Ein- 
heitskugel zu erstrecken. N, stellt das Quadrat des Normierungsfaktors dar. 

Wir legen uns nun die Frage vor, wie die Funktion u in der Umgebung eines 
im Endlichen liegenden Raumpunktes entwickelt werden kann. Eine Lösung, die 
die gleichen Symmetrieeigenschaften um den Koordinatenursprung aufweist wie (3), 
die aber im Gegensatz zu dieser Funktion im Koordinatenursprung endlich bleibt, 
wird durch 


Yr(9,9) ¥ 
0 


gegeben. Dabei wird 


2 
| В 
5 1% fi oo 2 
Palo) = (re) Ins wl) TE @ +) || er 
. PELO NA 


UEBER 


2 2 
durch die Besselsche Funktion Jn+% in der gleichen Weise dargestellt, wie sich C,(o) 
durch die Hankelsche Funktion Н@ ar ausdriickt. 


Bezeichnen wir nun die Entwicklung der Funktion u mit dem Quellpunkt in Q 


іп der Umgebung eines regulären Punktes P mit и(О, P) und in der Umgebung ihres 
Quellpunktes Q mit u(Q, Q) so wird: 


u(Q,P) — xS > Ym 39) a(Q,P Fre) 0, - (6) 
n= 0 m= —п 
u(Q, Q) = Y Ps Үт(9,ө) | AGO) 9 + a20) sep (2) 


wo die kleinen und die groflen Buchstaben a und 4 Entwicklungskoeffizienten 
bedeuten, die im allgemeinen komplexe Werte annehmen. 

Die Entwicklung (6) um einen vom Quellpunkt Q verschiedenen Raumpunkt P 
besteht nur aus einem regulären durch die Koeffizienten а?) charakterisierten Anteil. 
Die Entwicklung (7) um den Quellpunkt Q enthàlt hingegen auf'er dem regulüren 
durch die Koeffizienten a( 2) festgelegten Anteil auch einen singulären Anteil, der 
durch die Koeffizienten АО) bestimmt wird. Die Koffizienten 409 legen dabei 
die Intensitäten und Dis der einzelnen Multipolquellen fest. Durch ihre Angabe 
wird die Funktion u im Raumgebiete R eindeutig bestimmt. 
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Es sei noch darauf hingewiesen, daß die Achsen der Koordinatensysteme in 
den beiden Punkten Q und P keineswegs zu einander parallel sein müssen. 
Schließlich bemerken wir noch, daß wenn das Raumgebiet R ins Unendliche 
reicht die Funktion u hier durch die Reihe 
co 
u(Q, = X у у"(9 ә) 49: оо) 2107 e) (8) 
п-0т---п 
dargestellt werden kann. Die Glieder mit %,(0) werden ja durch Фе Sommerfeldsche 
Endlichkeits- und Ausstrahlungsbedingung ausgeschlossen, da sie ja nicht nur ins 
Unendliche fortschreitende sondern auch aus dem Unendlichen kommende Wellen 
enthalten. i 
82. Herleitung des verallgemeinerten Reziprozitätstheorems. Es seien im 
Raumgebiete R zwei Lösungen der Schwingungsgleichung u und и gegeben, 
deren Quellpunkte in den Raumpunkten 1 bzw. 2 gelegen sind, und die sonst alle die 
in 81 angegebenen Eigenschaften besitzen. Aus dem Greenschen Satze folgert man 
dann, daß 
du?) ә OU) 


n | df = 0, (9) 


Koo +F+Kı +K: 
wo das Flãchenintegral über die Flãchen F, ferner über zwei kleine Kugelflichen 
K, und K, um die Quellepunkte 1 und 2 und schliefilich — falls das Raumgebiet R 
ins Unendliche reicht — auch über eine unendlich ferne Kugelfläche К. zu erstrecken 
ist. Dabei bedeutet n die innere Normale des von allen diesen Flächen eingeschlossenen 
Raumes. 

Wir zeigen zunächst, daß das Flächenintegral (9) über die eventuell vorhandene 
Fläche K verschwindet. Hier kann ja sowohl für uf? als auch für u® die Entwicklung 
u(Q, oo) (8) verwendet werden, so daß wir erhalten: 


со 


Da die Differentiation nach о an den Kugelfunktionen nicht angreift, so erhalten wir 
dafür mit Hilfe von mL und no den Ausdruck: 


ЕР 152 s > Аа, D Ard Nn 4 E (=) А. 


п=0 m=—n п =0 т =—n' 0 
бе d [^ a} 
o de > 


der ersichtlich verschwindet. 


pim - 
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Ebenso verschwindet mit Rücksicht auf die Randbedingungen (1) das Integral (9), 
wenn wir es über die Flüchen F erstrecken. 

Es müssen also nur noch die Flüchenintegrale über die beiden Kugelflüchen К; 
und К, berechnet werden. Um zunächst das Integral über Кү auszuwerten, bemerken 
wir, daB in der Umgebung des Punktes 1 für и und и die Entwicklungen u(l, 1) 
bzw. u(2, l) zu verwenden sind. Es ist also: 


du) du) 
(1) же! (2) ТСХ 
Г. дп + on | JE 
к, 
10 
-+ f [oan _ г u(2,1) Zur dQ. ч 
© 
К, 


Mit Hilfe von (5) erhalten wir dafür 


Y z 2 y, d (2) 
Nam AGD а 5227 fes 
E Y E of z о do \e/)’ 


n=0 т---п 


da ein analoger Ausdruck, in dem ausschließlich nur die ф, (о) auftreten, identisch 
verschwindet. Mit Rücksicht auf die bekannte Relation 


6, (0) v. (e) — v. (0)5,(0) = — 
ergibt sich daher für das Integral (10) über die Kugelfläche K, der Ausdruck 


© +n 
l 
k 


n=0 m=—n 


(1,1) 40210 
№, ALY aD, (11) 


Bei der Auswertung des Integrales über K, ist u® = u(1, 2) und u® = u(2, 2) 
zu setzen. Auf dem gleichen Wege erhält man dann 


Yu) Jud) 
). 
Г. (pee == — цд) = у 3 NS AQ ад (12) 


n=0 m=—n 


K: 


Aus (9), (11) und (12) ergibt sich schließlich die gesuchte, Verallgemeinerung 
des Reziprozitätstheorems 


3 у Nom AGD, =) у М AQ, a, (13) 


n=0 т---п n=0 т---п 
8 3. Diskussion. Gleichung (13) stellt eine Beziehung zwischen den m 
koeffizienten Pe Pau und 4552 dar, die die Singularitáten der Lösungen u und u? 
in ihren Quellpunkten bestimmen, und den Entwicklungskoeffizienten а0 und 2059) 
dieser Funktionen in den Raumpunkten 2 bezw. 1, in denen sich diese "Funktionen 
regulär verhalten. Die Entwicklungskoeffizienten a’ und a(22, die das Verhalten 
des regulären Anteils der Funktionen u“ und u in ihren Quellpunkten festlegen, 


spielen in dem verallgemeinerten Reziprozitätstheorem keine Rolle. 
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Auf jeder Seite der Beziehung (13) treten diejenigen Koeffizienten а und A 
auf, die sich auf die Entwicklungen in ein und demselben Raumpunkt beziehen. 
Dabei ist z.B. auf der linken Seite der Gleichung (13) der mit dem Index n,m ver- 
sehene Koeffizient a2”, der zur Entwicklung der Funktion u? in dem Raum- 
punkte 1 gehört, mit dem mit dem Index n,—m versehenen Koefhzienten AGD, multi- 
pliziert, der das singuläre Verhalten der Funktion u“ in dem gleichen Raumpunkte 
festlegt. Die Analogie mit dem Maxwell-Bettischen Reziprozitätstheorem der Elasti- 
zitätstheorie für Kräfte, die an einzelnen Punkte angreifen, ist augenscheinlich. 

Wahlen wir die Funktion и nun speziell so, daß in ihrem Quellpunkte 2 
nur eine einzige Multipolquelle vorhanden ist etwa die, die den Indizes n,,—my ent- 
spricht. Es ist dann nur АСЕ + 0, während alle anderen 422 verschwinden. 
Das Reziprozititstheorem (13) lautet dann insbesondere 


2,2) (152) — 1,1) 2,1 
Nd AC Anom D 3 N, ad ык) а (14) 


п=0 m=—n 


Setzen wir weiter voraus, daß die Funktion u‘ bekannt ist. Der Koeffizient Ж 2) 


о-- то? 


der ihr singuläres Verhalten im Quellpunkte 2 bestimmt, und die Koeffizienten “абы D 


der Entwicklung von и) in dem Raumpunkte 1 sind dann ebenfalls als hel 
anzusehen. Aus (14) kann daher aX d.h. der die Indizes ny, my aufweisende Koeffizient 
in der Entwicklung von u im Ре 2 berechnet werden. Die die Intensitäten und 
Phasen der Multipolquellen im Punkte 1 bestimmenden Koeffizienten A? sind 
ja gegeben. Indem man n, m, alle zulässigen Werte durchlaufen läßt (ny = 1, 2, . 
ee no), kann man so die Entwicklung der Funktion и 

will 


mit den im Punkte 1 
ürlich vorgeschriebenen Multipolquellen angeben. Zur Feststellung des Funk- 
tionsverlaufs d.h. der Entwicklung von u“ in der Umgebung des Punktes 2 muß 
man somit eine unendliche Folge von Funktionen и?) mit je einer einzigen durch die 
Indizes ny, — то charakterisierten Multipolquelle verwenden. Bei diesen Überlegungen 
spielt die Funktion и) die Rolle einer verallgemeinerten Greenschen Funktion. 

Eine noch einfachere Gestalt nimmt das Reziprozitütstheorem (13) an, falls 
auch in der Funktion u‘? nur eine einzige, etwa durch die Indizes n,,—m, charakteri- 
sierte Multipolquelle auftritt. In diesem Falle wird. 


(2,2 12 _ 1,1 1 
N; ANS Te a OM "| pits Hore en t (15) 
Verwendet man normierte Kugelfunktionen (М, „= №, „ = 1) und setzt 
voraus, daß die beiden Koeffizienten А0) und 402) die die Intensitäten und 


путі пото? 


Phasen der Multipolquellen der beiden Funktionen u“ und u® bestimmen, einander 
gleich sind (trotzdem die beiden Multipolquellen von verschiedener Ordnung sind) 
so geht (15) über in die Beziehung 
5 | айыы, = аб (16) 
Dabei Mio 
aum den Entwicklungskoeffizienten der Funktion u® (mit dem Quellpunkt 
im Raumpunkte 1) in dem Raumpunkt 2, dessen Indizes ту, m, durch die Indizes 


= «уе 
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по, — то der Multipolquelle der Funktion u‘? (mit dem Quellpunkt 2) gegeben werden 
und 

OMA den Entwicklungskoeffizienten der Funktion u(? in dem Raumpunkte 1, 
dessen Indizes пу, m, durch die Indizes ту, —m, der Multipolquelle der Funktion и 
bestimmt sind. 

Handelt es sich darum den Funktionswert der Funktion u® im Punkte 2 an- 
zugeben, so reicht dazu die Kenntnis von aq) hin. Zu diesem Zweck genügt somit, 
wie kompliziert auch die Quellen im Punkte 1 sein mögen, nur eine einzige Greensche 
Funktion и), mit einer einfachen kugelsymmetrischen Quelle im Raumpunkte 2 
In diesem Falle geht ja (14) über in 


œo Hn 
№ AGP ai = У, M, NATO, af. 
n=0 m=—n 
Welche Entwicklungskoeffizienten a» von der Greenschen Funktion и in dem 
Raumpunkte 1 dabei gebraucht werden, hängt von der Art der in dem Punkte 1 
bei der Funktion u® vorgeschriebenen Multipolquellen ab. 

Ist auch in dem Raumpunkte 1 nur eine kugelsymmetrische Quelle vorhanden, 

so geht (16) in die Relation 
ayy = а 
über, die nichts anderes bedeutet, als das gewóhnliche Reziprozitätstheorem von der 
Vertauschbarkeit von Quell- und Aufpunkt. 

Die Beziehung (16) kann als eine direkte Verallgemeinerung dieses Reziprozitäts- 
theorems für den Fall angesehen werden, wo die Quellpunkte der beiden verall- 
gemeinerten Greenschen Funktionen durch einzelne Multipolquellen verschiedener 
Ordnungen gegeben werden. 

Der Fall, wo jede der beiden im verallgemeinerten Reziprozitätstheorem ins 
Spiel tretenden Funktionen in mehreren Raumpunkten Multipolquellen besitzt, 
kann mit Rücksicht auf den linearen Charakter der Schwingungsleichung unmittelbar 
auf den hier behandelten Fall zurückgeführt werden. Er stellt somit nur eine triviale 
Verallgemeinerung dieses Falles dar. 


КРАТКОЕ СОДЕРЖАНИЕ 


В. Рубинович, Об обобщении теоремы взаимности для решений уравнения 
колебаний Aut ёи = О с мультипольными истоками 


Теорема о возможности перестановки точки истока и. точки, в которой мы 
вычисляем значение функции Грина уравнения колебаний Aut+k2u=0, обобщена 
для двух решений этого уравнения ul) и u(2), При этом предполагается, что 
каждая из функций ul) и иб?) имеет только одну особую точку находящуюся B точке 
пространства 1 или 2 и определенную любыми мультипольными истоками. По- 
данное в этой работе обобщение является билинейным соотношением между сле- 
дующими двумя системами коэффициентов: 


190 | А. Rubinowicz 
ОО ООО E E л Гм эээ = == 
1. коэффициентов, представляющих напряженности и фазы мультиполей 
обоих функций ul) и u в их точках истока 1 и 2, a также 
2. коэффициентов в разложении этих функий ul) и ul) на сферические и бес- 
селевые функции в точках пространства 1 и 2, где эти функции регулярны. 
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SOLUTION OF THE FLUCTUATION PROBLEM IN THE 
CASCADE THEORY BY MEANS OF THE G— EQUATIONS 
OF JÁNOSSY WITHOUT IONIZATION LOSS 


By Jan LoPuszaANsKI 
Boleslav Bierut University, Institute of Theoretical Physics, Wrocław 
(Received June 29, 1954) 


A new method is given for solving the G-equations of Jánossy. The plane of the 
independent variables (x, €) is divided into bands: 1/(п + 1< = < 1/n (n — natural 
number). The solution for е > 1 may be easily found. Taking this solution into account, 
it is possible to linearize and solve the G-equantions for в — 1/2. In each subsequent 
band (in the direction of increasing n) the equations can be linearized and solved in the 
same manner provided one knows the solution for the preceding band. Thus one obtains 
recurrent formulae together with explicit expressions for them for nucleon and electron- 
photon cascades. Graphical representations of the solution can be given. It is also shown 
that this method permits to obtain immediately the factorial moments of the distribution. 


1. Introduction 


Jánossy (1950a) has given a mathematical formalism describing the development 
of the nucleon resp. electron-photon cascades. In his paper the case of approximation 
A (without ionization loss) as well as the case of approximationt B (with ionization 
loss) is investigated. The fundamental equations of this formalism are called G-equa- 
tions. This formalism was developped in a series of papers (Jánossy and Messel 1950, 
1951, Jánossy 1950b, Messel 1951a, 1951b) which dealt almost entirely with the calcu- 
lation of the two first moments of the distribution. Little attention was paid to the 
distribution function itself. Later Messel together with Potts and Gardner solved the 
problem completely in a series of recent papers (Messel and Gardner 1951, 1952, 
Messel 1952, Messel and Potts 1952a, 1952b, 1952c, 1952d). The starting point were 
the linear equations of the so called general number distribution function and not 
the non-linear G-equations. Their method gives an indirect way of obtaining the 
solution of the G-equations. 

In the present paper we give a straightforward method yielding solutions of 
the G-equations for nucleon or electron photon cascades as well as the factorial moments 
of the distribution. We confine ourselves to approximation A. The case of approxi- 
mation B will be dealt with in a later paper. 
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It is very difficult to get numerical results from our solutions. We hope, however, 
that this solution will enable us to push forward the qualitative considerations of 
the distribution as well as to find asymptotical formulae for large depth of the absorber. 


2. The nucleon cascade 


We confine ourselves to the case of a nucleon cascade in homogeneous nuclear 
matter. The dynamical element in this cascade process is the cross-section for nucle- 
on-nucleon collision. We postulate with Janossy and Heitler (1949a, 1949b) that the 
collision cross-section is the same for p-p, p-n, and n-n collisions and that it is a homo- 
geneous function of the primary energy Еу and the secondary energies Е’ and E” of 
the colliding particles. We denote the cross section by 


w (Ey, Е’, E") dE’ dE" = о (=, =”) de’, de”, € = = DE = ЖЫ (1) 

It satisfies the relations 
DES Ж ap etes (2a) 
w (=, €”) = 0 for р 407 (2b) 


The total cross section 
171 


| f (е, e) de’ de" = а 
0 


0 


is a constant. The collisions are non-elastic. We do not distinguish between protons 
and neutrons. 
The fundamental equation: (G - equation) reads 


11 
8G 
5 (ina) абы = |е (5.x) (Eux) we, e) de ar, (3) 
moe 00 


where 
„Св, и, x) = Ф (e, n, x) и", 


Ф (в, п, х) denoting the probability to encounter n nucleons of the cascade with 
energy > gE, and an arbitrary number of nucleons with energy < sE, after the depth 
x, the cascade being due to one primary particle of energy E,. The boundary condi- 
tions are 


u for E 
1 fors > 1. (4) 


For & > 1 equations (3) and (4) will be of course satisfied by 


660) =| 


G (в, u, x) = 1. | (5) 
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If we know (5) we can proceed to solve (3) and (4) for e> 1/,. We introduce the aux- 
iliary function 


Г, (ғ, и, x) = 1— G(e, и, x) 
and put Г, inte (3); thus we find 


Een, o- ffarar ut ae n (= ) «n (=) 
- ffa de" ш en =") Г, (5) Г, | +! 
0 0 А 4 


We can write the last term on the right hand side schematically as follows 

Led POSEE | Чыл 711 

piso = Г [л Abts! 
The first three integrals vanish, as in the subintegral product at least one Г, vanishes. 
The fourth integral vanishes too on account of (2a). Equation (6) reduces thus to 


i 
Potato = fw an (5). (7) 
0 


1 . 
where w(e) = Ј de”! {w (e’, €”) + w (£", €)}. By means of (7) one may extend 
0 


the definition od Г» for all values of =. Equation (7) is similar to the equation for 
the average number of nucleons with energy > &E, (Messel, 1951a). To solve (7) 
we apply the well know method of Mellin transforms (abbreviated hereafter to 
M.T.): 


1 
ree Г, (eu, x) de = he (s, u, x) 
б 


We find 
: с +100 
Гг = S fe £i р (LL ш) е я (Tum pra). (8) 
271i 5 
where S(e, x) denotes the average number. We choose >> 0. То evaluate the integral 
we can make use of the numerical results obtained by Messel (19513). 
To solve our problem for e > !/, we introduce for this interval Г, = 1— С (eu, x) 
and apply the same ANA as in the case of the Г, function. We obtain thus 


the equation 


а Г, (s) J^ 5)» (2) de’ = — ІШЕ Дт ‚(= w (='=”) dede" 
(9) 
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This equation is similar to the equation for the second factorial moment (Messel, 
1951a). Generally, for e>1/n we get for Г, = 1—G an analogous equation to (9): 


n 


1 171 
950 апо fr (5) weyae = = ffr (S) (5) «e. eae ae. 
9x € E E 
4 2% (10) 


where Г, 1 is assumed to be а known function. By means of the М.Т. (Messel, 1951a) 
we get the solution of (10) for => Үп. 


С (в, и, x) = 1 —(1—u) S (e, x) | (11) 


x i" fai 
, "Hn , [4 95 € "o 
+ faz fan ff ae de” ш (é', є (- ea е к 
0 0 00 


This recurrent formula enables us to express G in terms of the first moment and 
the cross-section only. 

In a similar manner we can deal with the case of the nucleon cascade in a finite 
absorber. If we denote by w, (&1,... е) de,...de, the cross-section for nucleons coming 


out of the bombarded nucleus, than the equation for Г, reads 


1 


аго f[n(5)s«e-Ress 0 


where 


ый 1 1 
Rene) 3 | а-а fp PE (=) Wy (E45 :.., Ey) 
{2 j= 
0 0 


y 22 e 
E 1 1 
TOT 
a= lo | wder... de, 
RI PB 0 


The expressions for w, are given in the paper of Janossy and Messel (1951). Although 
equation (12) is linear and we can solve it successively, this task looks forbidding. 


and 


3. The electron-photon cascade in full screening approximation 
We denote by w® and w® the well known Heitler-Bethe expressions for brems- 
strahlung and pair production cross-section respectively. We have w® (г) = 0 for e> 1 
1 
and a, = | ш9 (e) de. 
б 


The fundamental equations read now 
: : 


9 


ea) (1—2) 6-12 09‏ 860+ رو 


€ 
0 


Fluctuation Problem in the Cascade Theory 195 


where 
GO (e, и, 1, из, x) = ХЕ У DE eae Т КО 
nmn 
and the superscript denotes for i = 1 that the primary was an electron, for 
i = 2 — a photon. 
The boundary conditions become 
и, Юг e< 1 


| we ЛА ds 


G9 (e, u, из, 0) = | 


Applying the same method as in section 2 with the auxiliary functions ГӘ, the equations 
(13) become for e > Ш 


1 
аге (2) , ; € € 
n pa E @) („/ , (C=) RIRES (1) 
ІЗ + TP (6) — fu (eae [r$ (5) + ro (12) 
0 


1 
= — f их (e’) de’ DEP H Lo , = à $ (i = 1,2). (15) 
0 


We solve (15) by means of the M.T., following Jánossy and Messel (1950) and, taking 


into account (14) find finally, for €> a 
n 


2 
GO (e, uy, uy 2) -1— У S99 (6а) (1-ш) 
k—1 


x 2 1 1 р 
de’ 95") € (am) | € 
у = = (m) | — — 
+ fa er all | 8n е die B ix Medi аз 


an i — uus :) ; (16) 


where 59” (e, x) denote the first moments. 


4. Discussion of the results 


Every term of our solutions (11) and (16) may be represented graphically in 
9S 
form of a „tree“. The „branches“ of these „trees“ correspond to (35) ог S, the 


knots“ correspond to w in the subintegral functions. Each ,,tree“ represents a certain 
cascade type. These graphs resemble the graphs of Feynman in the field theory. 
It will be noted that for = > п there can appear at most (n— 1) particles. Consequently 
for => п the terms represented by graphs with more than (n—1) end ,branches" 
are equal to zero. The „branch“ as e.g. 542 correspond to the following scheme: an 
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electron of energy >=, emits a photon of energy > € E, and retains an energy <ek,. 
The number of possible types of „trees“ will be estimated in a subsequent note. 

The solutions (11) and (16) enable us to find in a simple manner the factorial 
moments of the distribution. Indeed, (11) and (16) are represented in form of a single 
and double power series of (u—1) ог (ш —1) and (u,—1) respectively, and thus the 
coefficients of these series are the factorial moments. 

The solution of the G-equations of Jánossy obtained by the indirect method 
of Messel and Potts can be shown to be consistent with our solution. The 
corresponding calculations are rather tedious. 


The author wishes to express his sincere thanks to Professor J. Rzewuski 
and Professor R.S. Ingarden for valuable discussions and advices as well as to 
Komisja Popierania Twórczości Naukowej i Artystycznej for support of this work 
by a grand. 


КРАТКОЕ СОДЕРЖАНИЕ 


Лопушанский, Решение проблемы флуктуаций в теории ливней с помо- 
щью С — уравнений Яноши без ионизационной потери 
Предлагается новый метод решения G — уравнений Яноши. Площадь неза- 
висимых переменных (X, $) разделена на полосы: 1/(n+1)<¢=1/n (n — натураль- 
ное число). Решение при є->1 не представляет особых трудностей. Yura это реше- 
ние, можно произвести линеаризацию и решить С-уравнения при g-!/» Во BCA- 
кой последующей полосе (в направлении возрастающего п) уравнения могут 
быть линеаризированы и решены таким же образом, в предположении что мы 
знаем решение для предшествующей полосы. Таким образом мы получаем ре- 
куррентные формулы вместе с их явным видом в случае нуклонных и электронно- 
фотонных ливней. Могут быть поданы графические представления решения. 
Показано также, что этот метод дает возможность получить непосредственно 
факториальные моменты распределения. 
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À covariant configuration theory of bosons is developed on the example of photons. 
The first part of the paper contains a configurational representation of the state vector. 
For this purpose the mixed representation of the Günther type is used in which N photons 
are treated configurationally while the electrons and the remaining photons may be treated 
in any other way one likes. In the second part of this paper a configurational equation 
is derived for photons, corresponding to the Salpeter — Bethe equation. 


1. Introduction. In the theory of elementary particles there prevails lately a tendf 
ency to explain the properties of these particles by introducing new .degrees o- 
freedom; this changes radically the equations of non-coupled fields (Yukawa 1950, 
Pais 1953). The discovery of the possibility of a connection between isotopic spin 
in the Pais theory and the internal degrees of freedom in non-local theories (Rayski 
1954) enhances the attractiveness of the above trend of research. 

It seems, however, that too little attention has been paid to another possibility, 
namely that of considering the elementary particles as bound states of simple particles.*) 
This possibility does not exclude а priori the possibility explored by the prevailing 
tendency: both may supplement each other. 

It would appear that the exploration of the second possibility qe be based 
on a covariant configurational formulation of the equations of motion also for bosons. 
Such formulation is given in this paper onthe example of photons. All arguments 
contained in this paper are applicable, mutatis mutandis, to x mesons. 

The author has deliberately chosen photons for the object of his investigation, 
in order to stress his opinion that photons do not present more esential difficulties 
in the configurational treatment than mesons. Let us recall in order to avoid misun- 
derstandings that tlie position operator of a particle (e.g. of a photon) has ex definitione 
eigenstates given in the Schrödinger representation by ö-functions. These states are 
not identical with the localized states in the meaning of the paper of Newton and 
Wigner (1949), which do not exist for photons. 

| *) The author understands that this idea is now Günther's point of view. 


(191) 
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2. Configurational representation of the state vector. We denote by a the 
complete set of commuting observables describing the state of a free photon. Let a, 
be the annihilation operator of the photon in state |a>, N,— the number of photons 
in state |a», and ¥ [0] — the state vector in the interaction picture of the physical 
system consisting of the electromagnetic field and the electron field. By al, а?,... 
we shall denote labels a for arbitrarily numbered states |a>. For shortness we shall 
write №, = М, (r= 1, 2, кү 

When N photons are present in the field, then 


уюу = 


Let us denote by a,,..., ay the labels a of states la occupied by photons 1, ..., N. 
Of course, in the set of N labels: (a4, ..., ay) some labels may occur more than once 
(as photons are bosons). 

It is well known (Fock 1932) that in the representation defined by the base vectors 
IN, IN... >; where N ZO E2) the state vectors Бо) has the 
components (1) 


V (NEN) [Gea E qe c OS PAU. 


у 
In the representation given by the base vectors |в. save On мпегека оо 
(i= 1,2,..., №) and N >0, 1,2, ... (and where it is agreed moreover that for 
N — 0 we have the vacuum state of free photons: 00. ..—) the components of the 
state vector Y[o] have the form 


1 
P (a; ... ам) [a] = <a... ex | To] >= үлі «00...|a, ... ам Plol>, (2) 
The numerical coefficient on the right hand side of this formula is the consequence 
of the normalization condition: 


1= «WW. = У <|¥ (N,N,...) [6] > = 3 Y <1%(а...ак) (11. 


№, М... N а...ам 


The standard kets and bras appearing in the above sums pertain to the degrees 
of freedom of the physical system under consideration other than those of the electro- 
magnetic field, which means that they pertain to the degrees of freedom of theelectron 
field. They must appear in the formula because the components (1) and TO, are still 
abstract state vectors with regard to the electron field. 

Let us assume that the operators of the set a commute with the energy of the 
free photon. We may then describe the electromagnetic field in the interaction picture 
by means of the operator 


Pa (9 = У a Pra (х), (3) 


a 
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where P,a (x) = p (x) exp (—iE, t/h), Qux) хша Rand ur 2 oer es 
the vector index; E, is the energy of the photon in the state |a>. The functions 
Pua (х) describe the states |a> in the Schrödinger representation of the photon. 

We shall now represent the state vector ¥ [o] in the Schrödinger representation 
of the photons in a relativistic covariant way. We denote by x, ..., Xy the four-dimen- 
sional positions of the N photons placed on the space-like hypersurface о, and by 


My, ..., Ay the vector indices of these photons. Let H denote the Hamiltonian of the 


free photon, then |xu> = exp (iHt/h)|xu> are the eigenvectors of the position 
operator of a photon in the interaction picture. We have further 


ш> = У lao $, (9), 


a 


and hence 


bi Ay Un > = 3 |a ... ам > Фиш (дд)... Ф*имам (2n). 
an 
Consequently, in the representation given by the base |x, ш... Xylty>, where N = 
— 0, 1, 2, ... (for N — 0 we have the state of the vacuum of the free photons) the 
components of the state vector ¥ [o] take, thanks to (2) and (3), the following form 


Ч. им (а...) [o] = XM Ду, |W [c] > 
1 Е 
= ум «00 ... | gu, (а) ++ Pan %y) P |01 >. (4) 


When for example the set a of commutig observables comprises the momentum 
and the linear polarization of the photon—then, using the normalization cube of 
volume V, we have in formula (3) 


l а 2 2 
Pua (X) = VE е, ee, wheree, = 1,  k,-—0- (5) 
V 
Instead of describing the photon field in the interaction picture by using the 
operator Ф, (x) given by (3) it will be more convenient to use the operator A, (x) 
(electromagnetic 4-potential). Let us note that we have in virtue of (5) 


A, (x) = Af (x) + Az (9, 
where ; 


А : 
АТ (x) = Y EJ Ga Qua (X) = + A, (x). (6) 


Here „+“ is proper for u = 1, 2, 3, and „—“ for u = 4 (we use a notation in which 
A, (х) = ig (x) ). Thanks to (3), (5) and (6), we have for any operator © independent 


eorex: 


free = fararo- i) aa, (т) 
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because —(9/9x,) exp (kx) = = k exp (ik,x,) (in the notation used x, = 1% = ict, 
k = ikọ, where ky = |k| = À). 

Using (4) and (7) and a flat hypersurface о, we may write for any two arbitrary 
states Y [о] and V? [o] 


Гар- Фа, EG SD RE Te] (8) 
, 2\N 9 9 
= Mene ... XN)[0] (5 9 E BER FD ea … ÆN)[0] 


(x; = (xu) (4 = 1, 2, ..., N)), where we have introduced the notation 


1 
Fa, NA О == > «чу | А} (x) 27 A} (xn) P [o] 


Үл! 
= 16 2200.21: At). A) (>, (9) 


VN! 
and the points x, ..., xy are lying on the space-like hypersurface c. The symbol :(...): 
denotes here the ordered product of operators, in which all the creation operators 
are placed to the left of the annihilation operators (Wick 1950). 
Formula (8) may be covariantly transcribed as follows 


* 2 x 
n ... 930,4 PC ... хм№)[0] fa dy, ... dy FO an -.Xw)[o] (10) 


where d,— — i (2/0x, — /ax,)/2 (cf. Rzewuski 1951) and (d30,) = (dx, dx, аж, 
аха dx, dx,, dx, dx, ахз, —idx,dx,dx;). 
Formula (10) gives the rules of normalization of function (9): 
А "M 
= уз da, de, E, uy ... Xx)[o] iss) dy, ... do Fa a + eno] >. 
We observe that the configuration function (9) is not the amplitude of proba- 
bility density for findig photons at the point жі, ..., xy on о (is not a wave function), 


in which it differs from (4). However, on the basis of (10), it may serve to calculate 
the scalar products of states, as follows 


<-> TED 
: EIN 
гэ? fan. ... 30, „<F won XN) а 8 d, . MC D А -- бл...хм)[о]:>. 


Ав in the function given by (9) there appear the operators А (x), this function is 


A rwn 
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QI U ا ا و‎ in 


more convenient than the wave function (4) or the configuration function introduced 


for photons by Landau and Peierls (1930). 


Functions (4) and (9) are obviously still abstract state vectors with respect to the 
electron field. The standard kets and bras appearing on the right hand sides of the 
normíalization condition and formula (11) stand in relation to this field. 


In the following considerations we shall use a different representation, namely 
the mixed representation in which a certain fixed number N of photons will be 
described configurationally, while the remaining photons will be described by means 
of occupation numbers (cf. Günther 1954). The mixed representation is defined when 
N is fixed by the base |N; N; ... x, дү... хушу. In this representation the com- 
ponents of the state vector ¥ [о] have the form 


n a (V1: № os Ху... xy)[0] = < ММ, ..., жил + хуи | Y[o] > (12) 


— «NN, ... | eu (n) <. Фик(им) [6] > = «NN, ... | Puuna --- XN)[o] >, 


m ! 
where the following notation has been introduced (not to be mistaken for the function 


(4) )): 4 
Уы -- хм) Lo] = ум gu) --- Puyn) P [o]. (13) 


The points x, ..., xy lie on the space-like hypersurface 

In view of the completeness of the system IN iN, ... > vector (13) is the configura- 
tional wave function of the N photons which were picked out by us. On the other 
hand, inthe interaction picture the same vector plays the role of the abstract state vector 
with regard to the electron field as well as with regard to the „non-configurational“ 
photons. 

Instead of the configurational wave function (13) we shall use now the following 
configuration function (not to be mistaken for the function (9)!) 


E a aad bai. = АХ бы)... А, Cn) F [o]. (14) 


This function satisfies the following normalization condition 


l= <YY> 


2\N 
= fa 0,,.. dg ON < Sr ... Xy) [0] (2) ага, Fun (а... XN)[0]2-- 
0 . 
The rule for the calculation of scalar products for any states with N „configurational“ 
photons is given by the formula 


LPO YO > 


! a z 
= J ds oni dy Ory FAY, (а... x) Lo] 2 d, dy, РОЧ MES wA 
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Here the standard ket and standard bra on the right hand side pertain equally to 
the electron field and to the »non-configurational« photons. 

In passing to the Heisenberg picture, we drop the functional dependence on o 
from the configuration function (14): 


Fa..un (x, ‚з х) = [c, = oo] Fun (хі 22; ху) [c] 


1 
= ум An, (ол)... Aux (xy) P [— ое], (16) 
where 
О (x) = 0-5, - оо] О(х) U (б, — оо], 
«e MA qu 
Tm U^! [в, — œ] V [o] = У [— оо]. 
Here 
ihe s = H(x) U[o, — оо], U[— оо, —00] = 1 (18) 
and 
H (3) = р) 4.6), ju) = iec: yy la): (19) 


We use the antysymmetrized current and the notation (у,) = (— ifa, В), w(x) 


= pe. 
Observe that (Schwinger 1948) 


ead rs 9 
Әх, A, (x) = Оо, — оо] Dx, Аи (x) Шын (6, — оо], 


from which it follows that in formula (15) we can instead of function (14) use 
function (16). 
Differentiating functionally the relation 


Fa. ux waning) [o] E Ufo, = со) Frun (26 к) 


Ы 


with regard to the hypersurface c at points x different from x, (i = 1,2, ..., № 
we get in virtue of (18), | 


ihc ÖFu...un (хі Ge XN) [c] 

до(х) 
From this we can see that in the case of the physical system containing permanently 
exactly N photons, the configuration function (14) is a »pure« function of the points 
х1, +++) Xy (is not a functional of о), i.e. does not change when с is deformed at points 
different from the 4-positions x, ..., xy of the N photons. Therefore in such cases 


= A(x) Раим (ae xn)le]- 


In (о, — co] Faun (%4-+-%y) = Fi un&---&N) 


because U[— oo, — oo] — 1. The state of the system of N photons characterized 
by the above condition will be called a »stable« state. 
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Generally, however, function (12) is a functional of с because of the possibility 
of creation or annihilation of electrons and „non-configurational“ photons at any point of 
the surface o. Function (16) is, on the contrary, always a »pure« function of the points 
Аа 

In formula (16) the points x, ..., xy lie on the space-like hypersurface с. We 
can therefore transcribe this formula as follows 


1 
NE, PRG EE ось eo == VNI TA (He Aun (х) ¥ [— со], (20) 


where the symbol T(...) is a chronological product as defined by Wick (1950). We 
adopt formula (20) as the definition in the Heisenberg picture of the configuration 
function of N photons for any relative position of the points 34, ..., x (cf. Günther 
1954). 

The component of the vector (20) corresponding to the true vacuu of electrons 
and photons (i.e. the vacuum of the interacting particles) describes a system comprising 


N photons. It has the form (cf. Gell-Mann and Low 1951) 


Tra eR Van y | Dan @ AN) > 


1 
^ VN! 


because the relation between the true vacuum state | V> and the state of the vacuum 
of the free particles lo is 


|V > = Ищю, =] | 0>, 


vacuum of E vacuum of » 

free photons// | free electrons // ' 
In particular, we have for a stable state 

«V | о, — оо] Fu..un (...®)> = «V | Е.м (Ху... хм) > = Faun (Vo. x). 


Formula (21) can be written, due to (17) and (18), as follows (cf. Gell-Mann and 
Low 1951, Dyson 1949): 


о | Доо, —oo] T(A;; (x... 4, (xn) ?I-»]>, (21) 
H H 


vaccuum of N 
free photons/ 


ОО 


Ри... “М (Уху. . ЖАЗ 


= <, oo] U" el (аа) U[os ,09].--A gn (хм) Оол оо] [оо] > 


+00 
E P | T = i) | dxy...dxyT(H (x1)... H (en) Ap, (ху)... Agen) оо] >. 
M n=0 кы (22) 
In the above formula the space-like hypersurfaces 9, ..., 0 pass without intersecting 
through the points x,, ..., xy respectively. 
Formula (22) Жа 1 the construction of the configuration function of N photons 
made of known matrix elements corresponding to the irreducible or reducible Feynman 
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diagrams, ordered according to powers of the coupling constant. We may restate 
the construction, given explicitly by formula (22), by an integral equation with a kernel 
containing only terms corresponding to irreducible Feynman diagrams. In this way 
we get an equation for photons corresponding to the Salpeter — Bethe equation. 
We will find this equation in the case of the system comprising permanently two 
photons, i.e. for the stable states of two photons, restricting ourselves to the equation 
with a kernel consisting only of terms corresponding to the lowest irreducible diagrams 
of interaction of two photons. Other terms can be easily taken into consideration. 

3. Configurational equation of two photons. Calculating the right side of 
formula (22) by Wick’s method (1950) in the case N = 2, we find, thanks to (9), 
that in the expansion according to powers of the coupling constant the terms giving 
the interaction between the two photons appear for the first time in the fourth order. 
After calculating the transition amplitude we see that these terms correspond to the 
diagrams shown in Fig. 1, and each of the terms appears 12 times. 


N NUS. ZEN VEN SS ZEN ^ 
/ Ss 2%” NE 7 N #2 S 22 55 
Fig. 1 


Each of the first four diagrams, and each of the last two, bear equally on the 
transition amplitude in view of the symmetry of the configuration function of two 
photons. 

If we consider only the above terms, formula (22) takes the following form 

F(a) = Ее, (Рза) + | dr du, d, dz, D" (x — ху) D't! (из — ху) 

X Си, (12320) Fou (VX3%4), | (23) 
where 


Сима, (%1 Xa Xg %4) 


e? i 
E i 2 tr [2y,,S (а — x9)yu,S (а — жз)уи S 4 (X3 — ма). S + (4 — ху) 


+ Yu (п — хз) کے‎ (з — №) VpS (хә — 24) у, S — ху)]. (24) 
The function 


1 
FY (Иж) = < 00 | А (x) A À (x) 7 [—co] > (25) 


ү? 
is the configuration function of two free photons. In formula (23) the retarded function 


D,(x—x) t»t 


D(x — х) = 15 — (26) 
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appears as the result of the relation 
A (x) A (x) = <0 | T(A,(x) Aj (x')) | 0 = hcó, Па — x). 
The recurrent dot over the operators denotes their contraction (Wick 1950). The 
S,(x—x) and D (x — x’) are the Feynman functions given by the relations 
Papa) = «0 | Тір (ура) |O> = Sax — x), 
A (а) АД) = «0 | TCA, Q) A (^) |0> = вед „О. (x — =). 


They satisfy the equations 
9 
(» Оа + +) 9+(х == x’) = — ió(x == s L1 D „(х == х) => ió(x — ж). 
u 


Let us note that the function D''(x — x’) satisfies the equation 
O Die — x) = ió(x — x’) (27) 
and that it has the following property 
At*(x) x sooner than a 


2 fas Aj *(x')d’, D' (x — x) = | 0 x later than o сз) 


(where d; = Ц0/Әх, — 9/9x,)]2. 
From formula (23) we obtain the integral equation by substituting the function 
of the free photons in the integral by the unknown function: 


+00 
Fun (7) = Е ә а(х) + il dx, dx dx, dx, DI Ge х) DEC —52*) (29) 


X os аза (21232314) Е, Lalla (Ухх). 


This is the configurational integral equation describing in the first approximation 
the stable states of two interacting photons. It restates the construction of the 
configuration function of two photons made of terms ыер; to the reducible 
diagrams consisting of the diagrams shown in fig. 1. 

We find the probability amplitude for the system of two photons by means 
of formula (15) (where N — 2). Let us calculate, for example, the probability amplitude 
for finding — in the stable two-photon state У [c] on the hypersurface с — of two 
free photons characterized by ФФ. By making use of (15), (29) and (28), we get 


«vo Pol = = fo, 430», < TA (x1%2) (F .) dd, Fan (xg) [o] > 


МЕ D) (Vars) (2 | d, dy Py, (Ух, ху) 
L. 
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ӘЛЕ 
= | do, do, FO) (Ух, хз) (2) d,d, F®, ( Vx, хә) 


a 


a +00 
+ fa dx, [ FO (Рах) Gum (халал) Fou, (Vosa). (80) 


For the scattering processes of photons on photons, we move с to + оо. In considering 
the case of single scattering, we substitute for F, „, (Ухх) the initial function of 
two free photons: Е, (Их хз). 

From (29) we derive the integral equation for the investigation of bound states 
of two bosons by omitting the free term ЕО, (Vx1x2) (cf. Salpeter and Bethe 1951). 

The integral equation (29) can be reduced to two integro-differential equations 
of the second order, or to one integro-differential equation of the fourth order. Thus 
in virtue of (27), we have 

+00 
AF pu (Ра) — i Í da, das dx, (х. — ж) Gr Cara ES) FT Vx. x,), 
55) 


+00 
ОР, ЕЛЕР) = 1 f ds dx dx DI (xy a xj) Green) FRERE) 
or 
+оо 
0, OF, „(Их x9) = — [ аа, СЛО) ЕЗ (К): 
—00 


The above integro-differential equations are obviously richer in solutions than 
the integral equation (29). 

The author is indebted to Dr M. Günther and Prof. J. Rayski for valuable 
suggestions and discussions. 

Note added in proof. Recently a paper has been published by K. Hain (Z. Natur- 
forrsch., 9a, 495 (1954)) concerning bound states of x mesons. In that paper mesons 
are also treated in the configurational way, but the question of normalization of 
the configuration function is not investigated. 


КРАТКОЕ СОДЕРЖАНИЕ 


Круликовский. Конфиурационное представление фотонов и уравнение 
их связанных состояний 

В этой работе автор формулирует на примере фотонов ковариантную кон- 
фигурационную теорию бозонов. Первая часть работы содержит конфигурацион- 
ное предстваление вектора состояния. При этом применяется смешанное пред- 
ставление, в котором М фотонов рассматривается конфигурационно, а электроны 
и остальные фотоны любым другим образом. Во второй части автор выводит 


-конфигураци i 
фигурационное уравнение для фотонов являющееся аналогом уравнения Саль- 
петера и Бете. 
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FORTPFLANZUNG VON SPRÜNGEN ELEKTROMAGNETISCHER 
FELDSTARKEN UND EINDEUTIGKEITSBEWEIS FÜR DAS 
ANFANGSWERTPROBLEM DER MAXWELLSCHEN 
GLEICHUNGEN 


Von A. RuBiNOWICZ 
Physikalisches Institut der Polnischen Akademie der Wissenschaften, Warschau 


(Eingegangen am 29 September 1954) 


Sprünge von elektromagnetischen Feldstürken, die auf mit beliebigen Geschwindigkeiten in der 
Richtung ihrer Normalen sich bewegenden Flüchen auftreten, nehmen an der Energiebilanz eines vier- 
dimensionalen Raum-Zeit-Gebietes im allgemeinen als Energie-Quellen oder Senken teil. Vom physi- 
kalischen Standpunkt kann man diese Energieänderungen durch elektrische und „magnetische“ Ströme 
erklären, die in den Sprungflächen fließen. Die Maxwellschen Gleichungen zeigen, daß man sich diese 
Flüchenstróme entstanden denken kann durch Zusammendrängen von Leitungsstrómen in die Sprung- 
flächen. Diese Ströme müssen verschwinden, damit die in Rede stehenden Sprünge keinen Beitrag zur 
Energiebilanz des elektromagnetischen Feldes liefern. Auf diese Weise ergeben sich Bedingungen für 
die physikalisch zulässigen Sprünge der elektromagnetischen Feldstürken, die u.a. besagen, daß solche 
Sprünge sich nur längs charakteristischer Flüchen fortpflanzen kónnen. 

Auch bei Vorhandensein von Sprüngen elektromagnetischer Feldstürken, die die angegebenen 
Bedingungen erfüllen, ist der Eindeutigkeitsbeweis für das Anfangswertproblem der Maxwellschen 
Gleichungen durchführbar und zwar in einer solchen Weise, daß er der endlichen Ausbreitungsge- 
schwindigkeit der elektromagnetischen Wellen Rechnung trägt. 


$1. Problemstellung 


Liegt eine Lósung eines bei einer physikalischen Fragestellung auftretenden 
mathematischen Problems vor, so ist man noch keineswegs sicher, daß sie vom 
physikalischen Standpunkt annehmbar ist. Die Funktion f(x — vt) stellt z.B. eine 
Lósung der Gleichung der schwingenden Saite 


lO 9% 


28 908 дж 
даг, auch wenn die Funktion f(z) an irgend einer Stelle einen Sprung aufweist. Eine 
solche unstetige Lósung ist aber physikalisch nicht realisierbar, da sie ja eine mit 
der Geschwindigkeit v erfolgende Fortpflanzug einer sich momentan 6ffnenden und 
wieder schlieBenden Rifstelle darstellt. 
(209) 
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Zur Entscheidung der Frage ob eine vorliegende Lósung eines in der Physik 
auftretenden mathematischen Problems physikalisch zu verwirklichen ist oder nicht, 
muß man für das betreffende Problem einen Eindeutigkeitsbeweis erbringen. Ein 
soleher Beweis zwingt uns ja, eine Gesamtheit von mathematischen Bedingungen 
und somit auch von ihnen entsprechenden physikalischen Bestimmungsstücken 
anzugeben, die zur eindeutigen Festlegung des mathematischen Problems sicherlich 
hinreichend, wenn auch vielleicht nicht alle notwendig sind. Ist nun eine Lósung 
gegeben, die alle diese Bedingungen erfüllt und fordert keine von ihnen physikalisch 
nicht realisierbare Eigenschaften, so kann man annehmen, daf) die gegebene Lósung 
physikalisch nicht widersinnig ist. Allerdings liefert der Eindeutigkeitsbeweis selbst 
keine Kriterien zur Feststellung, ob eine im Eindeutigkeitsbeweis auftretende mathe- 
matische Bedingung physikalisch zu verwirklichen ist. 

Wir wollen uns hier mit den Maxwellschen Gleichungen beschäftigen. Der 
Eindeutigkeitsbeweis für das Anfangswert- und auch das gemischte Problem der 
Maxwellschen Gleichungen (Cohn 1900, Riemann-Weber 1900, Rubinowicz 1926, 
Plebanski 1953) wurde bisher nur unter der Voraussetzung geführt, daß ihre Lösungen 
nur in ruhenden Flächen Sprünge aufweisen. Vom praktisch physikalischen Interesse 
sind aber auch Lösungen der Maxwellschen Gleichungen, die im Raum-Zeit-Gebiet 
Sprünge in dreidimensionalen, charakteristischen Flächen besitzen. Im dreidimensio- 
nalen Raume sind dies Sprünge, die in Flächen auftreten, welche sich mit Licht- 
geschwindigkeit in der Richtung ihrer Normalen bewegen. Solchen Sprüngen elektro- 
magnetischer Feldstärken begegnen wir z.B. in dem vom Verfasser (1954) behandelten 
Problem der Fortpflanzung unstetiger elektromagnetischer Signale in Wellenleitern. 
Wir werden daher im folgenden den Eindeutigkeitsbeweis auch für den Fall erbringen 
müssen, wo solche Sprünge vorhanden sind. 

Um die notwendige Allgemeinheit der durchzuführenden Überlegungen. zu 
wahren, werden wir jedoch im allgemeinen annehmen, daß die Unstetigkeitsflächen 
sich in der Richtung ihrer Normalen mit beliebigen Geschwindigkeiten bewegen, 
die speziell auch gleich der Lichtgeschwindigkeit sein können. 

Die Bedingungen, denen die Sprünge elektromagnetischer Feldstärken genügen 
müssen, werden wir aus dem Erhaltungssatz für die Energie ableiten. 

Da es uns in der vorliegenden Arbeit darauf ankommt physikalisch realisierbare 
Eigenschaften elektromagnetischer Felder von physikalisch nicht zu verwirklichenden 
zu unterscheiden, müssen wir im folgenden auch von magnetischen Ladungen und 
Strömen sprechen, eben um ihr Vorhandensein ausschließen zu können. 


82. Der Energiesatz für das elektromagnetische Feld 


Zur ТИЕР der Bedingungen für die Sprünge der elektromagnetischen Feld- 
_ Stárken und als Vorbereitung zum Beweis des Eindeutigkeitssatzes beschäftigen wir 
uns mit dem Energiesatz. Die Maxwellschen Gleichungen setzen wir zunächst in der 
 homogenen Gestalt (mit verschwindenden Ladungen und Strómen) 
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OE 4 = > < 
55 -+ М ey LL div e E = 0, 
T с 2.1) 
oH а M S 
ис; = — rot Е, div uH = 0 


voraus. Hier bedeutet т = ct. Zur Vereinfachung der Ausdrucksweise, soll jedoch т 
als die Zeit bezeichnet werden. Die Materialkonstanten e, u und ø sollen dabei stetige 
Funktionen der Ortskoordinaten sein, so daß wir es hier mit einem isotropen aber 
inhomogenen Halbleiter zu tun haben. 

Um das Raum-Zeit-Gebiet zu definieren, in dem wir die Lósung der Maxwellschen 
Gleichungen betrachten, bezeichnen wir mit 


f(x, У, 2» T) =0 (2.2) 


eine Fläche, die im dreidimensionalen Raume sich mit der Geschwindigkeit 


D. 


с 
= (2.3) 
in der. Richtung ihrer Normalen bewegt!. Sie ist eine Lösung der Differentialgleichung 
des Huygensschen Enveloppenprinzips 


lee ne, (2.4) 


wo m bei den zunächst folgenden Betrachtungen irgend eine im allgemeinen stetige 
Funktion der Raumkoordinaten bezeichnet, die eventuell auch von der Zeit abhängen 
kann. Die Größe m kann im Einklang mit (2.3) als ein Brechungsexponent gedeutet 
werden. 
Im Falle 
m? — eu (2.5) 


stellt die Gleichung (2.4) die Differentialgleichung der charakteristischen Flächen 
der Maxwellschen Gleichungen (2.1) dar. Bei den folgenden Überlegungen wird jedoch 
die Giltigkeit der Beziehung (2.5) zunächst noch nicht vorausgesetzt. 


Es sei 
fax, ¥, 2) = 0 Ä ` (26) 


eine geschlossene Flüche, die zur Zeit v, ein dreidimensionales Raumgebiet R, be- 
grenzt. Es gibt dann zwei Lósungen (2.2) der Differentialgleichung (2.4), die zur 
Zeit vy mit der Fläche (2.6) zusammenfallen. Die eine von ihnen umschliefit ein immer 
größer werdendes dreidimensionales Raumgebiet, die andere ein zunächst kleiner 


1 Wollte man den folgenden Betrachtungen Flüchen (2.2) zugrunde legen, die sich nicht in der 
Richtung ihrer Normalen sondern in beliebigen Richtungen bewegen, so würde dies die folgenden 
Betrachtungen unnótig komplizieren. Von praktischer Bedeutung sind ja doch nur Sprünge elektro- 
magnetischer Feldstärken in, charakteristischen Flächen und diese re ja zur Klasse der von uns 
betrachteten Flächen, 
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werdendes. Diese letztere geht aber schließlich sozusagen durch ein Umstülpen in 
eine Fliche über, die ebenfalls ins Unendliche abwandert. 

Ferner sei R} das dreidimensionale Raumgebiet, das zur Zeit т, von der sich 
zusammenziehenden Fläche (2.2) begrenzt wird und das durch eine stetige Umformung 
aus dem Raumgebiet R, hervorgeht. Dabei wird ausdrücklich vorausgesetzt, daß т; 
einen Zeitmoment vor dem völligen Umstülpen der sich zusammenziehenden Fläche 
bezeichnet. 

Schließlich sei G ein vierdimensionales, kegelstumpfartiges Raum-Zeit-Gebiet, 
das von den dreidimensionalen, in den Ebenen t = 7, und t = 7, liegenden Raum- 
gebieten Ку und К, sowie von dem Teil F der sich zusammenziehenden Fläche (2.2) 
begrentzt wird, der zwischen den Ebenen т = Tọ und т = 7, gelegen ist. | 

Um die Energiebilanz im Raum-Zeit-Gebiet С zu erhalten, integrieren wir über 
dieses Gebiet den differentiellen Energiesatz für die Maxwellschen Gleichungen (2.1): 


t v IW 
div Site ra = — Q (2.7) 


(= Ex w= eF tum, о-в). 
Die Anwendung des vierdimensionalen Gausßchen Satzes liefert dann die Relation 
(S, cos (N,x) + 5, cos (N,y) + 5, cos (N,2) + cW cos (N,t)) dF 
= — (оф. ; (2.8) 
б 


em 


Hier bedeutet dF ein dreidimensionales Flächenelement der Bn s 
von G und 


dv = dxdydzdt (2.9) 


ein Raumelement des vierdimensionalen Raum-Zeit-Gebietes x, y, z, v. N ist die 
äußere Normale an die Begrenzungsflichen von C. 

Um die Beziehung (2.8) zu deuten, müssen die hier auftretenden Flüchen- 
integrale berechnet werden. Wir bezeichnen sie alle insgesamt mit J und es seien 
11, und Ig die Anteile, die die Begrenzungsflüchen von С nämlich Ку, К, und F zul 
beisteuern. Es ist demnach 


| Ihn SED | (2.10) 
Für die drei hier rechts stehenden Integrale erhalten wir die nachstehenden Aus- 
drücke: 
1) In der Ebene Ry d. h. T = ty ist dF = dx dy dz und cos (N, x) = 
= cos (Му) = cos (N,z) = 0, cos (N,t) = —1. Daher wird 


h= — [е |, dedyde | (2.11) 
Ro { 
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2) In der Ebene R; d.h. t = т, wird analog, da hier cos (N, т) = 1 und alle 
übrigen Richtungskosinus verschwinden, 


Г = 
УКИ 


т, ах dy dz. (2.12) 


3) In der durch (2.2) dargestellten Fläche F ist es notwendig die Verhältnisse 
ein wenig näher zu betrachten. Die Richtungskosinus an diese Fläche werden mit 


Rücksicht auf die Differentialgleichung des Huygensschen Enveloppenprinzips 
(2.4) durch 


cos (N, x) — 


OF FB ARR ПЕТ Л ДИН 
КИ НУ тү uL. i DA = 
CEFR +SEE FSA OF mA 
gegeben. Die ersten drei dieser Ausdrücke lassen sich durch Richtungskosinus im 


x,y,z-Raum ausdrücken. Die Fläche (2.2) wird ja durch die dreidimensionale Ebene 
T = const in einer zweidimensionalen Fläche, die wir mit Ф(т) bezeichnen wollen, 


(2.13) 
cos: (№ т) = 


geschnitten. Sie wird durch (2.2) gegeben, falls wir hier 1— const setzen. Die Richtungs- 
kosinus der Normalen п an die Fläche Ф(т) sind mit Rücksicht auf (2.4) durch 


f fx 
i а. 2.14 
AFT mf ae. 
bestimmt. Es wurde dabei in (2.14) das negative Vorzeichen verwendet, damit die 
Normale n die Richtung angibt, in der sich die Fläche Ф(т) fortbewegt. N stellt ja 
die äußere Normale der Fläche F dar, so daß ihre Projektion auf den x,y,z-Raum 


die zu n entgegengesetzte Richtung aufweist. Mit Rücksicht auf (2.13) erhalten wir 
daher 


cos (n, x) — 


cos (Л, х) = —m cos (n,x) cos (№, т), ... . (2.15) 


Um noch das Flächenelement dF der Flüche (2.2) anzugeben bezeichnen wir . 
im dreidimensionalen x, y, z-Raum mit do das Flächenelement der zweidimensionalen | 


Fläche Ф(т). Es wird dann 
| do dt 


“sn (МЭ. 


dF 


Das dreidimensionale Flächenelement dodt kann nämlich als eine Projektion des 
dreidimensionalen Flüchenelementes dF angesehen werden, wobei der Winkel zwischen 
den beiden in Rede stehenden Flüchenelementen 7/2 — X (N,v) beträgt. Mit 
Rücksicht auf (2.13) wird daher 


: m pde ducc | | (2.16) 
А m cos (N, т) 
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Somit ergibt sich wegen (2.3) und (2.15) 
Ip = | dr f (—S, + vW)do. (2.17) 


To P(r) 


Die Beziehung (2.8) erhält somit mit Rücksicht auf (2.9), (2.10), (2.11), (2.12) 
und (2.17) die Gestalt 


| ~ 
- [v dx dy dz |w dx dy dz + = dr |5, + um do 


Ro Rı To Ф(т) 


-- BL ахау аат = 0. (2.18) 
G 


Diese Relation stellt die Energiebilanz für das vierdimensionale Raum-Zeit- 
Gebiet G dar?, da sie in der folgenden Weise interpretiert werden kann: 

Die zur Zeit Tọ im Raume Rọ enthaltene Energie wird verwendet um die nach- 
stehenden Energiebetráge zu liefern: 

a) die Energie, die zur Zeit , in dem dreidimensionalen Raume R, enthalten ist, 

b) die Energie, die durch die sich zusammenziehende Begrenzungsflüche Ф(т) 
in dem Zeitintervall von v, bis т, in der zur Bewegungsrichtung n von Ф(т) ent- 
gegengesetzten Richtung — n infolge der Energiestrómung (Integrand —5,) und 
infolge der Bewegung von Ф(т) (Integrand vW) durch diese Fläche hindurchgeht, 

c) die Energie, die zur Erzeugung der Jouleschen Wärme in dem Raum-Zeit- 
Gebiet G benótigt wird. 

Bei dem im 84 folgenden Beweis des Eindeutigkeitssatzes wird für uns der 
Fall von besonderer Wichtigkeit sein, wo alle in der Energiebilanz (2.18) auftretenden 
Integrale nicht negative Energiebetrüge darstellen. Aus der Tatsache, daü W>0 
und О > 0 ist, folgt daß die beiden ersten Integrale und das letzte Integral in (2.18) 
nicht negativ sein kónnen. Um uns Klarheit über das Vorzeichen des vorletzten 
Integrals zu verschaffen, gehen wir von der Ungleichung 


g; We E — ul)? >0 
aus, aus der mit Rücksicht auf (2.3) die Ungleichung 


2 MS Sen 
à; (ЕЁ t BH) > тс = EH 


folgt. Unter der Voraussetzung, daß 


— Vin >m ` (2.19) 


? Das Auftreten des Faktors 1/c in (2.18) erklärt sich daraus, daß wir hier nach v = ct und nicht 
nach ¢ integrieren. 
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ist, ergibt sich daraus nach (2.7) 

vW > | Sa 
Die Ungleichung (2.19) stellt daher eine hinreichende Bedingung dafür dar, daß 
auch das dritte Integral in der Energiebilanz (2.18) nicht negativ ist. 

Damit also alle in der Energiebilanz (2.18) auftretenden Integrale nicht negativ 
sind, muß somit mit Rücksicht auf (2.19), die durch (2.3) definierte Geschwindigkeit v 
nicht kleiner sein als 

с 


Cu. (2.20) 


Dies ist die Geschwindigkeit mit der sich elektromagnetische Wellen in einem nicht- 
leitenden Medium mit den Materialkonstanten £ und u ausbreiten. Diese Geschwindig- 


keit vz, der nach (2.19) т = ew entspricht, soll іп den folgenden Überlegungen 
kurz als die Lichtgeschwindigkeit bezeichnet werden. Einer mit dieser Geschwindig- 
keit im dreidimensionalen Raum in der Richtung ihrer Normalen sich bewegenden 
Fläche Ф(т) entspricht nach (2.5) im Raum-Zeit-Gebiet x,y,z, t eine charakteristische 
Fläche F der Maxwellschen Gleichungen (2.1). 

Das Raum-Zeit-Gebiet G, das von einer charakteristischen Flüche F mit begrenzt 
wird, soll in den folgenden Überlegungen mit C, bezeichnet werden. Es enthält in 
seinem Innern alle Raum-Zeit-Gebiete G, die durch eine mit einer gréferen Geschwin- 
digkeit v als die Lichtgeschwindigkeit v, erfolgende Zusammenziehung der Be- 
grenzungsfläche des dreidimensionalen Raumgebietes Ry entstehen. Diese Geschwin- 
digkeiten entsprechen solchen Brechungsexponenten m, die der Ungleichung (2.19) 
genügen. 

Die Energebilanz (2.18) haben wir unter der Voraussetzung abgeleitet, daß die 
elektrischen Materialkonstanten ви und o in dem ganzen Raumgebiet К, stetig 
sind. Man kann sich aber leicht überzeugen, daß sie auch richtig bleibt, falls diese 
Konstanten in gegebenen ruhenden Flächen Sprünge aufweisen. Jeder solchen Fläche 
entspricht nümlich im vierdimensionalen Raum-Zeit-Gebiet eine dreidimensionale 
Zylinderfliche, deren Erzeugende parallel der Zeitachse sind. Bei der Anwendung 
von (2.8) ist das hier auftretende Flächenintegral über die beiden Seiten dieser Zylinder- 
fläche zu erstrecken. Da jedoch die Normale N an diese Zylinderflächen senkrecht 
zur Zeitachse steht und daher cos (N,r) — 0 ist, so wird sie mit der Normalen 
an die Unstetigkeitsfliche im dreidimensionalen Raum identisch, so daf) in den 
in Rede stehenden Flächenintegralen über die Zylinderfläche nur die Normalkompo- 
nenten des Poyntingschen Vektors auftreten. Die Grenzbedingungen für das elektro- 
magnetische Feld sichern jedoch die Stetigkeit dieser Normalkomponenten, so daß 
die Summe der beiden Flächenintegrale über jede Zylinderfläche verschwindet. 

Auch das Vorhandensein von vollkommenen Leitern beeinträchtigt nicht die 
Gültigkeit der Energiebilanz (2.18), da die Normalkomponenten des Poyntingschen 
Vektors an den Oberflächen dieser Körper verschwinden und diese daher keine 
Beiträge zur Energiebilanz (2.18) liefern. 


§ 3. Der Energiesatz und die Fortpflanzung von Sprüngen 
des elektromagnetischen Feldes 


Die Richtigkeit der Energiebilanz (2.18) wird im allgemeinen in Frage gestellt, 
falls das elektromagnetische Feld im Raum-Zeit-Gebiet G irgend welche Spriinge 
an einer bewegten Fläche Ф*(т), d. В. an einer Fläche (2.2) aufweist, die sich mit einer 
Geschwindigkeit v* = c/m* in der Richtung ihrer äußeren oder inneren Normalen 
bewegt. Dabei nehmen wir an, daß m* und daher v* von m bzw. v verschieden sein 
kónnen. 

Diejenige Seite der Ф*(т) im Raum-Zeit-Gebiet G entsprechenden Fläche F* 
von der die Normale N ausgeht, móge mit dem Index 1, die Seite aber nach der 
hin sie gerichtet ist mit dem Index 2 bezeichnet werden. Es seien ferner 


E, = E, Н, = Н und E, = Е + AE, Н = Н+ AH 
die Grenzwerte der elektromagnetischen Feldstärken, wenn wir uns der Flüche F* 
von der Seite 1 bzw. 2 nähern. Die Vektoren 


E, — E, = AE, H, — H, = AH (3.1) 

stellen dann die Sprünge der elektromagnetischen Feldstärken an der Sprungflüche F* 

und daher auch an der ihr entsprechenden in der Richtung ihrer Normalen sich 
bewegenden Fläche Ф*(т) dar. 

Ist eine solche Sprungfläche F* innerhalb des Raum-Zeit-Gebietes G gelegen, 

so kann man in (2.8) bei der Anwendung des Gaußschen Satzes nicht über F* hin- 


wegintegrieren und muß daher zur Energiebilanz (2.18) noch das Integral 


Ti 


= | dr fic +), — (— S, + "Рас (3.2) 


to D*(r) 
hinzuzufügen. Dabei wurde die in 82 mit v bezeichnete Geschwindigkeit durch v* 
ersetzt. Die Normalenrichtung n wurde in (3.2) entgegengesetzt gleich der Projektion 
von N also gleich der der Seite 1 der Fläche F entsprechenden Normalenrichtung 
gewählt. Die Ausdrücke (—S, + v*W), und (—S; + v*W), bezeichnen dabei die 
Grenzwerte des Ausdruckes —S, + v*W wenn wir uns der Flache F* von der Seite 1 
bzw. 2 nähern. 


Um die Differenz der in (3.2) auftretenden Grenzwerte zu berechnen, bemerken 
wir daf nach (2.7) 
(Sr + 0*W) = — 7 n(( E + AE) x (Н + AH) 
l > > > > 
+= (e + AE)? + и + AB) (3.3) 


= (—S, + »*W), — JE + АЕ?) — J,(H + AH), 
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wo 
FE an cH Ab (3.4a) 
Ast > > > и > 
EX A nc А (8.45) 
с т, 
bezeichnet. 


Die Energiebilanz (2.18) ist somit nach (3.2) und (3.3) durch den Beitrag der 
Sprungfläche F* 


= jl dr f (ЕЈ + H' J;) do (3.5) 


To Ф”(т) 


zu ergánzen. Hier bezeichnen 
E LI AEP und HHL АН? 


die Mittelwerte der zu beiden Seiten der Sprungfläche Ф*(т) auftretenden Grenzwerte 
der elektromagnetischen Feldstärken. 


Es ist klar, daß die Größen 7 und if. als flachenhafte, elektrische bzw. »magneti- 
sche« Ströme angesehen werden können, die in der in der Richtung ihrer Normalen 
sich bewegenden Sprungfläche ®*(r) fließen. Der Ausdruck (3.5) kann dann nämlich 
als die Arbeit interpretiert werden, die an diesen Strómen seitens des elektromagne- 
tischen Feldes geleistet wird. 

Falls die Stróme jf und J? nicht verschwinden, liefert der Ausdruck (3.5) einen 
im allgemeinen nicht verschwindenden Beitrag zur Energiebilanz (2.18). Das Ver- 
schwinden der Ströme (3.4a und b) stellt daher eine hinreichende Bedingung dafür 
dar, daß in der Sprungfläche F* und somit auch in der Fläche Ф*(т) keine Energie- 
quellen vorhanden sind. 

Nehmen wir an, daß keine elektrischen oder magnetischen Flächenströme in 
einer in der Richtung ihrer Normalen fortschreitenden Sprungfläche ®*(r) fließen, 
so müssen wir voraussetzen, daß nach (3.4a und b) die Sprünge der elektromagneti- 
schen Feldstärken in Ф*(т) die Bedingungen ; 


AH x n= — AE, (3.62) 
m RR 
пх AE = £ AH (3.6b) 


erfüllen. 
Aus diesen beiden Bedingungen kénnen wir sogleich gewisse Folgerungen ziehen. 


Aus (3.62) und (3.6b) folgt zunächst durch skalare Multiplikation mit n 
>. n4É-0, n4H- OQ. (3.7) 
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Multiplizieren wir hingegen (3.6a) vektoriell mit n, so erhalten wir mit Rücksicht auf 


(3.7): n x (AH x n) = 


und daher mit Hilfe von (3.6b) 


(1-4) ado (3.8) 


Eine analoge Beziehung kann man für den Sprung AE der elektrischen Feld- 


stárke ableiten. Damit also im Falle verschwindender Stróme j. und 18 Sprünge 
elektromagnetischer Feldstärken auftreten kónnen, muf) somit die Beziehung 


mt^ enu (3.9) 
bestehen. 
Dies bedeutet: Eine elektromagnetische Sprungfläche in der keine elektrischen 


oder magnetischen Stróme ], und 7А fließen, muß sich mit der Lichtgeschwindig- 
keit v, (2.20) bewegen, mit der sich ebene elektromagnetische Wellen in einem nicht- 
leitenden, homogenen Medium mit den elektromagnetischen Materialkonstaten & 
und u ausbreiten. Die entsprechende Fläche F* muß also eine charakteristische 
Fläche sein. 

Die Beziehung (3.7) besagt dabei, daß in einer solchen mit der Lichtgeschwindig- 
keit v; (2.20) fortschreitenden Sprungfläche Ф*(т) nur die Transversalkomponenten 
der elektromagnetischen Feldstärken von Null verschiedene Sprünge aufweisen 
dürfen. 

Im Falle (3.9) nehmen die Beziehungen (3,6a und b) die Gestalt an 


AB х 5а 
(3.10) 
#х abe ғай. 


Werden in (3. 10) die Sprungwerte AE bzw. AH durch die elektromagnetischen 
Feldstárken E bzw. H ersetzt, so ergeben sich zwei Beziehungen, die von den elektro- 
magnetischen Feldstärken in ebenen elektromagnetischen Wellen erfüllt werden, 
die in einem homogenen nichtleitenden Medium mit den elektromagnetischen Material- 
konstanten & und y sich ausbreiten. Man muß jedoch beachten, daß die Beziehungen 
(3.10) auch in nichthomogenen Halbleitern für jede sich mit der Lichtgeschwindig- 
keit v; (2.20) in der Richtung ihrer Normalen bewegende Fläche Ф%(т) verwendbar 
sind. 

Bemerkenswert erscheint die Tatsache, daB. im Falle verschwindender Ströme /; 
und 5 in dem in den Sprungflüchen keine Energie, also auch keine Masse enthalten 


ist, die Bewegung dieser Flächen sich schon im wesentlichen aus den Maxwellschen 
Gleichungen ergibt. 
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8 4. Der Eindeutigkeitsbeweis für das Anfangswertproblem 


Es seien E,, H, und E,, H, zwei elektromagnetische Felder, die in dem in § 2 
definierten Raum-Zeit-Gebeit G; den inhomogenen Maxwellschen Gleichungen 


QE. Am, + А - 
Su du СЕ 1J) — rot div eE = 4ro, 

E Mies (41) 
в = — rot E, div uH = 0 


genügen. 0 und j sind gegebene Raumdichten der elektrischen Ladungen und Ströme. 
Beide Felder sollen zur Zeit ту in dem Raume R, die gleichen Anfangsbedingungen 
erfüllen, d. h. ihre Feldstärken die gleichen Werte annehmen. Wir setzen schließlich 
voraus, daß beide elektromagnetischen Felder in dem Raum-Zeit-Gebiet C, stetig, 
sind bis auf 

a) im dreidimensionalen Raume vorgegebene unbewegliche Flächen, in denen 


die Tangentialkomponenten von Ё, H, und rs Н, gegebene Sprünge erleiden, 

b) in G, befindliche Sprungflächen F* in denen die Sprungwerte der beiden 
Felder den Bedingungen (3. 4a und b) genügen, wobei längs der ganzen in G, ent- 
haltenen Fläche F* die Stromdichten J, und у. gegebene Werte annehmen. Die 
F* im dreidimensionalen Raume entsprechende Fläche Ф% (т) bewegt sich dabei 
mit der Lichtgeschwindigkeit о, (2.20) oder mit irgend einer anderen Geschwindigkeit 


v* — c[m*. Insbesondere kann dabei оош werden, daß ie = у; = 0 йз 
falls F* eine charakteristische Fläche ist. 

Überdies setzen wir voraus, daß die F entsprechende Fläche Ф (т) sich mit einer 
nicht kleineren Geschwindigkeit als die Lichtgeschwindigkeit v, bewegt, also die 
Ungleichung (2.19) erfüllt ist. 

Um zu zeigen, daß diese Bedingungen das elektromagnetische Feld in dem 
ganzen Raum-Zeit-Gebiet G, eindeutig bestimmen, betrachten wir das durch die 
Differenz dieser beiden Felder gegebene elektromagnetische Feld 

Е= Е.Е, Н= № — №. 
Dieses „Differenzfeld“ E, H genügt dann den homogenen Maxwellschen Gleichungen 
(2.1) und verschwindet zur Zeit то in dem Raume Rọ. Dieses Feld ist in dem ganzen 
Raum-Zeit-Gebiet G, stetig bis auf 

a) Flächen, die im dreidimensionalen Raume vorgegeben und unbeweglich 
sind, in denen insbesondere die Tangentialkomponenten der elektromagnetischen 
Feldstärken stetig sind, 

b) gegebene, sich in der Richtung ihrer АЕ bewegende Sprungflächen 
F*, in denen die dem Differenzfeld entsprechenden Stróme Үр und T. verschwindenÿ. 


УД 
3 Falls F* keine charakteristische Flüche ist, folgt aus (3.8) und der analogen Beziehung, in der AE 


statt AH auftritt, dag das Differenzfeld in F* stetig ist. 
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Е + ә ыыы 
Nach den in § 3 durchgeführten Betrachtungen gilt unter den obigen Voraus- 
setzungen für das Differenzfeld im Raum-Zeit-Gebiet С, die Energiebilanz (2.18), 
trotz der oben unter a) und b) angeführten, in G, enthaltenen Unstetigkeitsflächen 


der ursprünglich gegebenen Felder D H, und fo. Hp. 

Da aber zur Zeit v das Differenzfeld verschwindet, so verschwindet auch JW, , 
so daB in (2.18) eine verschwindende Summe von drei Integralen auftritt. Da diese 
Integrale, mit Rücksicht darauf, daß die Bedingung (2.19) erfüllt ist, alle nicht 
negativ sind, müssen sie auch einzeln verschwinden. Aus dem Verschwinden des 
Integrals über R, schließen wir auf das Verschwinden von |, und daher auch auf 
das Verschwinden des Differenzfeldes іп R, zur Zeit ті. Da wir für 7, in (2.18) auch 
jeden Zeitmoment т, der der Ungleichung v, St 5 т, genügt, benutzen können, 
so ist damit die Eindeutigkeit in dem ganzen Raum-Zeit-Gebiet G gewährleistet. 

Zur Durchführung des obigen Eindeutigkeitsbeweises war die Annahme uner- 
läßlich, daß das dritte in der Energiebilanz (2.18) des Differenzfeldes auftretende 
Integral einen nicht negativen Wert besitzt. Dies ist aber nur dann der Fall, wenn 
die Bedingung (2.19) erfüllt, d.h. die Fläche Ф (т) in der Richtung ihrer Normalen 
sich mit einer Geschwindigkeit v = c/m bewegt, die mindestens gleich der Licht- 
geschwindigkeit о; (2.20) ist. Physikalisch bedeutet dies, daß die Störungen, die 
von der Zeit ту an, von außen her in den Raum В, eindringen, sich mit keiner grós- 
seren Geschwindigkeit als der Lichtgeschwindigkeit v; fortpflanzen. Die Geschwin- 
digkeit v; ist somit, in Übereinstimmung mit den in $ 3 durchgeführten Uberlegun- 
gen über die Fortpflanzung von Sprungflüchen, als die Frontgeschwindigkeit der elektro- 
magnetischen Wellenbewegung anzusehen. Von diesem Gesichtspunkt aus ist auch 
die Tatsache zu verstehen, daß der Eindeutigkeitsbeweis auch noch durchführbar ist 
falls die Fläche Ф (т) sich mit einer größeren Geschwindigkeit als die Lichtgeschwin- 
digkeit v; bewegt. Eine solche Fläche schlieft ja stets nur einen Teil des dreidimen- 
sionalen Raumes ein, der von der mit der Geschwindigkeit v, sich bewegenden Fläche 
® (1) begrenzt wird. Man kann sich daher auf den Standpunkt stellen, daß der Ein- 
deutigkeitsbeweis von Cohn (1900) und Riemann-Weber (1900), bei dem der Anfangs- 
zustand in dem ganzen unendlichen Raume vorgegeben wird, als Grenzfall v — oo in 
unseren Überlegungen mit enthalten ist. 

Es sei noch bemerkt, daß der obige Eindeutigkeitsbeweis sich auch noch in 
dem Falle durchführen läßt, wo in dem Raumgebiet R, vollkommene Leiter anwe- 
send sind (an deren Oberflüche S, = 0 ist) oder auch der Anfangszustand nicht in 
der Ebene т = t, sondern auf einer beliebigen raumartigen Fläche vorgegeben ist. 
Die Bezeichnung ,,raumartig“ ist dabei nicht mit Rücksicht auf die konstante Vakuum- 
Lichtgeschwindigkeit c sondern in Hinblick auf die im allgemeinen variable Licht- 
geschwindigkeit v, (2.20) zu verstehen. 


§ 5. Ableitung der Bedingungen (3.4) aus den Maxwellschen Gleichungen 


Für alle in der Physik auftretenden partiellen Differentialgleichungen vom 
hyperbolischen Typ lassen sich Bedingungen ableiten für die Sprungwerte, die auf 


addu 2 rane 
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den sich bewegenden Flüchen Ф*(т) auftreten. Diese Bedingungen sind auch sicher 
lich für die Maxwellschen Gleichungen bekannt, wenn auch die Suche nach ihnen 
in der mir zur Verfügung stehenden Literatur und die Nachfrage bei Fachkollegen 
leider von einem negativen Erfolg begleitet war. Es soll daher hier in aller Kürze 
gezeigt werden, daß sie mit den Bedingungen (3.4a und b) übereinstimmen und daß 
die Deutung der Ausdrücke (3.6a und b) als elektrische bzw. magnetische Stróme 
auch vom Standpunkte der Maxwellschen Gleichungen gerechtfertigt ist. 

Zunüchst nehmen wir an, daf) der Flüche F* im dreidimensionalen Raume 
eine mit einer beliebigen Geschwindigkeit v* — c/m* in der Richtung ihrer Normalen 
fortschreitende Fläche Ф*(т) entspricht. Eine längs dieser Fläche F* in dem Raum- 
-Zeit-Gebiet auftretende Unstetigkeit kann als ein Grenzfall einer Folge von stetigen 
Fällen angesehen werden. Man kann dann z. B. die x — Komponente der ersten 
Maxwellschen Gleichung (4.1) d. h. 

Wy MATE c с ӘН, ӘН, 
er 1 с (obs = y — 22 


über ein vierdimensionales Raum-Zeit-Gebiet integrieren. Nach dem vierdimensio 


nalen Gaußschen Satz erhält man dann 


(ек. соз (N,t)dF + E foz + jx)dv 
G 


F 


= [a cos (N,y) — Н, cos (N,z)) dF. 
F 


Wir ziehen nun das Raum-Zeit-Gebiet G auf ein in der Sprungflüche F* liegendes 
Flächenelement dF zusammen. Dabei nehmen wir an, daß die Feldstärkekomponen- 
ten zwar endlich bleiben, aber in dF sich sprunghaft ändern und daß die Raumdichte / 
des elektrischen Stromes derart unendlich wird, daß auf dF ein Strom von einer 
endlichen Flüchendichte entsteht. Benutzen wir die gleichen Bezeichnungen wie in 
53, so erhalten wir bei diesem Grenzübergang : 


eAE, dF cos (М, т) + = lim [р dv 
G 

—AH, dF cos (N,y) — AH, dF cos (№). 
Dabei ist zu beachten, daß im Grenzfalle das Raumintegral des Stromes c E, keinen 
Beitrag zu (5.1) ergibt, da ø konstant ist und E, bei diesem Grenzübergang endlich 
bleibt und daß ferner der Grenzwert des Raumintegrals mit j, die x-Komponente des 
gesamten, in das dreidimensionale Flüchenelement dF zusammengedrängten Stromes 
liefert. Dividieren wir (5.1) durch m* dF cos (М, т), so erhalten wir mit Hilfe von 
(2.15) die Beziehung 


(5.1) 


ы. 4 Е Sui Jax = АН, cos (n, =) — AH; cos (n, у). (5.2) 
m c ; 
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Dabei bedeutet J,, mit Rücksicht auf (2.16) 


1 | х 3 
= i Lodo sce] , dv. 5.3 
Jex m* dF cos (N, т) ues fis % ас ат HE fi с) 
G G 


Die Beziehung (5.2) ist mit der x-Komponente der Vektorgleichung (3.4a) identisch. 
Ebenso folgert man aus div е E = 4 ло (vgl. Gl. (4.1)) die Beziehung 


nA ёЁ = — 4n P, (5.4) 
wo P, ühnlich wie (5.3), nämlich durch 
I 
= i 5.5 
ІР; ТУ: lim fe dv (5.5) 
G 


definiert ist und offenbar als eine flächenhafte wahre elektrische Ladungsdichte 
gedeutet werden kann. 

Um die Beziehung (3.4b) zu erhalten, mu8 man auf der linken Seite der zweiten 
Maxwellschen Gleichung (4.1) zu dem magnetischen Verschiebungsstrom noch ein 


Glied mit einem „magnetischen Leitungsstrom“ у, hinzufügen und den bei der ersten 
Maxwellschen Gleichung angegebenen Grenzübergang durchführen. 

Eine zur Gleichung (5.4) analoge Beziehung für den Sprung der Normalkom- 
ponenten von AuH erhält man, wenn man div uH gleich einer 4z-fachen „wahren 
magnetischen Ladung“ setzt. Auf diese Weise ergibt sich 


nA ull = — АарР,, (5.6) 


wo P, analog wie in (5.5) definiert ist. 
Die Maxwellschen Gleichungen ergeben somit die aus dem Energiesatz erhaltenen 
Beziehungen (3.4) für die Sprungwerte der elektromagnetischen Feldstärken und 


präzisieren die in $ 3 angeführte Deutung der Stróme jh und 1 als elektrische bzw. 
magnetische Leitungsstróme. AuBerdem ergeben sich aber noch die beiden Bezie- 
hungen (5.4) und (5.6), die als Definitionsgleichungen für die flachenhaften Ladungs- 
dichten P, und P,, aufzufassen sind, ebenso wie wir (3.4) als die Definitionsgleichungen 
für die flächenhaften Stromdichten J und УЫ ansehen Кбппеп. In mathematischer 
Hinsicht enthalten die Relationen (5.4) und (5. 6) nichts Neues. Falls wir aber die 


beiden Gleichungen (3.4a und b) skalar mit n multiplizieren und mit den Gleichun- 
gen (5.4) und (5.6) vergleichen, so erhalten wir die beiden Beziehungen 


v Р, = п], v Pies nj, (5.7) 
die die physikalische Aussage enthalten, daß die Normalkomponenten der Ströme 


үр und. Ts durch die mit der Geschwindigkeit v* = c/m* erfolgende Bewegung der 
Ek ken bzw. magnetischen Ladungen P, bzw. P,, entstehen, also Konvektions- 
stróme darstellen. Dies ist offenbar eine Folge der отаны А daß die Sprung- 
fläche Ф* (т) sich in der Richtung ihrer Normalen bewegt. 
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Die Bedingungen (3.4a und b), (5.4) und (5.6) für die Sprungwerte elektromagne- 
tischer Feldstärken an sich bewegenden Sprungflüchen kónnen auch aus den entspre- 
chenden Bedingungen an ruhenden Sprungflächen erhalten werden, wenn man die 
elektromagnetischen Feldstärken einer Lorentz Transformation unterwirft. Allerdings 
muß man dabei voraussetzen, daß die Geschwindigkeit v* mit der sich die Unstetig- 
keitsfläche bewegt, die Vakuumlichtgeschwindigkeit nicht überschreitet. Bezüglich 
der Transformartionsformeln vergleiche man etwa Paulis bekannten Enzyklopädie- 
artikel über die Relativitätstheorie (Pauli 1920, S. 655). 

Die Tatsache, daß wir die Bedingungen (3.4a und b) so einfach aus den Max- 
wellschen Gleichungen ableiten konnten, ist darauf zurückzuführen, daß wir die 
Überlegungen in einem vierdimensionalen Raum-Zeit-Gebiet durchgeführt haben. 
Ebenso kann man selbstverständlich auch bei anderen partiellen Differentialglei- 
chungen vom hyperbolischen Typus bei der Ableitung der „dynamischen“ Bedin- 
gungen für die Fortpflanzung von Sprungflächen vorgehen. Man kann z. B. leicht 
zeigen, daß sich im Falle der Elastizitätstheorie die richtigen dynamischen Bedin- 
gungen ergeben. 


8 6. Unstetige elektromagnetische Signale in Wellenleitern 


Das einfachste Beispiel für die Anwendung der obigen Kriterien stellt eine ebene 
elektromagnetische Welle in einem homogenen Nichtleiter dar. Mit den Bedingungen 
(3.10) ist der Fall vereinbar, daß die elektromagnetischen Transversalkomponenten, 
aus denen eine solche Welle besteht, irgendwelche Sprünge erleiden, da ja diese 
Feldkomponenten selbst und daher auch ihre Sprünge die Beziehungen (3.10) er- 
füllen. | 

Ein etwas komplizierteres Beispiel stellt die Fortpflanzung von unstetigen elektro- 
magnetischen Signalen in einem Wellenleiter dar. Im Falle von Wellen vom trans- 
versalen elektrischen Typus kann man Lósungen der Maxwellschen Gleichungen 
angeben (Rubinowicz 1954), bei denen die magnetische Longitudinalkomponente 
(H,) in der Wellenfront einen Sprung erleidet. Da dies das Auftreten von „magneti- 
schen Ladungen“ in der Wellenfront bedeuten würde, muß die in der Lösung der 
Maxwellschen Gleichungen verfügbare Phasenkonstante a so gewählt werden (a — 
= 1/2), daß in der Wellenfront die magnetische Longitudinalkomponente Н, ver- 
schwindet. Setzt man die in der angeführten Arbeit für die elektromagnetischen 
Feldstärken angegebenen Ausdrücke (2.8) in die Beziehungen (3.10) der vorliegenden 
Arbeit ein, so kónnen sie nur dann erfüllt werden, falls die a. a. O. auftretenden Ablei- 
tungen F,(z, t) und F,(z, t) in der Wellenfront der Relation Р, = — F,/v genügen. 
Nach der a. a. O. angegebenen Formel (4.4) ist dies aber tatsüchlich der Fall, 


КРАТКОЕ СОДЕРЖАНИЕ 


B. Рубинович, Распространение скачков напряженностей әлектромаінитноһо 
поля и доказательство однозначности для проблем с начальными условиями 
уравнений Максвелла. ue 

Скачки напряженностей электромагнитного поля появляющиеся на поверх- 
ности, движущиеся в направлении нормали с произвольной скоростью, участвуют 
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B балансе энергии четырехмерной области пространства-времени в качестве 
положительных и отрицательных источников энергии. Физической причиной 
этого обстоятельства можно считать электрические и ,,магнитные” токи теку- 
шие в поверхности разрыва. При этом уравнения Максвелла показывают, что 
эти токи возникают путем вдавливания в поверхность разрыва простран- 
ственных токов проводимости. Эти токи должны исчезать, чтобы скачки не 
участвовали в балансе энергии электромагнитного поля. Из сказанного вытекают 
условия для скачков напряженностей электромагнитного поля, допустимых 
с точки зрения физики. Эти условия, между прочим требуют, чтобы эти скачки 
распространялись только вдоль характеристических поверхностей. 

Также при наличии скачков напряженностей электромагнитного поля, кото- 
ры выполняют данные условия, можно провести доказательство однозначности 
для начальной проблемы уравнений Максвелла, учитивая конечную скорость 
распространения электромагнитных волн. 
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DER SATZ VON DER ERHALTUNG DES IMPULSES UND DIE 
FORTPFLANZUNG VON SPRUNGEN ELEKTROMAGNETISCHER 
FELDSTARKEN 


Von А. Rugınowicz 
Physikalisches Institut der Polnischen Akademie der Wissenschaften, Warschau 
(Eingegangen am 29. September 1954) 


Die Forderung, daf die mit Lichtgeschwindigkeit sich bewegenden Sprungflächen 
elektromagnetischer Feldstärken keinen Beitrag zur Impulsbilanz eines Raum-Zeit- 
Gebietes liefern, ist nur mit dem Minkowskischen nicht aber mit dem Abrahamschen 
Ansatz für den Energie-Impulstensor vereinbar. 


8 1. Problemstellung 


Neuerdings ist ein gesteigertes Interesse für die Frage nach der Gestalt des 
elektromagnetischen Energie-Impuls-Tensors festzustellen, sowohl im Falle isotroper 
als auch anisotroper Medien. Insbesondere hat eine Arbeit von M. v. Laue (1950, 
1952) eine lebhafte Diskussion! hervorgerufen. V. Laue hat sich nämlich gegen den 
bisher bevorzugten Abrahamschen Ansatz für den elektromagnetischen Energie- 
-Impulstensor ausgesprochen und sich für den Minkowskischen Ansatz entschieden. 
Wir wollen hier nicht näher auf diese sehr interessanten Erórterungen über den 
elektromagnetischen und den gesamten Energie-Impulstensor eingehen, die zum 
Teil auch Fragen von grundsätzlicher Bedeutung zur Sprache brachten. Das Ziel 
der vorliegenden Note ist vielmehr, nur einen Beitrag zur Frage zu liefern, für 
welchen von den beiden Ansätzen man sich im Falle isotroper Medien zu entscheiden 
hat, falls man vor die Wahl zwischen dem Abrahamschen und dem Minkowskischen 
Ansatz gestellt wird. Als Kriterium wollen wir dabei den Beitrag von elektromagne- 
tischen Sprungflächen zur Impulsbilanz eines Raum-Zeit-Gebietes verwenden. 

In der vorangehenden Arbeit? wurde gezeigt, daß sobald in einem Raum-Zeit- 
-Gebiet С bewegte Sprungflüchen Ф*(т) elektromagnetischer Feldstärken auftreten, 


1 Wir verweisen hier auf die Arbeiten von Balazs (1953), Beck (1953), Marx (1953), Marx und 
Gyórgyi (1954) und Ott (1953). 

2 In der vorliegenden Arbeit benutzen wir die Bezeichnungen der vorangehenden Arbeit und 
kennzeichnen die dort angegebenen Formeln durch ein vorgesetztes A. 
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sie an der Energiebilanz von G beteiligt sind. Diese Tatsache wird dadurch bedingt, 
daß im allgemeinen i in bew egten Sprungflüchen elektrische und magnetische flächen- 


hafte Stromdichten I, und J. auftreten. Interessant ist speziell der Fall, wo diese 
Stromdichten verschwinden und daher die Sprünge elektromagnetischer Feldstärken 
in Ф* (т) die Bedingungen (A, 3.10) erfüllen. In diesem Falle nimmt die Ф* (т) 
im Raum-Zeit-Gebiet G entsprechende Fläche F* an der Energiebilanz dieses 
Gebietes nicht teil und die Sprungfläche Ф* (т) bewegt sich mit der Lichtgeschwin- 
digkeit v, (A, 2.20). 

Analoge Überlegungen, wie wir sie in A durchgeführt haben, kann man auch 
im Falle anderer Erhaltungssätze anstellen. Wir werden uns in der vorliegenden 
Arbeit mit dem Impuls-Erhaltungssatz beschäftigen. Wir wollen zeigen, daß nur im 
Falle des Minkowskischen Ansatzes für den Energie-Impulstensor eine elektromagne- 
tische Sprungfläche, in der die Stromdichte y und » verschwinden, keinen Beitrag 
zur Impulsbilanz eines Raum-Zeit-Gebietes С liefert. Dieses Argument spricht ein- 
deutig gegen den Abrahamschen Ansatz. 


$ 2. Der Satz von der Erhaltung des Impulses für stetige elektromagnetische Felder 


Um sowohl den Minkowskischen als auch den Abrahamschen Ansatz zu umfassen, 
nehmen wir an, daß der Energie-Impulstensor durch das Schema 


| | м — 
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Wird die Beziehung (2.1) durch die Vakuum-Lichtgeschwindigkeit c dividiert 
und über das in A, $2 definierte Raum-Zeit-Gebiet С integriert, so erhalten wir 
mit Hilfe des Gaußschen Satzes und bei Benutzung der in A, $2 angegebenen Formeln 
für die Impulsbilanz in G mit Rücksicht auf v — c/m die Beziehung 


E , 4х dy dz — ШЕН) "ngu. fines s)as 


Ts D(t) 


1 E 2.3 
-1 са 29) 
с 
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Dabei setzen wir zunächst voraus, daß das in v = c/m auftretende m noch einen 
beliebigen Zahlenwert besitzt, also nicht durch (A,3.9) gegeben wird. Die F ent- 
sprechende Fläche Ф(т) bewegt sich also in der Richtung ihrer Normalen mit einer- 
beliebigen Geschwindigkeit v — c/m. 

Abgesehen vom Vorzeichen, haben die in (2.3) auftretenden Integrale die nach- 
stehenden Bedeutungen: 

Die beiden ersten Integrale auf der linken Seite von (2.3) bezeichnen die zu 
den Zeiten Tọ und т, in den Raumgebieten К, und К, enthaltenen Impulse. 


Der Integrand des dritten Integrals auf der linken Seite von (2.3) setzt sich aus 
"m und 25 5 zusammen. Mit Rücksicht darauf, daß der mit dem negativen Vor- 


zeicher versehene Maxwellsche Spannungstensor Т die Impulsströmung angibt, 
bedeutet das Integral iPS T ndo die in der Zeiteinheit erfolgende Impuls strómung 
in der zur Bewegungsrichtung der Fläche Ф(т) entgegengesetzten Richtung —n. 


Das Integral f. Фа) of S do ist hingegen dem Impulsbetrag gleich, den das durch die 
с 


Fläche Ф(т) eingeschlossene Raumgebiet in der Zeiteinheit dadurch verliert, 4аВ 
die in der Richtung ihrer Normalen sich bewegende Fläche Ф(т) mit der Geschwindig- 


keit v über Raumpunkte hinwegstreicht, in denen eine Impulsdichte > S vorhanden 


ist. Das gesamte dritte in (2.3) auf der linken Seite auftretende Integral gibt somit 
den gesamten Impulsbetrag an, der im Zeitintervall von Tg bis 7, durch die sich 
zusammenziehende Fläche Ф(т) hindurch herausströmt. gen * 

Mit Rücksicht darauf, daß die elektromagnetische Kraftdichte К {=K м oder К) 
als der zeitliche Zuwachs der mechanischen Impulsdichte angesehen werden kann, 
ist das in (2.3) auf der rechten Seite auftretende Integral als der gesamte im Raum- 
-Zeit-Gebiet G infolge der Einwirkung des elektromagnetischen Feldes stattfindende 
Zuwachs des mechanischen Impulses zu deuten. 

Im ganzen können wir somit die Beziehung (2.3) als die Impulsbilanz des Raum- 
Zeit-Gebietes G in der nachstehenden Weise interpretieren: ; 
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Der zur Zeit Tọ im Raumgebiet Rọ enthaltene elektromagnetische Impuls wird 
zur Deckung der nachstehenden Impulsbeträge verwendet: 

a) des Impulsbetrages, der zur Zeit т, in dem Raumgebiet Кі enthalten ist, 

b) des Impulsbetrages, der durch die Flüche F aus dem Raum-Zeit-Gebiet G 
heraustritt, und 

c) des elektromagnetischen Impulsbetrages der im Raum-Zeit-Gebiet С in 
mechanischen Impuls umgewandelt wird. 


$3. Anteil einer elektromagnetischen Sprungfläche an der Impulsbilanz eines 
Raum-Zeit-Gebietes 


Nehmen wir an, daß auf einer mit der Geschwindigkeit v* = c/m* sich bewegenden 
Fläche Ф%(т), die einer im Innern des Raum-Zeit-Gebietes С verlaufenden Fläche F* 


entspricht, Sprünge der elektromagnetischen Feldstürken E und H auftreten. Da 
über die Sprungfläche F* nicht hinwegintegriert werden kann, so ist dann auf der 
-linken Seite der Impulsbilanz (2.3) noch der Impulsbetrag 


-4 ја Е Р S (Mi = nEs) lae (3.1) 


тг) 


hinzuzufügen. Wegen 
Tr=— гг (ERE) — n Е?) Exe = (2H (nH) — nH?) 
erhält man mit Rücksicht auf (A,2.7), (A,3.1), be sowie v* = c[m* die Beziehung 
(ss e 28) = (29 жо 5) ti (i+ 44 (nAE) 
zs) can e os 2) «(ars at 
i). (3.2) 
Führen wir mittels (A,3.4), (A,5.4) und (A,5.6) die flächenhaften elektrischen und 


magnetischen Ladungsdichten P, und Р, sowie die Stromdichten 4» und үй ein, 
so ergibt sich mit Rücksicht auf (3.2) für den Impulsbetrag (3,1) der Ausdruck: 


1 з >; >, ik >. > [t^ € u = 
iios Mod Mid СЫА ES pe 
Lf dh "UP PES ger т PETER EUM ai) 


чи 
quee (2+2 ETP AE) x af | as (3.3) 


Dabei bedeuten E' — E 4- ЛЕ, H = + AHJ2 die mittleren Werte der elektro- 
magnetischen Feldstärken in der Sprungfläche ®*(r). 
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Die einzelnen Ausdrücke, die in dem Integranden des Integrals (3.3) auftreten, 
kónnen als flüchenhafte Kraftdichten gedeutet werden. Dies ist bei den beiden ersten 
Ausdrücken unmittelbar klar. Daß dies auch für den vierten Ausdruck in (3.3) zu- 
trifft, ergibt sich, falls man ihn mit dem Ausdruck j x Ble vergleicht, der die räumliche 
Kraftdichte angibt, die auf die elektrische Stromdichte j seitens eines magnetischen 
Feldes mit der Induktion B ausgeübt wird. Der dritte Ausdruck in (3.3) entspricht 
offenbar dem magnetischen Analogon des vierten Ausdruckes. 

Allerdings sind im Falle des Abrahamschen Ansatzes (p = 1) beim Impulssatz 
die Stromdichten nicht wie beim Energiesatz allein durch 9 und Ja gegeben. Zu Ui 
und Hm treten. nun in den runden Klammern von (3.3) additiv Ausdrücke hinzu, 


die zu AE bzw. AH proportional sind und die ebenfalls als Stromdichten gedeutet 
werden kónnen. Da jedoch diese zusätzlichen Stromdichten in dem Beitrag (A,3.5) 
der Sprungfläche Ф*(т) zur Energiebilanz (A,2.18) nicht auftreten, so muß man sie 
wenigstens in diesem Sinne als wattlose Stromdichten ansehen. 

Im Falle des Minkowskischen Ansatzes (p — eu) für den Energie-Impulstensor 
verschwinden in (3.3) diese wattlosen Stromdichten. Die in (3.3) verbleibenden 


Kraftdichten uf, Es uH und e cE x fe sind dann vollständig der Kraftdichte 
с с 


ix Bic bzw. ihrem magnetischen Gegenstiick analog. 
Die im Falle des Abrahamschen Ansatzes (p= 1) in (3.3) zusätzlich auftretenden 
1 eu—1- l eu —1 
4л 4л т” 
dem Ausdruck 


AE x H entsprechen offenbar 


AR = y o^ hang FJ 
Te на SB Hy Ex By, 
ga 20% 4лс 


3.4 
ot ot (80) 


der nach (2.2) im Abrahamschen Fall zusätzlich zur Minkowskischen Kraftdichte 
K м hinzutritt. 


Verschwindet der Ausdruck (3.3) nicht, so bedeutet dies, daß der zur Zeit To 
in dem Raumgebiet Rọ enthaltene Impuls auch noch zur Deckung des Impulsbe- 
trages (3.3) verwendet werden muß, der auf dem innerhalb G gelegenen Teil der 
Flüche F* enthalten ist. Wir wollen aber hier nicht die Schwierigkeiten besprechen, 


die durch diese Tatsache im Falle nichtverschwindender Flüchenstróme Fa und ‚> 
bedingt werden. 


Wir betrachten nur den Fall p = 7% = 0, wo die Sprungfläche Ф*(т) 
sich mit der Lichtgeschwindigkeit v, (A,2.20) fortbewegt. Da in diesem Falle nach 
(A,5.4) und (A,5.6) auch die Flächenladungen P, und P, verschwinden, so ver- 
schwindet bei Zugrundelegung des Minkowskischen Ansatzes der gesamte Impule- 


betrag (3.3). Die Sprungfläche Ф*(т) hat somit in diesem Falle weder ап der Energie- 
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noch an der Impuls-Bilanz einen Anteil, was vom physikalischen Standpunkte durch- 
aus annehmbar erscheint. 

Schwierigkeiten entstehen aber, wenn wir uns auf den Abrahamschen Stand- 
punkt stellen. Wegen p = 1 verschwindet in diesem Falle nicht der Impulsbetrag (3.3), 


selbst wenn die Flächenstrôme Je ТР und die Flächenladungen P, P, verschwinden. 
In diesem Falle nimmt also die Sprungfliche ®*(r) an der Energiebilanz nicht teil, 
obwohl sie an der Impulsbilanz beteiligt ist. Solchen Füllen begegnen wir in der 
Physik z. B. beim Zusammenstoß eines elastischen Spielballs mit einer starren Wand 
von grosser Masse, wo nur em Impuls- aber kein Energie-Austausch stattfindet. 
Trotzdem fällt es nicht leicht sich vorzustellen, daß die Sprungfläche ®*(r) sich 
so wie ein Spielball oder eine starre Wand in dem betrachteten Beispiel verhält. 
Schwierigkeiten scheint vor allem die Beantwortung der Frage zu bereiten ob in 
der strom- und ladungsfreien Sprungfliche der Impuls nur aufgestappelt oder in 
einen nichtelektromagnetischen z. B. mechanischen Impuls oder den Impuls irgend 
eines anderen Feldes (Mesonen?) umgewandelt wird. 

Man kónnte zur Rechtfertigung des Abrahamschen Ansatzes hóchstens die 
Tatsache anführen, daß der ImpulsfluB im elektromagnetischen Felde auch andere 
schwer zu akzeptierende Merkwürdigkeiten aufweist. So findet z. B. in einem elektro- 
statischen Felde eine Impulsströmung statt (Тп £ 0), obzwar die Impulsdichte 
(wegen Ee 0) verschwindet. 

Will man aber dieses Argument nicht gelten lassen, so wird man mit M. v. Laue 
(1950) den Minkowskischen Ansatz bevorzugen, falls man vor die Wahl zwischen 
dem Minkowskischen und Abrahamschen Tensor gestellt wird. Als Kriterium benützt 
dabei v. Laue die Forderung, daß die mit der Lichtgeschwindigkeit v, identische 
Strahlgeschwindigkeit einer ebenen Welle dem Einsteinschen Additionstheorem 
der Geschwindigkeiten genügen soll. Während aber das hier vorgebrachte Argument 
auch im Falle inhomogener Medien verwendbar ist, ist das v. Lauesche nur auf 
homogene Medien beschrünkt, da nur in einem solchen ebene elektromagnetische 
Wellen auftreten kónnen. 


Entscheidet man sich aber dennoch für den Abrahamschen Ansatz, lehnt aber 
das Vorhandensein von flächenhaften Impulsbeiträgen ab, die in sich mit der Licht- 
geschwindigkeit v, bewegenden Sprungflüchen ®*(r) auftreten, so bleibt wohl nichts 
anderes übrig als anzunehmen, daß solche elektromagnetische Sprungflächen in 
der Natur überhaupt nicht auftreten. Dabei kónnte man etwa darauf hinweisen, daB 
die elektrischen Ladungen und magnetischen Momente in der Natur stets mit einer 
Masse verknüpft sind. Mit Rücksicht auf die Tragheit der letzteren kónnen aber keine 
exakt sprunghaften Anderungen der elektromagnetischen Feldstärken erzeugt werden. 

Aber selbst wenn man diesen Standpunkt einnimmt, muB man zugeben, daB 
elektromagnetische Sprungflüchen, die sich mit der Lichtgeschwindigkeit vj be- 
Mio ein sehr nützlicher Grenzbegriff sind, dessen Verwendbarkeit auBer Frage 
steht. 
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Selbstverständlich verschwindet im Falle der Vakuum-Elektrodynamik der 
Unterschied zwischen den Minkowskischen und dem Abrahamschen Ansatz für den 
Energie-Impulstensor und daher auch die oben besprochene Schwierigkeit. 

Das hier behandelte Beispiel der Maxwellschen Gleichungen läßt vermuten, 
daß man die dynamischen Bedingungen für die Fortpflanzung von Sprungflächen 
wohl auch in anderen Fällen als Prüfstein für die richtige Formulierung von Erhaltungs- 
sätzen wird benutzen können. Im Falle der Elektrodynamik könnte man z. B. noch 


die Rolle der bewegten Sprungflächen im Falle des Impulsmoment-Erhaltungssatzes 
untersuchen. 


‚КРАТКОЕ СОДЕРЖАНИЕ 
B. Рубинович, Теорема о сохранении импульса и распространение скачков 
` направления элекромеанитно1о поля 

Требование, чтобы поверхности разрывов электромагнитных полей движу- 
щиеся со скоростью света, не участвовали в балансе импульса данной области про- 
странства-времени можно выполнить только в предположении, что тензор энер- 
гин-импульса имеет вид. поданный Минковским, а не вид Siege Pome Абра- 
гамом. 
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It is shown that in the relativistic equations of motion for a particle moving in 
external potentials some additional forces of purely relativistic origin appear. These 
forces are repulsive for scalar and pseudoscalar potentials in spite of the apparently 
attractive character of the potential; they prevent the particle from coming too near 
to the singularity of the potential. An adiabatic approximation of the interaction potential 
for two nucleons interacting through the pseudoscalar field has been obtained. An interesting 
feature of this potential is a strong repulsive term acting at short distances between the 
nucleons. | 


Introduction 


It is well known that the adiabatic approach in the meson theory of nuclear 
forces leads to high singularities in the non-central part of ‘the static interaction 
potential. But the non-relativistic Schródinger equation with an attractive potential 
containing a singularity of the type 1/7" (n > 2) gives no reasonable stationary solu- 
tions. This fact was very often understood as an argument against all those meson 
theories which lead to the appearence of high singularities. The proof of the equiv- 
alence of pseudoscalar and pseudovector couplings in the theory of pseudoscalar 
meson field pointed to the importance of relativistic effects in the meson theory of 
nuclear forces. However, the conclusion drawn from this fact that the adiabatic 
approach is in general inadequate seems rather too extreme. Though a great number 
of theoretical physicists are working nowadays with the nonadiabatic theory of nuclear 
forces, the nonadiabatic methods are much more complicated and difficult to handle. 
Thus, any attempt to develop the adiabatic approach and make it more exact seems 
fully justified. In this attempt the notion of the »equivalent« potential, which accounts 
for the main relativistic effects inherent in the relativistic form of the equations of 
motion, may be of great value. | 

One of the most ihteresting problems of the recent theory of nuclear forces is 
the problem of the repulsive forces acting at short distances between nucleons. 
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The existence of a repulsive core was postulated by Jastrow (1951) on the basis of 
some purely phenomenological considerations. Jastrow's assumption was supported 
by Levy (1952) in his non-adiabatic theory of two-nucleon systems. However, no 
satisfactory explanation of this effect was given. It is not clear for instance, whether 
the higher-order-in-g? contributions to the interaction potential cannot destroy the 
repulsive core. The question arises, therefore, whether it is necessary to have an 
evidently repulsive potential in the relativistic equations of motion, or whether the effect 
of repulsion may be hidden in the relativistic form of these equations in spite of 
the apparently attractive character of the potential. In this case the effect of the repul- 
sion could be independent of the sign of the higher-order-in-g? contributions. The 
results of this paper support rather this second alternative. It has been shown that 
the effect of repulsion at small distances follows already from the relativistic form 
of the equations of motion with a scalar or pseudoscalar potential, irrespective of 
its sign and origin. For the sake of simplicity, we restrict ourselves for the present 
to the relativistic equations of a particle in an external potential. The self-field of the 
particle is neglected. 


Equivalent potentials in the classical equations of motion 


Let us consider classical relativistic equations of motion for a particle in an 
external scalar potential ® and an electromagnetic field described by a four-vector 
potential A, = (41,45, Аз, 44 = ig): 


dé, ob = 
BENE] СЕД 
E + gà) A == | 8 9х, == eF ap ds | oe (1) 


m being the rest-mass of the particle, e and g its electric and mesonic charges, £, (s) 
denotes the four-vector of the. position of the particle, c the velocity of light. Further 


TUO EET / v? 
Жолы 22 Gap ЕС 1— = d | (2) 


d£, 


0ш- Vou, UR => MESSE Ха = 2,44 = ich. 
dt ч 


In (1) and the following equations the standard summation convention is employed. 
Greek and Latin indices take respectively the values 1, 2, 3, 4 and 1, 2,3. 
Equations (1) depend in a complicated manner on the velocity of the particle 
and on the potentials Ф and A,. The character of the forces acting on the particle 
cannot be easily understood from (1). The nonrelativistic equations of motion have 
a more obvious and intuitive form. Had we considered the velocity of light c as large 
enough to allow in (1) the transition c — co, we would have obtained too simplified 
equations of motion. Therefore, in the next approximation, the terms linear in v/c have 
been retained. This second approximation is usually called nonrelativistic. The 
procedure outlined above reduces to an expansion of the equations of motion in 
powers of v/c and the neglect of terms proportional to (vjc)” for пр 2. 
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Such an expansion is justifiable only if during the whole motion v <c. The velocity 
of light appears also in some other expressions, e.g. g®/c?. Thus the question arises 
how to treat such terms in which the velocity of the particle does not appear. In 
transition to non-relativistic equations of motion these expressions are often neglected. 
This is however equivalent to an additional restriction |6Ф |< mc?. Both these 
restrictions reduce the domain of applicability of the non-relativistic equations of 
motion particularly strongly in the case of singular potentials; the non-relativistic 
approximation becomes then inapplicable in the neighbourhood of the singularity, 
ie. in that region where the most interesting relativistic effects occur. So we ask 
whether any other expansion leading to the non-relativistic equations of motion may 
be found, so as to describe to a good approximation the motion of the particle near 
the singularity. This turns out to be actually possible, if certain forces of relativistic 
origin are added to the ordinary forces. 

The above method of transition to non-relativistic approximation may easily 
be explained by means of the Hamilton-Jacobi equation (cf. Werle 1953) 


EE occa 
which is strictly equivalent to (1). By putting 
S == mot +S (2,20), (4) 
we reduce (3) to the form | 
Al че а= 
ey? (5) 


, 


Had we been able to neglect m 22 
тс“ ot 


have the form of the non-relativistic Hamilton-Jacobi equation. Thus the non-relativ- 


compared with unity, the equation (5) would 


istic approximation may be obtained, supposing 


1,067 
EENE | d th 6 
"cs otc [Ades (6) 
We assume, therefore, that the function S’ can be expanded in powers of A: 
S = So AS, + A3S4 4. | (7) 


We assume further that from all the operations appearing in (6) only (1/тс?)9 [9 

changes the order of the expression. According to (6) this is equivalent to multiplying 

by À . Thus the order of the expression is increased by unity by this operation. Substitu- 

ting (7) into (5) and comparing the terms of the same order, we get the set of equations 
5 


2 
e кырку + Vs + Vu = 0 (8) 
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where we have put for convenience A"S, = S,. In the lowest approximation we 
obtain the non-relativistic Hamilton-Jacobi equation (8) with the »equivalent« poten- 


tials 


2542 

eee Е: (10) 
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Vn = gb, ° (11) 


We call the potentials (10), (11) »equivalent« because they describe very well all 
the main relativistic effects inherent in the strict equation (3) in spite of the non- 
relativistic form of equation (8). The further equations (9) etc. account for more 
subtle relativistic corrections. 

The above method of solving (3) by means of expansion (7) together with (6) 
is particularly simple and clear in the case of static potentials. Then the energy integral 
exists and we can put 


8-5 а) + S" (a), (12) 


where 
950 _ i қ 
Now ме assume 
| S” = AS; (xj) ات‎ 25, (x) +... (14) 


thus reducing (7) to the expansion in powers of the parameter A = e/mc?. On substitu- 
ting (12) and (14) into (5) and comparing the coefficients of the same powers of A, 
a somewhat simpler set of equations is obtained: 


l " е ~ jn 9S, 

i [уже а) р Кама nee (15) 
Il = СҮР 1 p 1 9S, & 

ne (vs -< 3) E ee Ber (16) 


From the form of expansion (14) it may easily be seen that the non-relativistic equation 
(15) with the equivalent potentials (10) nad (11) describes the better the motion of the 
particle in static fields A, and ® the smaller [| is. In the limit when £—0, 
we get from (15) exactly the same orbits as from the exact equation (3). Thus the 
smallness of the non-relativistic energy le] compared with the rest energy mc? forms 
a condition for the equation (15) to be a good approximation of (3). On th» other hand 
the velocity of the particle can reach in some regions values near c. The magnitude 
of g Фс? is also unrestricted. In Figs. 1 and 2 a comparison is given of the region 


of applicability of the usual non-relativistic equation АБ = —g grad Ф for the 
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particle moving in a central static scalar potential @ — Ф(г) and of the equation 


тё = = — grad Г„ with the equivalent potential (11). In Fig. 1 the cross-lined area 
denotes the range of energies E and distances r at which both the conditions v «c, 


Е 4 ғ 
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Fig. 4 


Fig. 3 


= << m are fulfilled. This area does not reach the centre of forces. The applica- 


bility of the equation PET = — grad V„ includes on the contrary the whole region 
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of the centre. By means of the equivalent potentials the motion of the particle can 
be described correctly in the static case provided E | eco 
Thus we see that the motion of the particle described by the relativistic equation (1) 
may to a good approximation be also described by means of the non-relativistic 
equation if one adds to the standard Lorentz and meson forces two forces of relativistic 
origin 
Ж 6% — =? 


Pr om grad ф?, Ё, =з grad Ф?. (17) 


The force F, tends to draw the particle irrespective of the sign of the charge in the 
direction of increasing ГА values. The second force F, connected with the meson 
field on the contrary pushes the particle out of the area of increasing |Ф| values. 
These forces are very important in the neighbourhood of the singularities where (1 ог 
| Ф | becomes very large. In Figs З and 4 the shape of the equivalent potentials V, 
and V, for the case of central attractive static fields p and Ф containing a singularity 


2 
of the type — ue (n > 1) is represented. 
7” 


Let us examine the difference between the orbits obtained from the exact equa- 
tion (3) and the approximate equation (15). As example we shall take the motion of 
a particle with charge e in the Coulomb field y = a ‚ A, =  — 0. Then 

r 


Veri oe aici ef = ее» (18) 


Denoting the angular momentum about the centre by J we get from the exact equa- 
tion (3) the following expressions for the orbits of the particle in polar coordinates г, 9: 


a) for the case Jc > 8: 


l = mc? À 2 2 1 E 7 
r Bag p epe == i) cos [У T 


(19) 
b) for the case J?c? < 2; 
1 Р 2 Я x ^ 
ras (sign p) ие + 2A J?e? ( +2) coh |o a — | 
mel + À JB 
md : (20) 


From the approximate equation (15) we obtain: 


a) for the case J2c2 > 5: 


1 mc? 3 : 
7 = pags У 51-20) + 24 Pet cos E ie fe — E: (19) 
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b) for the case J?c? < B?: 


ке 2 А 
—= (sign В) a Ra 2 A В? (1 — 22) + 24 Jc? coh |o T A i- La mc T 


(20) 
Comparing (19) and (20) with (19) nad (20°) respectively it is readily seen that the 
differences in the shape of the orbits may be neglected provided р = |= тс? «cus 
Such a good agreement is also found in other cases for the motion in static electric and 
magnetic fields. From (19°) and (207) we see immediately that all the typical relativ- 
istic effects such as the perihelion motion on rosettes and deformed hyperbolas, 
spiral motion and falling on the centre even at J > 0 etc. are perfectly preserved. Had 
we employed in the non-relativistic equation instead of (18) the pure Coulomb poten- 
tial f/r, then all these effects would have been lost and we should have obtained only 
motions on conic curves. 

Now we shall consider the motion in a static scalar field Ф. Let us take as an 
example the Yukawa attractive potential Ф = — g?e "Jr. In the case of a Coulomb 
field the additional relativistic force was attractive. Now, the situation is reversed: 
at small distances a force occurs repelling the particle from the centre. This may be 
readily seen from the form of the equivalent potential 


eur 54 e~ ur 


Ин = — е? 


r Dm S aD 
The distance r from the centre, where the attractive force changes into repulsive 
one can be obtained from the condition dV,,/dr = 0. The stronger the attraction at 
large distances (expressed e.g. by the magnitude of the coupling constant g), the 
stronger is the repulsion and the earlier will it occur as the particle approaches the 
centre. For the Yukawa potential г depends on g and р. For not to high values of g? 
the value of г’ = g?[mc? increases approximately proportionally to 22. Assuming 
g 10е, u = 7.1012 cm”! (e denotes here the elementary electric charge), we. get 
r © 10714 cm. If |e| «&mc? then the minimum distance гу at which the particle can 
approach the centre is about r’/2. The depth of the potential hole described by (21) 
depends weakly on g? and is approximately equal to тс?]2. | 

The occurence of repulsive forces at small distances from the singularities of 
the field is disguised by the complicated form of the relativistic equations. We may 
convince ourselves of the reality of these forces apart from any approximations by 
means of a direct and more exact analysis of (1), i.e. for 4, — 0, from 


met + en = — {6 ein (22) 


It may be shown that when approaching the singularity of Ф monotonically, the 
direction of the acceleration is reversed at the point where the effective mass (m+ g®/c*) 
becomes negative. This is true for the particle which approaches the singularity of 
the field close enough, no matter whether this field is static or non-static. Thus, at 
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small distances, the repulsion must occur between any arbitrarily moving particles 
interacting with one another through the scalar field Ф. It follows from a discussion of 
the exact relativistic equation (22) that the effect of repulsion appears at arbitrarily 
high energies. Therefore, we can try to change the form of the potentials in order 
to apply the non-relativistic equations of motion at arbitrary energies E. For a static 
scalar field an exact expression for the equivalent potential was recently found by 
Marx and Szamosi (1954). In this case the energy integral following from the fourth 


equation (22) exists: 
(23) 
By means of this relation the complicated expression (mc? + g®)/V1 — v?/c? may 


be eliminated from the equations of motion. Thus we see that for the case of a static 
scalar field there exists an equivalent potential of the form 


2mc? 


Vs aes (о * A (24) 


which inserted into the equation 


os (al 
dt? = IX хь-ер (25) 


describes exactly the relativistic motion of the particle. The dependence of the 
potential upon the energy should be noted. When the total energy E = mc?(1 + À) 
increases, the potential decreases, the shape remaining unchanged. For small |A] 
(i.e. small Б! the potential (24) acquires the approximate form (11). 


The falling on the centre in the classical equations of motion 


Let us discuss the motion of a particle in a static field y or ® from the point 
of view of the stability of the motion. We shall assume that the potential р or Ф which 


. . X . Li . Р 2 
attracts the particle is central and contains a singularity of the type — 0 (n2 1). 
т" 


In the case of an attractive potential p(r), which actually represents the fourth compo- 

nent of the vector Æ, we conclude on the basis of (20) that even for n = 1 the particle 

falls on the centre independently of the magnitude of energy if its angular momentum J 

is less than ШЕ For n > 1 the falling on the centre always occurs irrespective of 

the value of angular momentum and energy. This is readily seen from the form of 

b: i of the Hamilton-Jacobi equation (3) for a static central potential p = g(r): 
p= Ф = 0: ' 


"AR 
s=- + 0+ f = (Е—еру#— À — mie dr (26) 
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The motion takes place within the limits To <F & гу, where the expression under 
the root is positive. If п > 1, then in a certain neighbourhood of the point г = Othe 
expression under the root is always positive. It means that the particle always falls 
in its motion on the centre. 


2 
The motion in a scalar field with the same singularity — ut is quite different 
т" 


from that given above. At small distances from the centre the repulsive force always 
prevails over the attractive one. Therefore the effect of falling on the centre does not 
occur at all. This can be seen also from the exact solution of (3) for a particle moving 
in a central static scalar field Ф = ®(r): 


2 2 
s--ne pef наара (27) 


c? 


In a certain neighbourhood of the point r — 0 the expression under the root in (27): 
E 
f( E) e а n (me? go) (28) 


is always negative, so that the particle cannot reach distances smaller than rọ. The 
extreme values ry and r, are found to be the real and positive roots of the equation 


Tir, E) 0: 
All the results obtained here depend solely on the transformation properties 
of the fields. Therefore they are valid for any vector and scalar field. Thus not only 


the type of the singularities but also the transformation properties of the potential 


determine the main features of the motion at small distances from the singularity. 


Singular potentials in the Klein-Gordon equation 


Let us now investigate the question, whether the phenomena discussed above 
have their analogue in the Klein-Gordon equation of the quantum theory: 


д : e +\* o 
[= СЕС + с? (ау à) е (ne +4 о (29) 
Now we put 
wae Y ir 720) 
and asume the existence of an expansion of Ÿ in powers of some small parameter 
1: 
P= QA WAY... . (30 


We assume further that the operation (ih/mc?)9/9t is equivalent to multiplication by A’, 
it increases thus the order by unity. For the sake of brevity we do not intend to repeat 
the reasoning given already for the classical equations and therefore we shall restrict 
ourselves to the static potentials and stationary solutions. Then we can put 
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where A = e[mc?. On substituting (30) and (32) into (29) and comparing terms of 
the same order in À, we get the following set of equations 


2 
IL (av = a) EYE mel =0, (33) 


2 
ы | ihv - à) H Ve + Vin — 1 y= 5 E e0) 0. (34) 
So in the lowest approximation we get the non-relativistic Schrödinger equation (33) 
with the same equivalent potentials (10) and (11). In this approximation we have 
to do with exactly the same additional forces as in the classical equations of motion. 

Let us take, as an example, the motion of a particle with a charge e in the Coulomb 
potential ф = — Ze/r, A, = ® = 0. Then the Schrödinger equation (33) will take 
the form 


2 2 264 
E h A Ze T ^e E y. T: EP, (35) 


2m r 2mc? r? 


which.describes very well all the relativistic effects that are present in (29) provided 
В <mc?. The additional attractive term Z?e*/2mc?r? splits the levels with different 
quantum numbers J of angular momentum, i.e. it removes the known degeneracy which 
exists for a pure Coulomb potential. Additional attraction increases the density 
of probability of finding the particle near the centre. E.g. for [ = 0 the wave function Ч 
contains а singularity of the type r(Vi- 4-1)^ where а = 2Ze?|hc (i.e. approx. г" 
for small Z). For Z large enough the particle falls on the centre and no reasonably 
determined stationary states exist. For Г = 0 the critical value of Z equals hc/2e? = 
137/2. The solutions corresponding to Z > 137/2, 1 = 0 oscilate strongly in the 
neighbourhood of the point r — 0: 


V = 1-2 (e(E) cos A lnr + (Е) sin A In r}. (36) 

Due to this behaviour the eigen-values as well as eigen-functions of the energy cannot 
be found unambiguously. 

If the potential ep contains an attractive singularity of the type — у?/п? (n > 1) 

then for arbitrary values of the angular momentum / no reasonable stationary states 


exist at all. This is readily seen from the approximate equation (35) or more directly 
from the exact equation (29), i.e. 


о 7? : | 
ls ih = TM B — h?c?A + п V = 0. (37) 


The radial part of the stationary solutions of (37) has the following asymptotic form 
for г > 0 (Case 1950) 


War. fa (E) cos <=; + сз (E) sin À (38) 
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Disregarding the phase difference 2/2 which is not essential, both independent solu- 
tions occuring in (38) have a similar strongly oscilating character. Therefore the 
eigen-values and eigen-functions of the energy cannot be distinguished in a way 
reasonable from the physical point of view (Case 1950). This difficulty is very often 
interpreted as resulting solely from the occurrence of a singularity in the attractive 
field. This is not true, as we may well judge from the fact that the additional relativ- 
istic terms appearing in (33) have opposite signs in spite of the same attractive 
singularity in e? and gO. 

Let us examine the asymptotic behaviour of the wave function Yin an attractive 
scalar potential g® containing an attractive singularity of the type —y2/r"(n > 1). 
Equation (29) takes now the form 


2 eus 
Ie ee м (ne — wo (39) 


The asymptotic form of the stationary solutions of (39) is quite different from that 
given in (38): 
A 


т --— 
War le He x 4> 0. (40) 


In order to have quadratically integrable solutions we must require с, (E) = 0. This 
condition together with an appropriate condition at infinity determine the eigen- 
values and eigen-functions of the energy in an unambiguous way. It follows that 
due to the relativistic repulsion the eigen-problem for the motion of the particle 
in a scalar static field is always well defined irrespective of the height of the singularity. 
The repulsion at small distances from the centre manifests itself in a rapid decrease 
of МІР — r"”2e-24”=1, The higher the singularity of the field, the stronger is 
the repulsion at short distances and the quicker the decrease of the probability of 
finding the particle near the centre. 

Thus we see that with the Klein-Gordon equation the same additional repulsive 
or attractive forces appear. The asymptotic behaviour of the wave function near the 
singularity depends not only on the type of the singularity but also on the trans- 
formation properties of the fields. 


The Dirac equation with a scalar and vector potential 


Now we shall consider the Dirac equation for a particle moving in given external 
vector potential À, and scalar potential ®. The effect of relativistic attraction or 
repulsion appears also in the present case. The existence of these effects may be shown 
by iterating the Dirac equation 
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in the conventional manner by means of the operator 


= (% —th A + ca ул + о, (mc? + 0) 3 (42) 


In (41) and in the following the standard notation 01,02» 0з›0, а, for the Dirac matrices 
is used. The second order equation obtained by operating with (42) has an identical 
part with the Klein-Gordon equation (29) and some other terms connected with the 
existence of spin. In these terms the so called Dirac »mixing« matrices 01 and б» 
are still present. They make a direct transformation to the non-relativistic equation 
of motion impossible!. Thereore, we shall make a more consistent transition to the 
non-relativistic form, which does not introduce any additional solutions. After putting 


үл 1 Vs 
к LIRE Ya | отс 
p = А = (+ ве n 3% = W arca d Vs lo) 
Ya 0 0 
we can write equation (41) in the form of the two following equations 
550 Spe ir 
ih x — ep — ЕФ | y = со | -ihV— 4 X: (44) 
| ih Û, + 2met — ep + 06) z= (пу 2) (45) 


Here only usual Pauli spin matrices с, appear. From (44) and (45) we can eliminate y 
` or x, thus obtaining an equation of the second order. In the standard method the 
expressions ep and g® are neglected compared with mc?. In the case of singular 
potentials this procedure is equivalent to the exclusion of the region in the neighbour- 
hood of the singularity. Therefore we must eliminate у retaining all the terms. For 
the sake of simplicity we shall consider separately the motion in a central static 
potential of the vector type eg(r) and then in the scalar potential g@(r). Assuming 
that the- function y can be expanded in powers of À, we obtain for the stationary 
. solutions in the first case 


h2 572 Ta 
ere V, — U: pe i hs = 270 (46) 


L being the orbital momentum operator 


L=rx (—ihy). (47) 
Proceeding similarly with the scalar potential, we obtain 
h? 9 1-5 
назы, ero, (48) 


1 The influence of the relativistic correction appearing in the iterated Dirac equation upon the 


scattering of nucleons by a scalar static potential has been calculated by J. Plebaüski and J. Sawicki 
(unpublished paper). 
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where 


he dp. 1 hg аф 1 
an dr 2022. еф? Um 2m dr ` 2mc? + gD ` ea 


Equations (46) and (48) have the form of the non-relativistic Pauli equations. The 
familiar equivalent potentials V, and V,, occur there. They lead to the relativistic 
attraction or repulsion at small distances from the singularities. Moreover there are 
some other additional terms, e.g. a term containing the spin-orbit coupling. The 
difference in sign of all these relativistic corrections in (46) and (48) is to be noted, 
the structure of individual terms being the same. It is worth while mentioning that 


the aL term in (48) is distinguished not only by the sign before U,, which is opposite 
to the sign standing before U, (for г > го). But this term is also much greater in (48) 
than in (46) due to the structure of the denominator of U,, which equals zero at the 
point at which 2Ф = — 2mc?, contrarily to the denominator in U, which for the case 
of the attractive potential (ep < 0) increases monotonically while approaching the 
centre. 

A detailed analysis of the character -of the forces at small distances from the 
singularities on the basis of (46) and (48) is rather complicated. We can make this 
analysis more straightforward by discussing the asymptotic form of the wave functions 
in the neighbourhood of the singularities. Let us first consider the Dirac equation 
with the Coulomb potential ep = — Ze?}r, А, = Ф = 0. The stationary soluiions 
are well known for this case, and we can directly find how the additional relativistic 
attraction manifests itself. We find that this effect causes an increase in the density of 
probability of finding the particle in the neighbourhood of the centre. For the least 
possible. value of the quantum number of total angular momentum J — 1/2 the 
wave function ¥ contains a singularity of the type тҮ1—#—1 where b = Ze?[hc. It follows 
that for J = 1/2 and Z > hcle? = 137 we have no reasonable stationary solutions, 
because the solutions of the Dirac equation contain a superposition of two. highly 
oscillating functions cos À In r and sin A In г. 

Let us now consider the case in which the function eg(r) or g(r) contains a still | 
higher singularity of the same type —7?[r" (n > 1). After the separation of the angular 
part, the Dirac equation for the radial part of a stationary solution is reduced to the 
set of equations (cf. Case 1950) 


ШЕРІ? D Ж. (50) 
ағ r | ү 

dw ‘w y 

EI A AN -(1+ва? P)» (51) 
dr r 


The upper sign refers to the »vector« potential ep and the lower to the scalar poten- 
tial g®. For the sake of simplicity a system of units is employed in which = 1, c = 1, 
— ]. The quantum ‘number х equals + (J+ 1), the constant B denotes the 
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total energy in units of mc? i.e. B = 1 + À. Putting и = г/у and w = r” fa and 
assuming approximately B = 1, we get 


ЕКІ NIE 3 
dr = у ca Je | dr = |2 T T r Um ( ) 
Eliminating from (52) the function fọ, we obtain 
Де A 
dr? 25 r 4, 4 E zE E А 0. (53) 


Having а solution f, of (53), we can obtain f; directly fron? the first equation (52). 
Therefore the asymptotic behaviour of f} determines that of fọ. If the singular po- 
tential is the fourth component of à vector, then the asymptotic form of f} contains 
a superposition of two strongly oscillating functions 


TEA ҒА 
(Een mi e E)e 7 ni. (54) 
For the scalar potential we find оп the contrary that the asymptotic form of f, contains 


a superposition of the form! 
= Sr 55 
du m oi SS 


This great difference in the asymptotic behaviour is due to the difference in sign 


in (51). It follows from (54) that for an attractive highly singular potential which 
is the fourth component of a vector no reasonable stationary states exist for arbitrary J. 
But for the scalar potential with the same singularity the eigen-value problem is well 
defined because we must require c,(£) = 0. This condition together with an appro- 
priate requirement at infinity determines the eigen-values and eigen-functions un- 
ambiguously. The relativistic repulsion at short distances manifests itself in the 
extremely rapid decrease of the density of probability of finding the particle at 
a distance r from the centre. The higher the singularity, the quicker is the vanishing 
of М ~ eal", 


The Dirac equation with a pseudoscalar potential 


This case is of particular interest, as according to the present state of our know- 
ledge the nucleons interact mainly through the pseudoscalar meson field. The Dirac 
equation for a particle moving in a pseudoscalar potential ® with pseudoscalar 
coupling has the following form 


TL 
Img cies ee me + ge, о-о (56) 


This equation contains the potential Ф’ multiplied by the »mixing« matrix оу. Therefore 
it is often stated that the interaction of nucleons with a pseudoscalar field with pseudo- 
scalar coupling vanishes in the static approximation. This reasoning seems rather 


1) Note added in proof; After finishing of this work I found out some interesting papers by 
P.E. Kunin and J.M. Taksar published in Latwijas PSR Zinatnu Akademijas Vestis 11 (1951). 
which contain some results overlapping with those of mine. 
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superficial and delusive. Thanks to (43) we can write equation (56) in the form of 
the following set of two equations 


PEN E NE 
ai À mec: 74 = ісФіу; (57) 
(in Eod + 2а) ر‎ == == md oD! 
24 X ica y x + igd*y. (58) 
It is even simpler to eliminate X in the present case of a static pseudoscalar potential 
than in the case of a vector or scalar potential. For stationary states we get exactly 


Ec pn EU gyal Е À 
4ا‎ + p ers (1+ 4)» = 


the equation obtained having the form of a non-relativistic Pauli equation. The 
equivalent potential has now the form 
á he > 2 

Vin = эу; o yd +g ФА. (60) 
In (60) the term (g?/2mc?) Ð? accounting for the repulsion is present again. If O 
contains a singularity higher than 1/r, then it is immediately evident that the repulsive 
term will always prevail at short distances. At larger distances the spin term с. vo 
plays of course the decisive part. 

The potential ® can be created by another particle, which we denote by the 
number 2. The particle the motion of which we are just considering shall be denoted 
by the number 1. In the first approximation we neglect the retardation of the meson 
field by taking the well known »instantaneous« potential? 


= hg » et 


2mc r (61) 


On inserting (61) into (60), we obtain the following adiabatic approximation to the 
equivalent interaction potential between two nucleons 


, a? + > ет” a? uà e-?ur 
Vin = 3 (6r 68 S ROT — + 52000 7н: (62) 
The abbreviations used are ; 


3 (d, rigat) =» > Me 
pc AUR AN ۵= (1+2) 
3 3 h un 
E 522 ; Eie 
TE seat e ee отб 
Тһе 1/73 singularity appearing in the first part of the potential (62) is irrespectively 
of its sign not essential, since at small distances the repulsive term containing а 1/r4 


? Here we make use of the equivalence theorem of pseudoscalar and pseudovector couplings in 
the soures of the meson field. Cf. Dyson. (1948) and Case (1949). 
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singularity will always prevail. We thus see that the effect of relativistic repulsion 
changes completely the form of the interaction potential for small r. 

All these results can easily be generalized for the case of an arbitrary mixture 
of charged and neutral meson fields. E.g. for the symmetrical theory, we obtain 


in a similar way 


| a? 2 u, eur g?a? : 2 р: 64) 
V. = № Т(а) Т(:2) {oi Og + 512 В} - г 2 x Tai Таз) |. Q Due ( 
3.4 1 à 2mc = r 
= = 


It is clear that the above method of obtaining the equivalent potential for the 
problem of two nucleons interacting through a pseudoscalar meson field is not the 
best possible one. For a consistent transition to the two-body problem the nucleons 
should be treated from the beginning in a more symmetrical manner. It is probable, 
however, that the form of the equivalent potential (64) exhibits the main features of 
the interaction potential for two-body problem. Of course the effect of retardation 
or higher-order-in-g? processes can considerably change the shape of the repulsive 
term, the effect of relativistic repulsion being rather preserved. In any case the potential 
(64) illustrates well the fact that the repulsive forces occuring at small distances 
between nucleons result already from the relativistic form of the equations of motion. 


Conclusions 


The considerations contained in this paper show that in the relativistic equations 
of motion there appear some additional forces, which play a particularly important 
partein the neighbourhood of the singularities of the potential. Their character may 
be easily understood by passing from the complicated relativistic equations to the 
more familiar and obvious non-relativistic equations. This transition is to be done 
in such a manner as to get rid of all restrictions of the type v «£c, gd ne”, eu. 
which are untenable in the vicinity of the sigularity. A method of expansion of 
relativistic equations has been given, which in the first approximation leads to the 


non-relativistic equations satisfying these requirements. For the static potentials 
the method reduces to an expansion in powers of the parameter À = e/mc?. The forms 
of the equivalent potentials obtained describe very well the motion of the particle 
in the neighbourhood of the singularity provided | à | “1. It has been shown that, 
according to the transformation properties of the potential, in all considered cases 
some additional attractive or repulsive forces appear at small distances from the 
singularity. For the scalar and pseudoscalar potentials these forces are repulsive.The 
reality of this relativistic repulsion has been also shown by studying the asymptotic 
behaviour of the stationary solutions of the exact relativistic equations. 

In the end it is worth while to discuss briefly the importance of these relativistic 
forces. In the atomic theory additional relativistic attractive forces become essential 
at distances 101—105 times smaller than the atomic diameter. Thus they give only 
very small corrections to the results of the common non-relativistic atomic theory, 
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The case of nucleons interacting through a scalar or pseudoscalar potential is quite 
different. Rough estimations lead to the conclusion that the forces of relativistic 
repulsion begin to play an important role at distances not much less than the range 
of the common nuclear forces. Therefore, it is clear that they ought to play a very 
decisive part in the theory of nuclear forces. 


The author wishes to express his thanks to Prof. L. Infeld for many valuable 
suggestions concerning the problem of singular potentials and to Doc. J. Plebañski 
for an interesting discussion on the variational method of obtaining equation (1). 


КРАТКОЕ СОДЕРЖАНИЕ 


Верле, Синулярные потенциалы B релятивистских уравнениях движения 


Показано, что в релятивистских уравнениях движения частицы движу- 
шейся во внешних потенциалах появляются своеобразные добавочные силы чисто 
релятивистского происхождения. Эти силы являются отталкивающими для ска- 
лярного и псевдоскалярного потенциалов несмотря на притягивающий -- каза- 
лось бы -- характер потенциала. Эти силы не дают частице приблизиться слиш- 
ком к сингулярности потенциала. Получено адиабатическое приближение взаимо- 
действия двух нуклонов для случая псевдоскалярного поля. Интересной особен- 
ностью этого потенциала является сильно отталкивающий член действующий на 
малых расстояниях между нуклонами. 
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In a previous paper (Lopuszanski, 1954, hereafter referred to as I) an analytical 
solution of the fluctuation problem was given by means of the G-equations of Jánossy 
for nucleon and electron — photon cascades without ionization loss. The methods 
developed in I will now be applied to the case of the electron-photon cascade taking 
into account the ionization loss. The case of a nucleon cascade with ionization loss can 
be dealt with in the same manner. 

We shall make use of the assumption that the cross sections for bremsstrahlung 
and pair production denoted by w and w(? respectively, are homogeneous functions 
of the energies of the scattering and scattered particles. The collisions are elastic. 
We assume further that the loss of energy f suffered by the electron on a unit pau 
is constant. We have 

B for i = 1 (electron) 
be lo for i — 2 (photon) 


The fundamental стшп for the threshold energy E (Jánossy 1950) are 


ш 
(Е, E) + a, GO (Ey, E) + бук (En E) p 
E, 
fa GOXE', E) G&—»(E, — E' ,E) и® (ғ E , i=1,2 
0 0 
0 


where 


GY (Ey, E, uy, Up, 2) = d DOE,, E, пу, na, x) um шір, 
; 73,12 
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are the generating functions of the problem. For brevity these functions are denoted 
in equations (1) by GO (E, E). Further we have 


1. 
a = | w (e) de 
ù 4 
E, denotes the energy of the primary particle, x — the depth of the absorber. 


The boundary conditions are 
и; for ELS E, 
1 for E> Е, 


We introduce, as in I, the auxiliary function Г @—]— GO, With these functions 
equations (1) take the form 


GH (Eo, E, uy, Ug, 0) = | 


Iro oro 
— (i) 
ОА А en 8E, 
Е\ E 
‘D(H’ (3—1) = (i) 
- fu (Е', E) + '6—(E, — E', Е)) w p Е А 
Е, 
E i Га) (Е', E) T@— (E, — Е, E) wo (2 ) EE 
5 % 0 


Now the methods developed in I enable us to reduce the non-linear equations (2) to 
the infinite set of linear ah 


ш ғуаға 
‚Г? + В, a (TO(E', E) + PG—(E, — E' 0 
+ ai СЕ E ( ) + (Eo — Е’, E) ш E 3T 
Eo 
E se ГОЕ', E) ГЕТЕ, — Е, E) ШЕ: (3) 
0) Е, 


0 


valid for E> E,/n, where Г are the unknown functions and Г , are the solutions 
of our problem for E > E,/(n — 1). We regard the latter as given. 

The solution of (3) can be given in terms of the first moments of the distribution 
and of the cross sections. If we denote the four first moments by S“™(E,, E, x), (i, m= 
— 1,2) the solution of (3) reads 


2 
IK Eo, E, u, ua, 2) = У -S99(E,, E, x) (1 — um) 


m=1 


x 2 E Е, 
11 950») 11 , 1 Е 
NONE Ser (UE а-ә far gr um (5. (4) 
0 0 


pU ES E, Uy, Иә, &) ПОЕ“ == E", E, U], Иә, é) 
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To prove it, we must take into account that S$” satisfy the following equations 


9 S Gm) . 9 S Gm) 
er МЫНА ` а. бат) 2 
ga (Bos E) + 04S (Ey, E) + BF (Вы E) 


(5) 
E, 
= | {SGm(E", E) + SC-im(E, — E', Е)} и® a PE nete 
2 0 0 
with the boundary conditions 

: д. (ур БЕСТЕ 

(іт) = im A9 

ио | for E > Eo 


The expressions for S“” сап be deduced from the results obtained by Bhabha and 
Chakrabarty (1943, 1948) for the differential energy-spectrum. Also the results obtained 
by Messel and Potts (1952) yield immediately the S“” functions. 
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In 1947 F. Friedlander published an interesting contribution to our knowledge 
of the asymptotic behaviour of solutions of Maxwell’s equations (Friedlander 1947). 
He derived a vectorial intensity law (i. e. a law concerning the vectors E and H of the 
electric and magnetic field and not the scalar amplitude u of the scalar theory of light) 
and also a law of the variation of the direction of polarization. Both laws were given 
in the approximation of geometrical optics for inhomogeneous and unmagnetic 
(u = 1) but isotropic and continuous media. "The aim of the present note is to call 
attention to a small but essential mistake which occurred in Friedlander's intensity 
law and its derivation. 

Friedlander's intesity law has the form (1. c. formulae (11) and (11a)) 


€E? dS = const. 


mds- pen along a „tube“ of rays, (1) 


where € (x, y, z) denotes the dielectric constant, E and Н — the electric and magnetic 
vectors respectively, dS — the normal cross-section of the infinitesimal pencil (,,tube“) 
of rays. Formulae (1) are derived from the integral law of Friedlander (his formula 


(10a)) 
li IER 57 dS — 0, (2) 


where f is a phase function satisfying the eiconal equation 
| (grad f)? = e, (3) 
whose characteristics are determined by (т — a parameter) 


dx ӘР dy ӘУ dz of df _ 


Es , 


dv Эл’ у’ ad 


ӘЛ 712- of +, 2 5 
5 ox] 29x’ © dy} | 2 ду” ат \9z 294” 


(255) 


(4) 
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In (2) 9/8 v denotes the outward normal derivative to the closed surface of integration. 
Friedlander takes the integration in (2) over the walls of the ray tube and two normal 
end sections of it. He writes: „It follows, since 9f/9 v = 0 on the walls of the tube and 
290 = + e on the ends (by (3) and (4)), that the quantity &Е?45 is constant along 
the tube“ (we apply to this quotation our notations and numeration of formulae). 
It is easy to see that 9//9 v = + e does not follow from (3) and (4) but 

9f 


== Ye-n (5) 


on the ends of the tube (п — the index of refraction). Therefore the intensity law has 


the form 
Ve Ep*dS = п р^ ао = const. 
2 2 
Наз E 
Ve ^ 
It is easy to check the correctness of (1?) on simple examples. If = = e(x) and the тау 
tube has a direction of the x—axis, 45 = const., by adding formulae (1°) we get 


along a ray tube. (17) 


1 
Ve (ғ E? + Н?) = const. (6) 
€ 

But the expression in the bracket is equal (up to some multiplicative constant) to the 
energy density W of the electromagnetic field. Therefore we have 


ТЕ == = — const. (7) 
€ 


Formula (7) is in accordance with conservation of energy of the electromagnetic field 
since n is inversely proportional to the light velocity in the given medium. From (1) 
it would follow that W = const, which is in contradiction to the Poynting-Umov law. 

We see from (1’) that (for dS = const.) the absolute value of E decreases with 
increasing n and the absolute value of H increases (it would not be so for Friedlander's 
law (1)). Therefore by increasing п the energy of the electromagnetic field passes from 
the electric into the magnetic form. Since only the electric form contributes to the 
light registration (by means of the eye, the photographic plate or other light indicators) 
the light intensity is defined only by the E-vector 


= урус ЕЗ | (8) 


(see Born 1933, p. 30), y is the Poynting vector; we have rejected the unessential con- 
stant с/4л, where c is the light velocity in vacuum. From (17) we have 


‘I 45 = const 


or I = const if dS = const, which expresses the constancy of light flow in the ray tube. 


REFERENCES 
Born, M. Optik, Berlin, 1933. | 
Friedlander, Е. С., Proc. Cambr. Phil. Soc., 43, 284 (1947). 


Vol. XIV (1955) ACTA PHYSICA POLONICA Fasc. 3 


NUMERICAL DETERMINATION OF THE GROUND STATE 
OF THE MESONIC ATOM 


By МАСІЕЈ SUFFCZYNSKI 
Institute of Theoretical Physics, University of Warsaw 
(Received November 29, 1954) 


The distribution of electrical charge inside heavy nuclei influences considerably 
the energy levels of the mesonic atoms. The influence is most marked for the ground 
state and for high values of the atomic charge. Various charge distributions were 
considered: uniform, exponential, and gaussian distributions, as well as the smoothed 
uniform one which according to recent reports is the most probable. The energy 
levels were determined numerically for the greater part of the charge distributions, the 
work of Hill and Ford (1954) being the most exhaustive; see also Fitch and Rainwater 
(1953), as well as Cooper and Henley (1953). 


re 10 20 30-10 "cm 


Fig. 1. Potential energy V(r) for a — 1,744. 1071? cm, that is a — 0,935. 


In May 1954 I began to determine the energy levels of the mesonic atom 
for an exponential charge distribution which seemed promising after the experiments 
on the scattering of fast electrons on heavy nuclei of Hofstadter, Fechter, and McIntyre 
(1953). It is only the simplicity of means by which this determination was done which 
encourages me to publish the details of this work; a preliminary report was already 
published in Nouvo Cimento (1954). 
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The charge distribution. The exponential distribution of atomic charge Ze has 
the form 
о (г) = Ze (2a?) exp ( —r[a), (1.1) 
where r is the distance from the centre, and a is the parameter of the distribution. 
The potential energy resulting from this distribution, obtained by solving the Poisson 
equation AV = ep, has the form 


y(n--ze E — era p : + м] (1.2) 
It has а finite minimum at the origin (see Fig. 1): 
V (0) = — Ze|(2a) i (1.3) 


The equations. Dirac equations are assumed for the fone Their radial parts 
are written, as follows 


dy , К 
= = (k 1) y, (E' Ив (2.1) 
d 
VE (E439 ру (E — V + 2 u dch. 22 - 
' It is convenient to introduce the dimensionless variables 
; : : F 
x = rja DE t3) 
a = alû, where û = (си), — — iu Vee 


| En 


| | eg) 
Equations (2.1) and Q2). for А 5 = Li CR for the | m xi ipo pe be 
, "e Mo - I i. P 1 M LME IE Ro E 
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The boundary conditions. On physical grounds, we have for x = 0 


xy, = 0, xp, = 0. (3.1) 
Assuming in the neighbourhood of x = 0 the developments 
‚ уб) = x (ay + а, x + ag x? + ...), % 5 0), (3.2) 
z (x) = Zo + x! (by + b, x + b x? ...), b, zz 0, (3.3) 
v (x) = + v4 x + vx? + ..., | (3.4) 
4 
j 
2 
Ye 
1 
INN 5 7 9 (ii CRC (727) 
pd 
1 Ж in 
"td 
- Z 
2 
3 ( 
4} 291 % ! + 1 
` BE Wo 9i ) f : E 
| b rd (1% sm) CNE Му апе иан 4 
е“. 5 777. ? 
£ 2. 93 us Герба 2 
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The initial value of z can be chosen arbitrary. For convenience it was taken 


equal to 4 in the present work; thus the initial values were 
у = 0, z= 4, RT 0, z == 0 (3.6) 


The boundary condition at infinity require both y and z to vanish there, because 
the proper y, and y, behave asymptotically like vanishing exponential functions. 

Numerical integration. The point by point integration of eqs. (2.1) and 
(2.2' ) was done for different trial values of E. If such a value was to high both y and z 
were going down too low and z crossed the x axis; for too low trial values y crossed 
the x axis. 

To begin with four initial values of the functions were needed; they were 
found by the standard formulae given by Collatz (1951) on page 43. The formulae 
were iterated four times or more. Usually more than three, sometimes even six, 
decimals recurred in the initial values. 


The initial values being found, Milne's method of integration was used (Milne 


1949, page 135). The formula of Milne 


Ah 0. ; ; 28 4529 
а= zo + эу (да 2 +223) + —d9 — (4.1) 
served as ,predictor and the Simpson rule 
he , А 31?) 
za = zo + چ‎ (ш + 4n + za) — 90” (4.2) 


as „corrector“. For y, which changed more abruptly than 2, no predictor was needed. 
The corrector was repeated when necessary. Weddle's rule was used to check the 
calculations (Milne 1953, page 49). 


The interval h was chosen equal to 1/4 from x = 0 to x = 3 and to 1/2 for greater x. 
Later it turned out that even h = 1 for x > 11 was sufficient to distinguish between 
the energy trial values. The curves for the best trial values crossed the x axis at about 
x — 21. Only two or three decimals recurred then because of the spreading in the 
course of the point by point integration. 


Milne's method was found convenient for use on the calculating machine employed, 
(a Mercedes Euclid „mit Multiplikatorwerk und Speicher“ with 2x8 digits). 


Results. The atomic constants, that is the fine structure constant, electron rest 
mass energy and Compton wave length À,, were taken from Cohen and DuMond (1953). 
The mass of the -meson was taken equal to 207 electron masses after Barkas (Smith, 
Birnbaum, and Barkas 1953). No correction for reduced mass was made, so that the 
unit of length was 


À = Мио) = 4,207 = 1,865 - 10-13 om. 
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The ground state was determined for the Pb atom (Z — 82) for three values of 
the parameter a. The energies found are given in the table below: 


V(O) E' 
MeV MeV 
p 31,646 —10,260 — E' < —10,255 
— 33,846 — 10,609 < Е’ < — 10,602 
— 84,811 —10,757 < E' < —10,750 
10,20 
-1040 


17 18-107 "cm 


Fig. 3. Dependence of the energy eigenvalue E' on the parameter a. 


The functions y and z are shown for a — 0,935 in Fig. 2. It is seen from Fig. 3 that 
by a 1% change of the parameter a the energy of the ground state changes by 0,5%. 
Higher states were not determined. The necessary integrations demanded a greater 
amount of work because of the greater range of the corresponding wave functions. 


I thank Professors J. Werle and M. Danysz for helpful discussions. 
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Taft and Apker (1951, see also Seitz 1954, p. 87) observed the photoelectric 
emission of electrons from F'—centers in RbJ at a temperature of about 85? K; 


the maximum ofthe F-center absorption band for RbJ lay at about 7600 À, correspond- 
ing to a photon energy of about 1,6 eV. Initially no photoelectric emission was obser- 
ved when a sample containing F-centers was irradiated with photons of energy 1,9 


eV (6500 A). After prolonged irradiation, however, the yield curve developed a new 
branch extending to as low energies as 1,6 eV. This branch could be removed by 
heating the sample to room temperature or by irradiating it with 1 eV photons, which 
are known to lie in the F’-band of RbJ. In the present note, the author wishes to men- 
tion some of his own observations which seem to confirm the interpretation of Taft 
and Apker's results as emission from F'-centers. | 

It is well known that the optical absorption band of F-centers in the case of NaCl 


lies in the vicinity of 4700 A (2,6 eV). In this optical region a photoelectric emission 
without preliminary irradiation can be observed when a photon counter is used. 
(Sujak 1953). Irradiation at hy~ 2,6 eV can produce F'— centers besides the above 
mentioned photoelectric emission. The absorption band of F'—centers in NaCl lies 
in the vicinity of 7170 À (1,7 eV). When using a photon counter, the production of 
F'— centers during irradiation with light from the F— center absorption band must 
show itself by an enhanced emission of photoelectrons. In fact, this has been observed 
(see Fig 1) by the author in his experimental arrangement described in 1953. 

On the other hand, the F'—centers must be also detectable by photoemission 
of electrons when irradiating with light of the F' —band region, after preliminary obser- 
vation of the photoelectric emission using light of the F— band region (4700 À — max.). 
This is to be expected, of course, if the temperature of the examined sample (NaCl) 
guarantees the stability in time of the F’—centers. This applies to the room tempera- 
ture of about 18°C in which the author's observations have been made. 

Initially, when irradiated with light of the F’ absorption band the photoelectric 
emission was very low. This small photoelectric emission (which remains rather steady 

d (265) 
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during the time of irradiation) seems to be due to the long-wave branch of the F-band. 
The light used was not strictly monochromatic: 7150 + 600 À and besides the 
total energy of the photons of the F’— band was much larger than the energy corre- 
sponding to the light of the F—band. 

When the NaCl layer containing the F-centers was first. irradiated with photons 
of the F-band! and then with photons of the Е —band (7170 À) an unsteady enhanced 
emission was observed. This enhanced level of photoelectric emission seems to be 


room temp 
F | F 
71704 


OR PRESS 10 15 20 4 
[min]t — TAA 


а à Ad. 


Fig. 1. Time dependence of photoelectric emission of NaCl containing F-centers when irradiated 
with light a to the F and F’ absorption bands 


due to electrons which are ejected from the photoionized F’—centers. Such a pheno- | 
menon was not observed with KCl samples at room temperature, as was to be ‘expected: = 
F’—centers in KC] are not stable at this temperature (Pohl 1938). Further investiges © 
tions of temperature and time de j idence are needed to be able to give a more icc e 
interpretation of the о mechanism of photon counters. nie ^^ 
sathodes are шеш 
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We confine ourselves for the sake of simplicity to the case of nucleon cascade, 
but the results given below may be easily generalized to the case of the electron — 
photon cascade. The development of the nucleon cascade is described by Jánossy's 
G-equation (Jánossy 1950) 


tht 
9 G (e, u, x) егі day DEC, E e у 
66.55 ep ih w (&', =") de аа (Sus) 6 Eux) (1) 
0 0 


where G is the generating function for the probabilities p (n/e, x) to find n nucleons 
of energy >= at the depth x, w(e’, е”) being the cross section for nucleon-nucleon 
collisions normalized to unity. The solution of this equation reads for = > Үл (Lopu- 
szañski 1954) | 


121 


& e | ае de" 
Г„(в, u, x) = (1 — и) $, (e, x) — ЕЕЕ КЕТТІҢ IE 


0 0 


СЕЕ [Еще] = ТЕЙТ " (2) 


Ife, и, x) = 1— С (е, и, x) fore > 


where 


alr 


and S (ex) denotes the first moment of the distribution. 

It is evident that for each threshold energy e>0 the cascade must die out. This 
follows immediately from (1) and (2). 

We introduce now the conditional probability 


i. || 25.63 p (n | £, x) 
DI сл Т шүү САГ 
У’ pml ea) 
m=l 


(265) 
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to find n nucleons of energy > ғ at depth x if at least one nucleon of energy > = 
at depth x exists. It is by no means evident what is the limit of this expression for x 
tending to infinity. We show that for S,*(e, x) = 5, (ғ, Х)/Г (e, 0, x) 
Бот, (Е) = оге б (3) 
xX +00 


The only supposition, we make, is that lim д S,* (e, x)/9x exists. iom (3) it follows 


x00 


than immediately that 
lim p* (Це, x) = 1 


xao 


lim p* (пе, x) =0 п мә 4 


x = 00 


We outline only the proof of (3). Making use of tbe differential equation for 5% (ғ, x), 
we get after some calculations the formula SL 


(а) Zu ее.) | 
ed [зы 5 Tes Les  d[ Іш 


which holds either when 


lim ES (e, x) = const ПЕ ог а 2-0 (6) 
от when 
" = TDI ay 
lim r(z a | dor 0<е «exl MEC 
| == Кең ул | | Шел? T 
16) occurs the rst) folo, еше for E «1 there is alvays Sj (62) =. 


Y- 


Probability Distribution Function of Cosmic Cascades 267 


We have derived this properties of Г from the properties of the first moment 
investigated by Urbanik (1954). Taking ©? into account, we can evaluate the second 
term on the right-hand side of formula (8) and show that 


lim ШЕР — ZA) 
€ E 
€ | de’ ағ” E r(&.0 2 0 


q OPER Г (e, 0, x) F 


The second theorem which we will prove is 


1574 САНАР AM 
HUM Ires! == (0) for = = 2e > 0. (10) 


Indeed (10) is equivalent to (3) as (3) follows from (10) in the same manner as it follows 
from (7) and conversely it may be easily shown by (8) that (3) implies (10). 

The interpretation.of the two theorems (3) and (10) may be stated as follows: 
for large depths of the absorber if we find at least one nucleon of energy above ғ we may 
be sure that there is only one nucleon and its energy lies between e and 2e. In other 
words nucleons of energies above ғ are not able to multiply. 

A more detailed treatment of this problem will be published later. 


The author wishes to express his sincere thanks to Professor J. Rzewuski and 


Mgr K. Urbanik for helpful discussions. 
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We have shown in the previous letter that 


lim p* (Ile, x) = 1, impa л Лх) 02 (1) 
x +00 x00 
and 
ТАМ ЕДІЛ ЕЙ ; 


Formulae (1) give the first approximation of the distribution function for large x. 
To get the second approximation, we must confine ourselves to conditional probabil- 


ities: 
uut UL MN for n 22 
X p (m | е, x) (3) 
m=2 


For further work we need the auxiliary theorem 


1 , 
lim От le; ) E 
CONES GE 


where S, denotes the second factorial moment of the distribution. 
We use also the following equation (Lopuszanski 1954) 


x 1 "Ll 
25 є! D de de" 
г 0 0 


r(5 se) г ose) (5) 
€ € 


(4) 


"n 


| О for ғ > 4e" 


со fore 5 € 
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The function Г (е, и, x) divided by (1—и) gives the generating function for 


у Е, (Eux) 

Zen == m! дит | Ж” 

Dividing the equation for B | Г 
u 


function which tends to zero and making use of (4), we get after some rather compli- 
cated calculations the following result 


| received from (5) by a suitably chosen 
u=0 


l—u 


$2 
5 52 = - == 1. for e>0 (6) 
x p (m | 8, х) 
m= 2 
which means that 
lim p** (21е, x) = 1 and lim p** (nl, х) =0 for n >3 (0 
x +00 с 


In а similar manner we can proceed further in the case of higher approximations. 
We get finally the general result 


HAMM NUT SAGE IS EN o, 
Y eina juae (8) 
From (8) it follows that 
inui. 20 for e > 0 (9) 


x +00 SnlE, x) 


where S, denotes the n-th factorial moment. 
In terms of S,, 1" may be writen as follows 


P(e, u, 4) = (1 a) S (6,4) — (1 — щ# 5 5, (e, а) (1-4) LS, (6,2) — 


(10) 
For large x, (9), the series (10) becomes a Leibniz series, which is very rapidly conver- 
gent. Consequently two approximation steps are already sufficient to obtain a fairly 
good approximation of the probability that at great depth x there is at least one nucleon 
of energy > є: 


$ (в 2) — + Sa (e, а) < T (e, 0,9 < S, (e, a) (1) 


To obtain the expansion of Г іп terms of S, it is convenient to use the recurrent 
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formulae of the solution of the G-equation (Łopuszański 1954) 


x 1 11 
95 є' &'|de' de" 
Tale, и, x) = (1 — и) 5, (e, x) — |а (na 55525 
и) = (L — u) S (ea) ІШІ ne) ff w (5.7: 


0 0 0 


12 
А аш. (5,4 ( 


valid for e > 1/n. Substituting for Г, ,the expressions containing Г, 2 and so on, we 
obtain succesively the terms of expansion (10). 
The computations given above can be easily generalized to the case of the electron 
|. —photon cascade. i an 
A more exhaustive treatment of this problem will be published later. 
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